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ABSTRACT 


A  review  is  given  of  theoretical  studies  on  infrasound  generation  and 
propagation  through  the  atmosphere  which  were  carried  out  under  the  contract. 
These  studies  include  (1)  further  development  and  application  of  n  computer 
program  for  the  prediction  of  pressure  signatures  at  larpe  distances  from 
nuclear  explosions,  (2)  development  of  an  alternative  approximate  model  for 
waveform  synthesis  based  on  Lamb's  edpe  mode,  (3)  development  of  a  peometrical 
acoustics'  theory  Incorporating  nonlinear  effects,  dispersion,  and  wave 
distortion  at  rustics,  and  (4)  a  theoretical  model  for  the  prediction  of 
acoustic  pravity  wave  generation  by  rising  and  oscillating  fireballs. 

Numerical  studies  are  reviewed  which  indicate  the  dependence  of  far  field 
waveforms  on  energy  yield  and  burst  height.  Implications  of  tne  Lamb  edge 
mode  theory  include  a  new  method  for  estimating  energy  yield  from  waveforms 
and  an  explanation  of  amplitude  anomalies  in  terms  of  focusing  or  defocusing 
of  horizontal  ray  paths. 
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Chapter  I 


INTRODUCTION 


1.1  SCOPE  OK  THE  REPORT 

The  present  report  summarizes  investigations  carried  out  hy 
the  authors  during  the  years  1970-1972  on  the  propagation  of  low 
frequency  pressure  disturbances  under  Air  Force  Contract  No. 

F19(>2B-70-C~0008  with  the  Air  Force  Cambridge  Research  Laboratories, 

Bedford,  Massachusetts.  The  study  perform >d  was  theoretical  in 
nature . 

The  central  topic  of  this  study  was  the  generation  and 
propagation  of  infrasonic  waves  in  the  atmosphere,  The  principal 
emphasis  was  on  waves  from  man  made  nuclear  explosions  although 
certain  aspects  of  the  study  pertain  to  waves  generated  by  natural 
phenomena  including,  in  particular,  severe  weather. 

Specific  topics  considered  during  the  study  include  the 
following : 

1. )  The  further  development  and  application  of  a  computer  program 
INFRASONIC  WAVEFORMS  for  the  prediction  of  pressure  signatures  as 
would  be  detected  at  large  horizontal  distances  following  the  detonation 
of  a  nuclear  device  In  the  atmosphere.  The  original  version  of  this 
program  was  developed  by  two  of  the  authors  [Pierce  and  Po3ey,  1970] 
under  a  previous  Air  Force  Contract  (FI 9628-67-C-021 7) . 

2. )  The  development  of  an  alternative  and  considerably  simpler 
theoretical  model  for  waveform  synthesis  based  on  the  hypothesis  that 
dominant  features  of  the  early  arriving  low  frequency  part  of  the 
pressure  wave  are  carried  by  the  real  atmosphere's  counterpart  of 
Lamb's  edge  mode . 

3. )  The  development  of  a  preliminary  theoretical  model  applicable  to 
the  prediction  of  shorter  period  pressure  disturbances  at  intermediate 

to  moderate  distances  from  lower  yield  (1  to  500  kT)  nuclear  explosions. 

4. )  The  source  mechanisms  hy  which  fireball  rise  and  subsequent  near 
field  phenomena  may  excite  acoustic-gravity  waves. 

5. )  The  preliminary  development  of  a  theoretical  model  for  the  generation 
of  acoustic-gravity  waves  by  severe  weather. 

6. )  The  effects  of  atmospheric  turbulence  on  guided  infrasound  propagation. 
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During  the  course  of  the  research  reported  here,  the  study  of 
some  of  the  above  topics  reached  the  point  where  a  definable  phase 
of  the  research  was  sufficiently  complete  that  it  was  logical  to 
summarize  the  results  in  a  form  appropriate  for  publication.  Such 
results  have  been  given  in  seven  published  papers,  in  one  M.S.  thesis 
and  in  one  Ph.D.  thesis.  The  abstracts  of  these  are  listed  in  Chapter  II 
of  the  present  report. 

In  subsequent  chapters,  a  brief  summary  of  the  accomplishments 
under  the  contract  are  given.  These  include  what  the  authors  beliuve 
to  be  the  principal  highlights  of  the  work  already  published  under 
the  contract  plus  a  somewhat  more  detailed  description  of  the  work, 
some  of  which  is  only  partially  complete,  not  previously  reported.  In 
Chapter  VII,  some  recommendations  are  made  for  future  work  in  this  field. 


1.2  BACKGROUND  OF  THE  REPORT 


The  general  topics  of  infrasonic  wave  propagation,  generation, 
and  detection  have  been  of  considerable  interest  to  a  large  segment 
of  the  scientific  and  defense  communities  for  some  time.  A  recently 
published  bibliography  (the  existence  of  which  allows  us  to  omit 
extensive  citations  here)  lists  [Thomas,  Pierce,  Flinn,  and  Craine,  197l] 
over  600  titles,  most  of  which  are  directly  concerned  with  infrasound. 
Literature  pertaining  to  the  infrasonic  detection  of  nucleai  explosions 
constitutes  a  considerable  portion  of  these.  Earlier  work  by  Rayleigh  [l89o]. 
Lamb  [1908,19101,  G.  I.  Taylor  [l 929, 1936] ,  Pekeris  [1939,1938]  and 
Scorer  [19503,  among  others,  which  was  concerned  with  waves  from  the 
Krakatoa  eruption  [Symond,  1888^  and  from  the  great  Siberian  meteorite 
(Whipple,  1930]  is  also  directly  applicable  to  the  understanding  and 
interpretation  of  nuclear  explosion  waves. 


The  present  report  thus  merely  sunaarizes  a  continuation  of  a  small 
number  of  facets  of  a  lengthy  pattern  of  research  wnich  has  been  carried 
on  by  a  large  number  of  investigators  in  the  past.  In  a  more  restricted 
sense,  the  work  reported  here  represents  a  continuation  of  work  done  in 
two  previous  studies  performed  under  contract  for  Air  Force  Cambridge 
Research  Laboratories.  The  first  of  these  was  Air  Force  Contract 
No.  AF19(628)-3891  with  Avco  Corporation  during  1964-1966;  the  second  was 
Air  Force  Contract  No.  AF19628-67-C-0217  with  the  Massachusetts  Institute 
of  Technology  duting  1967-1969.  Summaries  of  this  earlier  work  may  be 
found  in  the  appropriate  final  reports  by  Pierce  and  Moo  [1967]  and  by 
Pierce  and  Posey  [l970]. 

One  of  the  principal  resilts  of  the  two  aforementioned  previous 
contracts  was  a  computer  program  INFRASONIC  WAVEFORMS;  the  deck  listing 
of  the  then  current  version  of  which  is  given  in  the  report  by  Pierce  and 
Posey  [l970] .  This  program  enables  one  to  compute  the  pressure  waveform 
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at  a  distant  point  following  the  detonation  of  a  nuclear  explosion 
in  the  atmosphere.  The  primary  limitation  on  the  program's  applicability 
to  realistic  situations  is  that  the  atmosphere  is  assumed  to  he  perfectly 
stratified.  However,  the  temperature  and  wind  profiles  may  be  arbitrarily 
specified.  The  general  theory  underlying  this  program  is  somewhat  similar 
to  that  developed  by  Harkrider  [l964]  but  differs  from  his  in  that  it 
incorporates  background  winds  and  in  that  it  has  a  different  source  model 
for  a  nuclear  explosion. 
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Chapter  II 


PAPERS  AND  THESES 


2.1  CONTENTS  OF  PAPERS  WRITTEN  DURING  THE  CONTRACT 

A.  P.  Pierce,  J.  W.  Posey,  and  F,.  F.  Iliff,  Variation  of  Nuclear 
Explosion  Generated  Acoustic-Gravity  Wave  Forms  with  Burst  Height 
and  with  Energy  Yield,  J.  Geophys.  Res.  ,  5025-5042  (1971)  . 

A  formulation  of  a  method  for  synthesizing  theoretical  infrasonlc 
far-field  pressure  wave  forms  generated  by  nuclear  explosions  is 
described  which  differr  from  that  described  by  Harkrider  in  1964 
primarily  in  the  method  by  :*Mch  the  source  presence  is  incorporated. 

The  rationale  of  the  model  Is  described  t nd  a  number  of  predicted 
wave  forms  created  by  the  detonation  of  nuclear  explosions  in 
temperature-stratified  and  wind-stratified  atmospheres  are 
presented  to  illustrate  the  effects  of  winds  and  the  effect  of  the 
varying  energy  yield  and  of  the  height  of  hurst.  Theoretical 
wave  forms  are  compared  with  a  numher  of  observed  wave  forms 
previously  exhibited  by  Harkrider.  Discrepancies  are  pointed  out 
anrng  the  predictions  of  Scorer,  of  Weston,  of  Hunt,  Palmer,  and 
Penney,  and  of  Harkrider  as  to  the  effects  of  height  of  burst. 

Our  formulation  predicts  the  early  por  ion  of  wave  forms  received 
on  the  ground  to  he  nearly  proportional  to  yield  and  to  increase 
slowly  with  height  of  burst  up  to  a  height  of  the  order  of  40  km 
(which  depends  on  yield)  and  then  to  decrease  rapidly  with  further 
increasing  height.  These  predictions  are  explained  in  simpler  terms 
if  the  early  portion  of  the  wave  form  is  assumed  to  be  carried  in 
a  single  composite  guided  mode  analogous  to  that  predicted  by  Lamb 
for  the  isothermal  atmosphere. 

A.  D.  Pierce  and  J.  W.  Posey,  Theory  of  the  Excitation  and  Propagation 
of  Lamb's  Atmospheric  Edge  Mode  from  Nuclear  Explosions,  Ceophys.  J.  R.  astr. 
Soc.  26  ,  341-368  (1971)  . 

A  relatively  simple  theoretical  model  of  near-surface  pressure 
pulse  propagation  from  nuclear  explosions  is  developed  from  the 
hypothesis  that  the  early  portion  of  the  pulse  travels  in  the  real 
atmosphere's  counterpart  of  Lamb's  edge  mode.  Various  concepts  of 
geometrical  acoustics  are  used  in  the  construction  of  the  model.  The 
pulse  travels  along  horizontal  ray  paths  which  may  be  refracted  by 
horizontal  variations  in  the  height-averaged  effective  soufcd  speed 
and  wind  velocity.  The  dispersion  of  the  pulse  as  it  propagates 
along  such  paths  is  governed  by  a  one-dimensional  wave  equation 
similar  to  that  derived  by  Korteweg and de  Vries  in  1895.  Coefficients 
are  found  by  comparison  with  the  dispersion  relation  derived  by 
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Garrett  in  1969  as  an  expansion  in  terms  of  a  parameter  c  describing 
deviations  of  the  atmospheric  profile  from  isothermal.  The  method 
of  incorporation  of  terms  governing  accumulative  far-field  non-linear 
effects  is  indicated.  Height  variations  in  propagating  variables  are 
approximated  to  zeroth  order  in  c  .  Variations  along  horizontal 
ray  paths  are  governed  by  the  principle  of  conservation  of  modal 
wave  action.  Hie  excitation  of  the  edge  mode  is  found  hy  matching 
the  general  form  of  the  far-field  quasi-geometrical  solution  to  the 
near-field  solution  for  a  point  energy  source  in  an  isothermal 
atmosphere.  Explicit  expressions  are  given  for  the  far-field  pressure 
disturbance.  Computations  agree  favourably  with  those  based  on  a 
multimode  theory  for  the  first  few  cycles  of  the  waveform.  The  edge 
mode  theory  leads  to  a  number  of  implications,  which  may  be  compared 
with  existing  data.  These  include  the  prediction  that  the  energy 
of  the  explosion  yield  may  be  estimated  fiom  a  knowledge  of  the  first 
peak-to-peak  period  and  first  peak-to-trough  pressure  amplitude  in 
far-field  records  with  very  little  information  concerning  atmospheric- 
structure  ,  and  the  prediction  that  anomalous  azimuthal  variations  in 
waveform  amplitudes  may  be  caused  by  focusing  or  defocuslng  of 
horizontal  ray  paths  for  the  edge  mode. 

J,  E.  Thomas,  A.  D.  Pierce,  E.  A.  Fltnn,  and  L.  B.  Craine,  Bibliography 

on  Infrasonic  Waves,  Geophys.  J.  R.  Astr.  Soc.  2j>,  399-426  (1971)  , 

This  bibliography  is  currei..  through  1970  and  is  composed  of  references 
related  to  the  eener at  ion ,  propsgat lor ,  and  detection  of  sound 
radiation  at  ground  level  and  ionospheric  heights.  The  frequency 
range  of  the  pressure  signals  considered  here  encompasses  those 
normally  associated  with  acoustics,  infraaonic,  acoustic-gravity 
and  gravity  waves. 

C.  A.  Moo  and  A.  0.  Pierce,  Generation  of  Anomalous  Ionospheric  Oscillations 

by  Thunderstorms,  in  Effects  of  Atmospheric  Acoratlc  Gravity  Waves  on 

Electromagnetic  Propagation,  AGARD  Conference  Proceedings  No.  115 

(Harford  House,  London,  1972)  pp.  17-1  to  17-6, 

Recent  experimental  evidence,  based  on  radio  HE  Doppler  sounding 
of  the  ionosphere,  shows  oscillations  of  the  upper  atmosphere 
during  periods  of  thunderstorm  activity.  These  oscillations  have 
periods  in  the  range  of  2  min  to  5  min,  frequently  for  many  hours 
duration.  The  coherence  of  the  oscillations  is  consistent  with 
the  interpretation  generally  given  that  they  are  caused  by  the 
passage  of  long  wavelength  infrasonic  waves.  There  are.  apparently 
no  similar  distinct  oscillations  with  the  same  pe  od  range  associated 
with  air  motion  in  the  troposphere  during  severe  weather.  However, 
convective  activity  is  known  to  generate  fluctuations  and  internal 
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waves  with  periods  near  and  above  Brunt-Vaisala  periods.  The 
present  paper  proposes  a  theory  for  the  generation  of  these 
ionospheric  .2-5  min  period  waves,  based  on  concepts  similar  to 
thoae  used  by  Lighthill  in  the  theory  of  aerodynamic  sound. 

That  the  periods  of  the  ionospheric  oscillations  are  substantially 
lower  than  tropospheric  Brunt -Va’i salt  periods  is  explained  as  being 
due  to  the  fact  that  the  predominant  radiated  frequencies  are 
obtained  as  a  summation,  via  non-linear  processes,  of  the 
frequencies  of  the  interacting  velocity  fields  generated  by  the 
weather  system.  These  radiated  frequencies  are  not  necessarily 
observed  in  the  near  field. 


A.  D.  Pierce,  A  Model  for  Acoustic-Gravity  Wave  Kxcitation  by  Buoyantly 
Rising  and  Oscillating  Air  Masses,  in  Effects  of  Atmospheric  Acoustic 
Gravity  Waves  on  Electromagnetic  Wave  Propagation,  AGARD  Conference 
Proceedings  No.  115  (Harford  House,  London,  1972),  pp.  4-6  to  4-10. 


A  somewhat  general  mathematical  model  is  developed  for  the  study 
of  the  excitation  of  acoustic-gravity  waves  by  rising  and  oscillating 
air  masses.  Sources  are  initially  described  by  distributions  of 
fluid  dynamic  quantities  over  a  moving  closed  surface.  Analysis 
then  indicates  that  insofar  as  wave  generation  is  concerned,  such 
surface  distributions  are  equivalent  to  concentrated  point  sources 
at  the  center  of  the  volume.  The  resulting  linearized  inhomogeneous 
wave  equations  are  derived  and  solved  in  terms  of  Green's  functions. 
The  case  of  an  isothermal  atmosphere  is  discussed  in  some  derail. 


J.  W.  Posey  and  A.  D.  Pierce,  Explosive  Excitation  of  Lamb's  Atmospheric 
Edge  Mode,  in  Effects  of  Atmospheric  Acoustic  Gravity  Waves  on  Electromagnetic 
Wave  Propagation,  AGARD  Conference  Proceedings  No.  115  (Harford  House, 

London,  1972),  pp.  12-1  to  12-8. 


It  has  been  previously  demonstrated  that  far-field  ground-level  pressure 
observations  of  explosively  generated  acoustic-gravity  waves  are  often 
dominated  by  the  Lamb  atmospheric  edge  mode  for  the  first  cycle  or  two. 

In  this  paper,  particular  attention  is  given  to  the  excitation  of  this 
mode  by  a  blast  wave  from  a  large  atmospheric  explosion.  It  is  found 
that  the  strength  of  the  excitation  is  strongly  dependent  upon  the  tail 
of  the  blast  wave.  However,  information  on  the  precise  form  of  the  blast 
wave  tail  is  unavailable,  so  that  any  analytical  representation  must  be 
somewhat  arbitrary. 


A  theoretical  development  shows  that, for  the  pure  Lamb  mode,  a  simple 
analytical  relation  exists  between  the  energy  of  the  source  and  the 
initial  amplitude  and  period  of  the  far-field  pressure  waveform.  This 
relation  is  compared  with  some  empirical  data  and  appears  to  he  in  fair 
agreement  with  yield  estimates  based  on  seismic  observations. 


-19- 


J,  W.  Posey  and  A.  I).  Pierce,  Estimation  cf  Nuclear  Explosion  Energies 
from  Microbarograph  Records,  Nature  232  ,  253  (1971) 

A  theoretical  relation  derived  from  the  edge  mode  theory  developed 
by  the  authors  is  computed  with  microbarograms  published  by  Donn 
and  Shaw  and  by  Harkrider.  Energy  yields  are  taken  as  estimated 
by  Bath.  The  comparison  of  theory  and  data  indicates  the  energy 
yield  may  be  estimated  from  a  knowledge  of  the  first  peak  to  peak 
amplitude  and  from  the  period  between  the  first  two  successive 
positive  peaks. 


2.2  CONTENTS  OF  THESES  WRITTEN  DURING  THE  CONTRACT 

J.  W.  Posey,  Application  of  Lamb  Edge  Mode  Theory  in  Lhe  Analysis  of 
Explosively  Generated  Infrasound,  Ph.D.  thesis.  Department  of  Mechanical 
Engineering,  Massachusetts  Institute  of  Technology,  August  1971  . 

The  propagation  theory  for  Lamb's  atmospheric  edge  mode  of  acoustic- 
gravity  waves  is  presented.  The  excitation  of  this  mode  by  lar  • 
atmospheric  explosiors  is  investigated  and  observations  of  explosively 
generated  infrasound  obtained  from  the  literature  are  examined  for  the 
presence  of  the  edge  mode.  The  theory  utilizes  Brethcrton's  and  Gar¬ 
rett*  3  previous  results  for  the  phase  velocity  and  dispersion  of  this 
mode  and  incorporates  various  techniques  from  geometrical  acoustics 
and  guided  wave  theory.  The  result  is  a  two  dimensional  model  for 
transient  pulse  propagation  over  the  surface  of  a  curved  earth  through 
an  atmosphere  which  is  not  necessarily  perfectly  stratified.  It  is 
found  that  the  excitation  of  acoustic-gravity  waves  by  the  Mast  wave 
from  a  nuclear  explosion  is  primarily  governed  by  the  tail  of  the 
blast  wave.  The  analytical  approximation  introduced  by  Gladstone  is 
selected  to  represent  the  pulse  primarily  because  Of  its  mathematical 
simplicity,  although  its  limitations  are  discussed  in  some  detail.  The 
form  of  the  far  field  edge  mode  pressure  waveform  is  examined,  and  the 
effects  of  some  possible  variations  in  parameters  of  the  medium  and  of 
the  source  are  investigated.  An  approximate  relation  giving  t’’e  source 
yield  (In  equivalent  k.llotons  of  TNT)  in  terms  of  the  first  peak- 

to-trough  amplitude  Pppj>  an^  tfme  between  the  first  two  peaks 

of  the  ground  level  Lamb  mode  pressure  waveform  is 

ykt-o.o3*  ’‘"“W','  T,.23/>m 


which  is  derived  for  bursts  in  the  troposphere  .  Here  p  is  in  Pbars 

FrT 
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and  T  »  in  see.  Comparison  of  the  theory  with  experiment  suggests 

I  ,  z 

that  the  first  cycle  or  two  in  each  of  the  relevant  emriricnl 
microbarograms ,  although  not  necessarily  the  part  of  the  record  with 
the  largest  amplitude,  is  often  dominated  by  the  Lamb  mode. 

L.  Tang,  Propagation  of  Guided  Acoustic  Gravity  Waves  in  a  Turbulent 
Atmosphere,  M.S.  thesis,  Department  of  Mechanical  Engineering,  Mnssachusei  s 
lnstitue  of  Technology,  June,  1971. 

The  effects  of  atmospheric  turbulence  on  the  propagation  of  Lamb 
mode  acoustic  gravity  waves  arc  studied.  A  general  variational 
principle  for  acoustic  gravity  waves  is  first  derived  and  then 
used  to  obtain  the  set  of  approximate  partial  differential  equations 
for  the  propagation  of  guided  I.an.h  mode  waves  in  an  inhomogenous 
atmosphere.  Tbit  set  of  equations  is  then  solved  for  the  scattering 
of  waves  by  turbulence  using  the  Born  approximation.  A  stochastic 
integral  is  subsequently  derived  to  describe  the  random  behavior 
of  the  scattered  wave.  With  the  assumption  that  the  correlation 
function  of  temperature  fluctuations  obeys  the  Kolmogorov  two-thirds 
law,  the  effective  cross  sectional  length  of  scattering  is  obtained. 
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Chapter  III 


THK  MODAI,  SYNTHESIS  OK 
INFRASONIC  WAVEFORMS 


3.1  STATUS  OF  INFRASONIC  WAVK FORMS 

INFRASONIC  WAVEFORMS  is  a  digital  computer  program,  the  first  version 
of  which  was  written  by  Fierce  and  Fosey  [ 1 97 Oj  under  Air  Force  Contract 
AF1 9628-67-C-021 7.  The  nature  of  the  program  is  such  that  it  produces, 
among  other  quantities,  a  synthesized  plot  of  pressure  veisus  time  (Fig.  1) 
which  would  correspond  to  the  infrasonic  waveform  signature  expected  at 
a  given  far  field  point  in  the  atmosphere  following  the  detonation  of  a 
nuclear  device  in  the  atmosphere.  Input  presupposes  the  atmosphere  to  he 
approximated  by  a  sequence  of  horizontal  layers  of  constant  temperature 
and  wind  velocity. 

In  the  time  following  its  original  writing,  INFRASONIC  WAVEFORMS  has 
been  adapted  for  use  at  a  number  of  research  laboratories.  These  include 
the  following: 

1.  Air  Force  Cambridge  Research  Laboratories 
(Mrs.  F,.  F.  Hi f f ) 

2.  National  Oceans  and  Atmospheres  Administration 
(Dr.  J.  M.  Young) 

3.  I.amont-Doherty  Geological  Observatory  of  Columbia  University 
(Dr.  W.  L.  Ponu  and  Dr.  N.  K.  Balachandran) 

4.  University  of  Idaho  (Dr,  J.  F.  Thomas) 

5.  Avco  Everett  Research  Laboratory 
(Mr.  G.  M,  Daniels) 

The  various  versions  of  the  program  in  operation  at  these  laboratories  differ 
primarily  in  the  manner  in  which  the  computer  interfaces  with  the  plotter 
for  the  plotting  of  waveforms. 

In  the  present  section,  we  describe  some  minor  modifications  which 
have  been  made  to  the  program  during  the  past  three  years.  We  also  give 
some  additional  guidelines  to  the  program's  use. 


3.1,1  Earth  Curvature  Effects 

In  the  deck  listing  of  the  subroutine  TMPT  which  is  given  on  pages 
287-290  of  the  report  by  Fierce  and  Fosey  [l970j,we  have  made  one  minor 


Preceding  page  blank 
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Figure  2.  Comparison  of  nuclear-explosion  generatec  acoustic-gravity  waves  upv 
crosswind,  and  downwind  from  the  source.  In  each  case,  the  observer  is  10,000 
from  the  source  in  a  subtropical  surtner  atmosphere.  The  contributing  rodes,  as 
as  their  sum,  are  shown  .for  each  direction. 


change  to  take  into  account  the  curvature  of  the  earth.  On  page  28B 
following  the  statement  NOFN  “  KF1N(N)  which  follows  FORTRAN  st  at  orient  11 
we  have  added  the  following  cards 


C  DETERMINF  THE  EARTH  CURVATURE  CORRECTION  TIMES  ROBS** (-0. 5) 

RAD  -  ROBS  /  6374 . 

CF  -  (  1.  /  (  6374 .  *  SIN(RAD)  )  )**0 .5 
C 

The  program  then  resumes  with  the  COMMENT  card  "THE  MODE  N  RESPONSE  IS 
FOUND...".  Then  FORTRAN  statement  51  is  changed  to  read 


51  TNINT(N,IT)  -  CF  *  TNINT(N.TT) 


such  that  the  factor  (ROBS) ** (-0. 5)  in  the  previous  version  is  replaced  by  CF. 

1  /2 

The  change  just  described  thus  replaces  the  amplitude  factor  (1/r) 
appropriate  to  cylindrical  spreading  by  a  factor  [resin(r/rc) ]~1/^  appropriate 
to  circular  spreading  of  waves  over  a  surface  of  a  sphere  of  radius  r 
Here  r  ■  6374  km  is  the  radius  of  the  earth  and  r  is  the  great  circle 
distance  of  the  observation  point  from  the  source.  The  theory  underlying 
the  incorporation  of  this  earth  curvature  correction  has  been  described 
previously  by  Weston  [l 961 J  and  W  Harkrider  J1964],  among  others.  One  may 
note  that  the  correction  is  minor  if  r  <<  r  .  It  corrects  the  amplitude 
upwards  by  about  IX  at  a  distance  of  5600  kme.  At  the  antipode  where 
r  -  nrp  the  theory  underlying  the  approximation  breaks  down  and  we  have 
not  yet  modified  the  program  to  take  this  case  into  account. 

The  program,  however,  with  relatively  minor  changes  can  still  he 
used  to  estimate  pressure  waveforms  of  waves  which  have  traveled  over  half 
way  around  the  globe,  i.e.  past  the  antipode,  providing  the  wave  has 
traveled  at  least,  say,  1000  km  past  the  antipode.  Following  the  theoretical 
development  of  Rrune,  Nafo,  and  Alsop  [l96l],  one  finds  that  the  appropriate 
modifications  needed  are  (1)  to  replace  sin  (r/re)  by  its  absolute  magnitude: 

(2)  to  shift  the  Fourier  transforms  of  each  modal  waveform  by  a  phase  of  w/2 
radians;  and  (3)  to  reinterpret  r  (ROBS  in  the  program)  as  the  total 
distance  the  wave  has  traveled  (whlcli  in  this  case  would  be  between  it  and 
2nr  ).  These  modifications  can  he  accomplished  by  replacing  the  above  definition 


of  CF  in  subroutine  TMPT  by 

CF  -  (  1.  /  (  6374.  *  ARS(SIN(RAD  ))))**0.5 


Then,  also  in  the  deck  listing  of  ITfPT,  between  lines  TMPT  149  and  TMPT  150 
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as  given  on  page  289  of  the  report  by  Pierce  ami  Posey  fl970],  one 
should  insert  a  card 

PH 2  -  PH 2  +  1.570796 


An  identical  card  should  also  be  Inc-irted  between  lines  TMPT  166  and 
TMl'T  167.  The  input  value  of  ROBS  should  also  be  selected  to  correspond 
to  the  longer  propagation  path. 

One  should  note  that,  in  regard  to  the  above  suggested  modifications 
fer  antipodal  arrivals,  the  underlying  theoretical  basis  presumes  the 
atmosphere  to  he  cyl indri cal ly  symmetric  about  an  axis  passing  through 
the  point  of  detonation.  Since  this  is  a  rather  stringent  idealization, 
considerable  caution  should  be  used  in  the  interpretation  of  the  results 
of  the  computations. 


3.1.2  Other  Modifications  in  the  Program 

Also,  in  regards  to  subroutine  TMPT,  a  minor  change  in  the  FORMAT 
statement  104  was  to  change  the  FORMAT  specification  F12.2  to  F12.4  to 
allow  more  significant  figures  to  be  printed  in  the  output  tabulation  of 
TNINT  (TOPT.TT). 

One  minor  error  in  the  deck  listing  of  the  subroutine  NXTPNT  was 
pointed  out  to  the  authors  by  T.  H.  Kuckertz  of  the  University  of  Idaho.  In 
FORTRAN  statement  26  of  this  subroutine,  the  nurabei  -5  should  have  been 
punched  +5.  As  best  we  can  tell,  this  error  will  never  affect  the  output 
of  the  program  since,  for  any  conceivable  atmosphere  model,  IPK  can  never 
be  5  at  a  point  when  statement  26  is  reached. 

In  subroutine  SOURCE,  the  present  version  at  M.l.T.  has  set  TAS  -  0.33 
instead  of  0.48  as  was  done  in  the  original  listing.  The  revised  source 
model  theory  described  here  in  Sec.  3.2  and  by  Pierce,  Posey,  and  II iff  fl971J 
suggests  that  this  should  be  the  appropriate  choice  if  the  source  is 
sufficiently  high  that  there  is  no  appreciable  nonlinear  interaction  between 
the  incident  shock  and  the  ground  reflected  shock. 

If  the  explosion  is  relatively  close  to  the  ground,  the  subroutine 
SOURCE  should  be  modified  following  a  suggestion  made  to  the  authors  by 
J.  Reed.  In  accordance  with  the  theory  outlined  in  the  next  section,  one 
should  replace  the  statement  computing  the  quantity  PAS  in  the  subroutine 
by 


PAS  -  ( (21 . 4E+3) /I • 0) *(1 . 61 ) 
and  the  statement  defining  TAS  by 
TAS  »  0.416 
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Just,  when  a  Riven  explosion  should  he  classified  as  near  ground 
or  above  ground  depends  primarily  on  the  ratio  of  the  height  of  hurst  h 
to  the  explosion  characteristic  length  ( K/ p0 )  3 / 3  where  F  is  the  energy 
yield  and  p0  is  ttie  atmospheric  pressure.  If  the  ratio  is  (treater  than, 
say  5,  the  above  Rround  source  model  would  appear  to  he  more  appropriate. 

If  it  is  less  t ban  1,  the  near  Rround  source  model  would  appear  appropriate 
For  intermediate  ranges  of  altitude,  lust  which  choice  is  hotter  remains 
a  topic  of  investiRation.  (Dr.  !i .  I,.  Murphy  and  his  colleagues  at  Mt  . 
Auburr  Research.  Associates  arc  currently  developing  an  improved  source 
model  which  takes  the  nonlinear  interaction  between  incident  and  reflected 
shocks  into  account.) 

Wc  would  suRRest  lor  the  present  that  one  take 
‘  5™  <'KT>'/3ft 


as  defining,  a  surface  burst,  where  Y„_  is  the  yield  in  KT. 

K  i 

3.1.3  Further  SuRRostions  for  Program  Users 

In  the  usage  of  the  propram  hy  oursel»PN  and  by  our  co] leagues  In 
other  laboratories  wc  have  become  aware  of  a  number  of  possible  pitfalls 
in  its  use.  In  order  to  increase  the  effectiveness  of  future  use  of  the 
program  it  would  seem  appropriate  to  mention  these  here.  The  present 
discussion  amends  that  of  Chapter  111  in  the  previous  report  by  Fierce 
and  Tosoy  ^1 970 J . 

While,  apriori,  one  should  not  expect  the  temperature  In  the  uppermost 
layer  (a  halt'spacc  in  the  theory)  to  affect  waveforms  observed  near  the 
ground,  the  actual  computations  may  indicate  otherwise.  This  would  appear 
to  he  due  to  the  neglect  of  leaking  modes  in  the  theoretical  formulation. 
Any  mode  or  portion  of  a  mode  which  1< aks  energy  into  the  upper  halfspace 
is  associated  with  a  complex  wavenumber  and  appears  formally  to  tie  damped 
with  propagation  in  the  tiorizont.nl  direction.  The  possible  values  of 
phase  velocity  vp  and  angular  frequency  ui  which  could  in  principle 
correspond  to  a  nonleaking  mode  is  limited  primarily  by  the  choice  of 
temperature  in  the  uppermost  halfspace.  The  forbidden  regions  where  all 
conceivable  modes  must  he  leaky  correspond  to  values  of  w  and  Vp  for 
which  ttie  normal  mode  dispersion  function  does  not  exist.  Such  domains 
in  ttie  w  vs  v  plane  correspond  to  the  regions  containing  X's  in 
computer  generated  plots, of  normal  mode  dispersion  waves  such  as  shown  in 
Fig".  3,3  and  3.4  of  the  previous  report  hy  Fierce  and  Posey.  Interestingly 
if  one  takes  the  atmosphere  to  he  of  some  very  large  hut  finite  height  li 
terminated  by  either  a  rigid  surface  or  a  free  surface,  the  modes  are  all 
non  leaky  (there's  nowhere  to  leak)  and  the  dispersion  curves  and  height 


profiles  for  inodes  which  propagate  primarily  in  the  lower  atmosphere 
are  nearly  the  same  in  both  cases.  If,  however,  one  terminates  the 
atmosphere  by  an  isothermal  halfspace  extending  above  h,  certain  portions 
of  the  lower  atmosphere  modes  will  become  leaky.  Since  the  computer  program 
discards  such  modes,  regardless  of  how  small  the  leakage  may  actually  be, 
it  is  difficult  to  say  whether  more  nearly  correct  results  will  be  obtained 
with  the  unbounded  atmosphere  model  than  with  one  terminated  above  by  either 
a  free  surface  or  a  half  space. 

A  specific  suggestion  for  circumventing  this  problem  without  actually 
revising  the  computational  procedure  to  include  leaky  modes  is  that  one 
always  choose  the  temperature  of  the  upper  halfspace  to  be  sufficiently  high 
that  oig  "Or-Dl^g/c  for  the  upper  halfspace  is  smaller  than  any 
frequency  one  suspects  should  be  dominant  in  the  synthesized  waveform.  Also, 
one  should  want  the  sound  speed  in  the  upper  half  space  to  exceed  the 
horizontal  phase  velocity  of  any  wavelet  which  contributes  substantially 
to  the  overall  waveform.  We  suggest  that  one  check  his  output  after  the 
computation  to  see  if  his  choice  of  a  model  atmosphere  satisfies  these 
criteria.  If  not,  he  should  either  revise  his  model  to  accomplish  this, 
subject  to  whatever  constraints  one  imposes  on  the  relation  of  the  model 
to  reality,  or  else  he  should  be  very  careful  in  the  interpretation  of 
hxS  computational  results. 

Other  problems  which  arise  in  the  use  of  the  program  are  associated 
with  the  truncation  of  the  integration  over  frequency  in  the  synthesis  of 
individual  nodal  waveforms.  When  one  chooses  his  input  parameters  0M1 . 

OM2,  VI,  and  V2,  he  implicitly  limits  the  frequency  range  for  the  Fourier 
inversions  to  an  interval  contained  between  0M1  and  OM2  radians.  The 
interval  Is  also  limited  by  the  fact  that  the  dispersion  ann/e  fo  the  non- 
leaking  portion  of  any  given  mode  may  not  exist  for  all  frequencies  between 
0M1  and  0M2  and  also  by  the  fact  that  its  phase  velocities  may  not  fall 
between  VI  and  V2  km/sec.  Thus  the  integration  represents  an  approximation 
and  one  should  make  some  assessment  as  to  whether  the  app  oximation  has 
any  validity. 

One  possible  simple  means  of  checking  whether  the  truncation 

approximation  is  applicable  is  to  examine  the  range  of  group  velocities 

printed  out  by  the  program  for  the  Individual  inodes.  Suppose  the  group 

velocities  vary  in  magnitude  between  v  ,  and  v  for  an 

8, min  g,max 

individual  mode.  Then  one  would  expect  the  contribution  from  that  mode 

to  be  of  significant  magnitude  at  best  only  for  times  between  r/V  and 

g  ,max 

r/v^  with  a  possible  extension  of  perhaps  one  or  two  cycles  on  either 

side  of  this  interval.  If  this  is  not  the  case,  then  the  approximation  Is 
probably  insufficient.  If  one  suspects  it  to  be  insufficient,  then  h<  should  next 
check  to  see  If  the  mode  in  question  gives  a  substantial  contribution  to 
the  overall  waveform  for  the  time  interval  of  Interest.  If  It  doesn't,  then 
one  probably  need  not  worry  about  the  error  introduced  by  truncation  in  as  far  as 
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that  mode's  contribution  is  concerned.  If  it  does,  then  one  should  alter 
his  input  parameters  to  improve  the  situation  either  by  extending  the  range 
of  frequencies  and  phase  velocities  considered  or  else  hv  changing  the 
assumed  form  of  the  upper  atmosphere  to  make  the  mode  noi  ‘akin,,  over  a 
wider  frequency  interval. 

Anyone  who  has  attempted  to  use  the  program  to  synthesize  waveforms 
at  relatively  close  in  distances  may  have  noted  that  the  individual  modal 
waveforms  may  not  be  causal,  i.e. ,  the  program  may  formally  predict  a 
disturbance  before  the  time  of  detonation,  which,  of  course,  is  meaningless. 
In  principle,  the  superposition  of  all  the  modes  should  result  in  a  causal 
waveform, hut  there  is  no  apriori  reason  in  the  mathematical  theory  why 
any  given  individual  waveform  should  so  be.  At  larger  distances,  this  lack 
of  modal  causality  is  not  a  significant  factor  since  the  noncausal  portions 
of  the  synthesized  waveform  tend  to  be  of  negligible  amplitude  compared 
to  the  portions  which  arrive  at  times  corresponding  to  the  group  velocity. 

Another  question  one  should  consider  is  whether  one  has  incorporated 
enough  modes  in  his  synthesis.  If  one  is  seeking  to  descrihe  the  disturbance 
in  a  given  time  interval  at  range  r  then  he  probably  need  only  consider 
those  modes  and  frequency  intervals  which  correspond  to  group  velocities 
r/t  where  t  falls  within  the  desired  time  interval.  In  general,  modal 
contributions  tend  to  decrease  with  increasing  mode  number,  so  if,  for 
example,  one  has  convinced  himself  that  contributions  from  the  ,  S  ^ 

modes  are  all  negligible,  then  the  same  should  be  true  foi  the  S^,  Sg ,  etc. 


3.2  SOURCE  MODELS 

The  source  model  for  TNFRASONIC  MODES  has  been  somewhat  modified  since 
the  writing  of  the  report  by  Pierce  and  Posey  [jl970] .  Here  we  outline  the 
revised  theory  for  the  present  source  model.  The  first  part  of  the  discussion 
follows  that  outlined  in  the  paper  by  Pierce,  Posey  and  Iliff  [l97l]. 


3.2.1  Above  Oround  Explosions 

The  basic  mathematical  model  is  the  linearized  equations  of  hydrodynamics 
for  a  temperature-stratified  and  wind-stratified  compressible  atmosphere.  The 
geometry  is  sketched  in  Fig. 2  .  The  presence  of  the  source  is  assumed  to 
be  taken  into  account  by  the  addition  of  a  source  term  to  that  equation 
corresponding  to  energy  conserva' 1  on .  Thus,  the  governing  equations  for 
acoustic  pressure  deviation  p  i"  om  ambient  pressure  pc,  acoustic  density 
deviation  p,  and  fluid  velocity  deviation  u,  are  taken  as 

pofDtu  +  -  -  9p  -  RPeup  (3.1a) 
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Figure  2.  Sketch  illustrating  the  various  idealizations  incorporated 
in  the  mathematical  formulation  of  the  problem  of  acoustic-gravity  wave 
propagation  from  a  nuclear  blast. 
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(3.1h) 


D  P  +  V.f()  u)  *0 
t  o 

2  «♦ 

[D  p  +  u*Vp  ]  -  c  [D  P  +  u-Vp  ] 

1  t  o  t  a 

*  4nc?f  (t)  i(?  ~  r"  )  (3. ic) 

E  o 

where 

D  -  (3/3t)  +  v-7 


is  the  time  derivative  corresponding  to  an  observer  moving  with  the 
ambient  wind  velocity  v.  The  quantity  “e  is  the  unit  v  ctor  in  the 

vertical  direction,  while  r0  denotes  the  source  location;  c  «=  £yPo7pq] 

is  the  speed  of  sound.  The  source  term  (and,  in  particular,  the  form  of  the 
function  fE(t)  on  the  right-hand  side  of  Eq .  (3.1c))  is  chosen  with 

the  objective  of  insuring  that,  at  larger  distances,  the  solution  of  Eqs. 

(3.1a,  b,  c) should  agree  with  what  might  be  expected  if  the  nonlinear  terms 
were  retained  arid  a  more  realistic  source  model  were  used. 

To  explain  our  choice  for  the  energy  source  function  f  (t),  the  basic 

t 

nonlinear  hydrodynamic  model  for  a  nuclear  explosion  [Taylor,  1950;  Brode, 

1955,  1968]  is  reviewed,  which  consists  of  an  initially  isothermal  sphere 
of  very  small  radius  in  an  unbounded  isothermal  atmosphere  with  negligible 
gravity.  The  initial  sphere  has  ambient  density  and  fluid  velocity,  but  is 
assumed  to  have  very  high  temperature  and  pressure.  The  total  initial 
potential  energy  (the  specific  heat  of  air  is  assumed  independent  of 

temperature)  inside  the  sphere  is  assumed  to  be  some  fixed  fraction  of  the 

total  energy  Y  (in  kilotons  of  TNT,  where  1  KT  *  4,2  *  10^  ergs) 

released  by  the  explosion.  According  to  The  Effects  of  Nuclear  Weapons  (ENW) 
[Glasstone,  1957,  1962],  this  initial  potential  energy  (which  should  be 
the  3ame  as  the  total  hydrodynamic  energy)  is  roughly  one-half  of  Y^  .  We 

assume  this  proportionality  factor  does  not  vary  with  yield  or  with  height 

of  burst.  It  would  appear  from  comments  in  sections  2.119,  7.25,  and  from 
Figure  1.22  in  ENW [Glasstone ,  1962]  that  this  is  a  good  approximation  for 
bursts  below  100,000  feet  (30  km)  and  a  fair  approximation  at  altitudes  up 
to  350,000  feet  (100  km),  although  no  quantitative  description  of  the 
variation  of  this  fraction  is  given.  (In  any  event,  the  results  of  thi3 
study  will  not  be  invalidated  by  such  a  variation,  provided  energy  yield 
is  always  interpreted  as  twice  the  hydrodynamic  energy  of  the  explosion.) 

The  pressure  wave  form  in  this  Idealized  basic  nonlinear  model  corresponds 
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to  hydrodynamic  scaling,  i.e.,  at  distance  R 


p  -  poF(R/cTy,  t/Ty) 


0.2) 


where  F  is  a  universal  function  of  two  dimensionless  quantities  and 

where  Ty  is  a  characteristic  time  that  may  he  defined  as  any  convenient 

constant  times  (Y„~/p  )*'3/c.  We  choose  T„  to  be  defined  such  that. 

Ki  o  i 

when  pq  and  c  correspond  to  a  standard  reference  pressure  (p  ^  “  1  atm) 

and  sound  speed  (cre£  “  331  m/sec)  and  when  is  1  KT,  then  Ty  is  eq  .1 

to  the  time  duration  t  ,  of  the  blast  wave  at  a  distance  R  (which  we 

ref  ref 

take  as  1  mile)  from  a  1-KT  explosion.  (As  explained  later,  the  appropriate 
numerical  choice  of  t  ^  would  appear  to  be  0.33  sec.)  Thus, 


(where,  to  maintain  dimensional  purity,  Y  is  here  interpreted  as  the  ratio 
of  the  explosion's  energy  to  a  reference  energy  of  4.2X1019  ergs). 

Since,  in  the  homogeneous  atmosphere  case  currently  being  considered, 
the  amplitude  of  F  may  be  expected  to  fall  off  nearly  inversely  with  R 
(spherical  spreading  in  linear  acoustics),  Eq. (3. 2) should  be  rewritten 

in  the  form  ' 


P  -  <aP)ref(P0/Pref)^TY/(creftref)](Rref/R) 


X*([t  -  R/c]/Ty,R/(cTy>)  (3.4) 

where,  at  large  distances  R,  ii  should  be  a  relatively  slowly  varying  function 
of  its  second  argument.  Here  (AP)ref  is  the  shock  overpressure  at  the 

reference  distance  R  ^(1  mile)  from  a  1-KT  air  burst  in  a  standard  reference 

atmosphere.  (Ap)  ^  is  taken  to  be  34  mbar;  the  function  is  dimensionless 

and  is  normalized  such  that  it  has  a  value  of  unity  immediately  behind  the 
shock. 
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With  reference  to  the  mnlinear  model  described  above,  a  suitable 

choice  for  f_,(t)  woulo  appear  to  be  that  furction  which  would  lead  to 
h 

a  solution  of  Eqs.(3.1a,  b,  c)  which  agrees  with  Eq-  (3.4)  at 

moderate  distances  from  the  source  when  gravity  is  neglected  and  when 

the  atmosphere  is  assumed  to  be  homogeneous.  The  ambient  atmosphere  for 

this  matching  process  is  taken  as  that  characteristic  of  the  burst  location. 

The  distance  selected  for  the  matching  is  {cT„/fc  „t  ,1}  (i.e.,  the  scaled 

Y  ref  ref  ’ 

counterpart  of  1  mile  from  1  KT  in  the  reference  atmosphere),  such  that 
cTy/R  is  identically  equal  to  cref cref ^Rref  • 

Since  the  linearized  equations,  vlth  the  neglect  of  gravity  and 
inhomogei  jities,  give 

p  -  R-1fE' (t  -  R/c)  (3.5) 

the  matching  of  the  two  solutions  under  the  circumstances  described 
above  gives 

fE(t)  -  BoPocTy2Q(t/TY)  (3.6) 

where  the  dimensionless  constant  is  given  by 

Bo  *  f(Ap)s/prefJ{Rref/(Creft„f]}  (3.7) 


(Its  value  for  the  choice  of  numeric*-!  values  described  previously  is 
0.40.)  The  function  Q(t/Ty)  is  simply  the  integral  of  ij>,  i.e., 


Q(0 


W.,  A  )  dC 

o 


(3.8) 


where  the  integration  constant  is  such  that  fg(t)  vanishes  at  sufficiently 
early  times.  The  quantity  Aq  is  just  [R  ^ /c  j. tref  J  (which  turns 
out  to  have  a  numerical  value  of  14). 

The  apparent  arbitrariness  in  the  selection  of  the  distance  R  is 

largely  immaterial  insofar  as  the  prediction  of  the  far-field  wave  form 
is  concerned.  It  turns  out  that  the  low-frequency  pulse  predicted  at  large 
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distances  (as  well  sight  be  expected  for  any  low-frequency  pulse)  la 

largely  governed  by  the  first  few  terns  in  the  power  series  expansion  In 

angular  frequency  w  of  the  Fourier  transforn  (over  tine)  of  f„(t)  . 

£. 

Thus,  a  hey  parameter  (which  may  be  considered  as  proportional  to  the 
ceroth  order  tarn  In  this  expansion)  would  be 


4*  |  fE(t)  dt  -  (y  -  1)(K/c2)E  (3.9) 

19 

where  E  -  4.2  X  10  Y^  ergs  is  the  total  energy  released  by  the  explosion 
and  where  K  is  a  dimensionless  constant  given  by 

K  -  4ir (y  -  l)"1(Ap)refRref(trefCref)2 

xj  Q(£)  d£/(4.2Xlo’9  ergs)  (3.10) 


It  is  assumed  here  that  fj.(t)  1«  integrable,  and,  in  particular,  that 

it  vanishes  as  t  «•  .  This  would  seem  to  be  a  necessary  requirement 
since  (at  leant  In  the  acoustic  limit  implied  by  Eq -  (3.5))  the  particle 

displacement  is  proportional  to  fg/R.  If  f£(t)  were  to  approach  a  non¬ 
zero  limit  as  t  this  would  tend  to  imply  that  the  net  amount  of  work 

done  in  expanding  or  compressing  the  gas  behind  the  shock  diverges  (rather 
than  approaches  a  limit)  as  the  shock  radius  Increases  without  limit,  which 
in  turn  implies  that  the  energy  associated  with  the  source  1b  unlimited. 

The  actual  near-field  blast  form  can  be  inferred  from  Brode's  [1955] 
calculations  with  the  hydrodynamic  energy  (which  Brode  denotes  as  E  ) 

interpreted  as  one-half  the  total  energy  released  by  the  explosion.  Brode's 
calculations  show  the  blast  beginning  with  a  shock  whose  overpressure 
decreases  faster  than  1/R,  followed  by  a  positive  phase  whose  duration  is 
nearly  constant  with  R  and  equal  to  1.22  <Etot/PCJ>1/3/c  •  Although 

very  little  information  can  be  inferred  from  Brode's  paper  concerning  the 
wave  form  after  the  negative  phase,  it  seems  (from  extrapolation  of  the 
calculated  wave  form)  that  there  is  a  very  3raall  oscillation  whose  period  is 
roughly  of  the  order  of  twice  the  negative  phase  and  whose  amplitude  decreases 
fairly  rapidly.  The  energy  feeding  this  oscillation  is  the  energy  dissipated 
at  the  Btart  of  the  positive  phase.  On  this  basis,  on«  can  argue  that  the 
value  of  the  constant  K,  as  computed  from  Eq. (3.10)  should  be  independent 

of  the  choice  of  Rref»  providing  that  t  j-,  and  Q(£)  correspond 
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to  whatever  choice  is  made.  The  rationale  cf  this  conclusion  is  that 
the  internal  energy  relative  to  ambient  is  p/('t  -  1)  per  unit  volume; 

if  p  is  taken  as  given  by  Eq  .(3.4),  and  if  R/cT^  in  the  second 

argument  of  is  approximated  by  Rref^c^y^*  and  if  the  total  internal 

energy  (  a  spatial  integral)  is  evaluated  at  the  instant  the  shock  radius 
is  the  scaled  counterpart  of  R  one  finds,  after  an  integration  by 

parts,  that  this  energy  at  large  t  approaches  (K) <E)  where  K  is 
defined  in  Eq.  (  3.10).  The  fact  that  the  major  contribution  to  the 
integral  comes  at  larger  R  Justifies  one  setting  R  ■  Rref  in  the 

second  argument  of  . 

Since  we  would  expect  this  limit  to  exist  (even  the  dissipated  energy 
must  eventually  be  transformed  Into  internal  energy)  we  must  conclude  that 
K  is  independent  of  R  ^  at  sufficiently  large  Rref*  Another  conclusion 

which  seems  warranted  is  that  K  should  not  exceed  one-half  since  the  increase 
in  internal  energy  cannot  exceed  the  total  hydrodynamic  energy  released  by 
the  explosion. 

Concerning  the  choices  of  Rref>  tref*  (AP)ref»  we  taRe  Rref  "  *  wile, 

since  this  seems  amply  large  for  a  1-KT  explosion,  the  overpressure  being  of 
the  order  of  3Z  ambient  at  this  distance.  The  actual  value  of  (Ap)^^  can 

then  be  determined  from  Brode's  calculations  to  be  34  mbar.  The  same  number 

car.  also  be  inferred  from  Figure  3.946  m  the  1957  edition  of  ENW  [Classtcne, 

l°r*J.  (The  figure  gives  the  overpressure  at  the  ground  at  a  distance  of 

from  a  1-KT  air  burst  as  1.00  psi.  Since  nonlinear  effects  are  small 

a.  .is  distance,  we  expect  the  overpressure.  If  the  presence  of  the  ground 

we  a.gl.ected,  to  be  one-half  of  this  or  0.50  psi  *  34  mbar.)  The  choice 

of  -  e£  as  inferred  from  Brode's  results  (with  an  extrapolation  based  on 

the  tfieory  of  weak  spherically  diverging  shocks  [Bethe  and  Fuchs,  1947; 

Landau  1945;  Whitham,  1950])  is  found  to  be  0.33  sec.  On  the  other  hand, 

Figure  i  96  in  the  1957  edition  of  ENW,  with  an  appropriate  extrapolation, 

gives  f  ef  “  0.48  sec.  Since  the  latter  is  not  documented  and  hardly 

cons.  ifes  an  adequately  reported  piece  of  research,  we  prefer  to  adopt 

the  number  based  on  Brode's  computations  of  t  .  *  0, 33  sec.  In  addition,  as 

rci 

discussed  below,  energy  considerations  would  tend  to  suggest  that  the  ENW 
figure  is  inappropriate. 

With  regard  to  the  functional  form  of  ^  ,  it  would  appear  adequate  to 
use  a  fairly  simple  analytic  approximation,  represented  by  the  function 
[Giasstone,  1962] 


*(C,  Ao)  -  (1  -  F,)e~KV(0 


(3.11) 


when  U(0  Is  the  Heaviside  step  function  and  F, 


is  measured  from  the  onset 
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of  the  shock.  While  this  function  has  the  disadvantages  of  giving  more 
'sag'  to  the  positive  phase  decay  than  implied  by  Brode's  fit  as  represented 
by  his  equations  22  and  23  and  of  hiving  a  negative  phase  of  unlimited 
duration,  it  does  have  the  desired  features  of  leading  to  an  fE(t),  which 

vanishes  at  sufficiently  large  t,  and  of  giving  a  finite  value  to  K  in 

Eq-  (3. 10),  although  it  is  not  clear  to  what  extent  its  predicted  Integral 
over  Q,  which  turns  out  to  be  1 ,  would  agree  with  that  for  the  true  wave 
form.  Considering  the  fact  that  the  integral  simply  denotes  the  center 
of  gravity  of  the  wave  form  iJ/(£  ,  A  )  with  respect  to  the  start  of  the 
blast  wave,  one  would  doubt  that  th?s  center  of  gravity  would  fall  too  far 
from  the  location  of  the  end  of  the  positive  phase  (which  is  st  £  •  1). 

One  ancillary  reason  for  adopting  Eq  •  (3.11)  is  that  the  same  function 

has  been  used  previously  by  Harkrider,  and  we  wish  to  compare  our  results 
with  his. 

With  t|/  given  by  Eq.  (3.11),  the  source  function  becomes 


f£(t)  -  (4k)-1 (y  -  l)(K)(E/c2) 

•(t/TY2)e“t/TYu(t)  (3.12) 

where  (with  tref  -  0.33  sec)  the  value  of  K  in  Eq.  0.10)  is  taken 

an  0.50.  In  contrast,  if  we  used  t  *  ■  0.46  sec,  we  would  have  found 

tre 1 

K  -  1.1.  Since,  as  pointed  out  previously,  there  are  strong  theoretical 
grounds  for  supposing  K  to  be  less  than  one-half,  the  choice  of 
t  f  -  0.33  sec  based  on  Brode's  calculations  appears  preferable.  The  time 
off gin  in  Eq.  (3.12)  was  chosen  so  that  the  shock  would  arrive  at  the 

matching  distance  at  a  time  equal  to  this  distance  divided  by  the  ambient 
sound  speed.  While  this  does  not  necessarily  coincide  with  the  time  of 
detonation  (due  to  finite  amplitude  effects),  the  choice  of  origin  is 
Immaterial  insofar  as  the  amplitude  and  relative  shape  of  the  received 
wave  form  are  concerned. 

A  shoitcowing  of  the  source  model  described  above  is  its  neglect  of 
gravitational  effects  and  of  sound  speed  and  density  gradients  during  the 
early  history  of  the  blast.  Apart  from  the  desire  for  a  relatively  simple 
model,  the  justification  for  this  neglect  is  that  the  influence  of  these 
effects  accumulates  at  a  relatively  slow  rate  and  does  not  become  appreciable 
until  relatively  late  during  the  history  of  the  formation  of  the  blast  wave. 
Since  we  are  interested  primarily  in  the  earliest  portion  (say,  the  first 
15  to  30  minutes)  of  the  pressure  pulse  which  propagates  through  the  lower 
atmosphere  (which  is  frequently  all  that  ie  detected  by  ground-level 
tmterobarographs) ,  we  may  view  the  early  portion  of  the  blast  wave  (including 
the  first  positive  and  negative  phases)  as  the  source  of  the  pulse. 
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It  would  appear  that  waves  Renerated  by  fireball  rise  due  to  buoyancy 
may  be  excluded  on  the  basis  that  they  would  arrive  too  late. 

From  the  discussion  accompanying  Eq.  (3.9^  it  would  also  appear 
that,  insofar  as  the  low-frequency  far-field  pulse  is  concerned,  the  only 
critical  assumptions  are  (i)  that  the  enerRy  radiated  by  the  near-field 
blast  be  radiated  isotropically  and  (ii)  that  the  total  additional  internal 
energy  in  the  atmosphere  caused  by  the  explosion  be  of  the  order  of  one- 
half  the  total  energy  yield.  It  would  appear  that  these  assumptions  would 
be  valid  provided  the  format  Ion  of  the  positive  phase  (which  spreads 
isotropically  and  which  includes  the  bulk  of  the  hydrodynamic  energy  at 
very  early  times)  should  not  be  altered  significantly  (in  terms  of  energy 
content)  by  gravity  and  inhomogeneities.  A  possible  quantitative  criterion 
is  that  the  time  for  the  positive  phase  to  form  (which  may  be  taken  as  being 
of  the  order  T^)  be  less  than  one-fourth  (since  the  positive  phase  resembles 
one-fourth  of  a  full  cycle  of  a  sine  wave)  of  the  characteristic  period  for 
acoustic-gravity  waves  in  the  atmosphere.  This  period  [Vals'ala,  1925; 

Brunt,  1927]  is  generally  greater  than  5  minutes.  Thus,  we  could  take 
T  <  1  minute,  as  a  possible  restriction.  This  would  restrict  the  model  to 
•  ^ 

yields  of  less  than  about  10  MT  at  sea  level  and  to  yields  of  less  than 
1  MT  at  65  km.  We  should  emphasize  that  the  prominence  of  periods  of  the 
order  of  the  brunt-Vaisala  period  in  the  far-field  wave  forms  does  not 
necessarily  indicate  the  inappropriateness  of  the  model,  as  these  periods 
can  arise  through  the  accumulative  lonR-range  dispersion  of  the  early 
portion  of  the  blast  wave  form. 


3.2.2  Near  Surface  Explosions 

The  following  describes  the  appropriate  source  model  for  near  surface 
explosions  which  was  developed  following  a  suggestion  by  J.  Reed  of  Sandia 
Corporation.  Suppose  Y  is  the  total  energy  yield  of  a  near  ground 
explosion.  Then  the  blast  wave  at  moderate  distances  might  appear  (from 
symmetry  considerations)  to  be  the  same  as  from  an  air  burst  of  yield  2Y 
when  the  presence  of  the  ground  is  neglected.  On  the  basis  of  this  assertion 
and  the  reasoning  outlined  in  the  previous  section,  we  would  infer  that 
the  blast  overpressure  (Ap)*^  and  positive  phase  duration  tref  of 

the  blast  wave  waveform  on  the  ground  at  a  distance  of  ^ref  from  the 

explosion  should  be  jointly  taken  as  (for  Y^  »  1  KT) 

1/3 

R  F  ■  2  miles 
ref 

* 

Apre^  *  34  mhar 


t 

re 


f 


-  0.33X2 


1/3 


sec 
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If  w*  wish  to  have  a  linear  source  model  such  as  is  implicit  in  the 


existing  theories,  we  would  have  to  interpret  ^Pref 


to  be  twice  that 


overpressure  which  would  be  observed  if  the  ground  were  not  present.  This 

means  we  would  want  Ap  ,  -  17  mbar.  If,  in  addition,  we  were  to  assume 

rcf  1/3 

that  spherical  spreading  holds  reasonably  well  between  1  and  2  miles,  then 

one  would  take  APref  "  17*2^3  when  one  backs  ip  to  R _ t  -  1  mile. 


ref 


All  this  suggests  that  the  program's  output  should  be  more  satisfactory 
for  near  surface  bursts  were  one  to  .jointly  set 


R  t  ■  1  mile 
ref 

Ap^j.  -  17*23^3  ■  21.4  mbar 

t  ,  -  0.32  *21/3  -  0.42  sec 
ref 


in  the  existing  computer  program.  The  method  in  which  this  is  done  is 
as  described  in  Sec.  3.1,2. 


3.3  VARIATION  OF  WAVEFORMS  WITH  ENERGY  YIELD  AND  WITH  BURST  HEIGHT 


One  of  the  principal  results  obtained  with  the  use  of  the  modal 
synthesis  program  during  the  past  reporting  period  was  a  clarification 
of  the  effects  of  burst  height  end  energy  yield  on  the  waveforms  generated 
by  nuclear  explosions  (Pierce,  Posey,  and  Iliff,  1971].  Literature 
published  previous  to  this  had  certain  discrepancies  which  were  apparently 
not  well  understood.  Scorer's  and  Hunt,  Palmer,  and  Penney's  results  would 
indicate  that  the  earliest  portion  of  the  pulse  received  at  large  distances 
(2000  km  or  Rreater)  should  have  amplitudes  directly  proportional  to 
energy  yield.  On  the  other  hand,  a  plot  (  which  we  omit  here  for  brevity) 
of  the  amplitudes  computed  numerically  by  Harkrider  ^Figure  15  in  his  paper) 


would  suggest  that  these  should  vary  with  yield 
above  5  MT. 


as 


,0.6 


for  yields 


The  two  relevant  sets  of  calculations  on  burnt  height  effects  exhibited 
by  Weston  [1962]  indicate  that  the  amplitudes  decrease  drastically  as  burpt 
height  is  increased  (although  this  fact  is  not  noted  explicit!)  in  Weston's 
text).  His  figures  indicate  the  amplitudes  of  the  earlier  portion  of  the 
pulse  to  be  down  by  factors  of  about  1/60  and  1/300  for  bursts  at  39  km  and 
76  km,  respectively,  from  the  amplitudes  for  a  burst  at  ground  level.  In 
contrast,  Harkrider's  altitude  effect  prediction  as  exhibited  in  Figure  17 
of  his  paper  (with  the  correction  suggested  by  Pierce  [1965  ])  indicates  that. 
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for  altitudes  up  to  17  km,  the  decrease  in  amplitude  with  burst  height 
is  much  smaller  than  one  would  expect  from  an  extrapolation  of  Weston's 
results.  Hnrkrlder's  results  would  suggest  the  net  decrease  to  only 
be  of  the  order  of  30%  when  the  explosion  is  raised  trom  2  to  17  km. 

Our  findings  by  use  of  INFRAS0N1C  WAVEFORMS  tend  to  more  nearly 
support  the  predictions  of  Scorer  and  of  Hunt,  Palmer,  and  Penney  as 
regards  yield-amplitude  proportionality  than  the  predictions  of  Harkrider's 

dependence.  As  regards  the  height  of  burst  dependence,  we  agree 
neither  with  Harkrider  nor  with  Weston.  We  find  a  small  increase  of 
amplitude  with  increasing  height  up  to  an  altitude  of  the  order  of  AO  km 
followed  by  a  decrease  with  increasing  altitude. 

In  Fig.  3,  a  number  of  waveforms  computed  for  fixed  height  of  burst 
apd  fixed  observer  location,  hut  for  variable  energy  yield,  are 
simultaneously  plotted.  In  each  case,  the  overpressure  divided  by  the 
yield  (in  megatons)  is  plotted  rather  than  the  overpressure.  Thus,  the 
near  coincidence  of  the  curves  during  the  early  portion  of  the  wave  form 
indicates  that  this  portion  of  the  wave  form  is  nearly  directly  proportional 
to  energy  yield. 

That  tills  proportionality  should  hold  in  the  range  of  lower  yields 
is  evident  from  the  general  theoretical  formulation.  If,  during  the  early 
portion  of  the  wave  form,  the  bulk  of  the  contribution  ho  the  dominant 
modes  comes  primarily  from  lower  frequencies  where  mT  <<•  l,  then  the 

3‘ 

sum  of  the  modal  waveforms  will  appear  proportional  to  T  (and  hence  to  E) 

at  tl.  ise  timi-s. 

An  estimate  of  the  range  of  yields  over  which  this  proportionality 


should  hold  may  be  obtained  if  one  takes 


in  the  condition 


to  be  the  frequency  corresponding  to  the  period  between  the  first  two  prominent 
peaks  in  the  wave  form.  Since  this  period  is  typically  about  A  minutes 
(see,  for  example,  the  data  exhibited  by  Donn  and  Shaw  [1967))  one  requires 
Ty  <<  AO  sec,  or,  for  sea-level  bursts,  F  <<  500  MT.  The  condition  becomes 

much  more  stringent  when  the  heig  of  hurst  is  raised. 

It  is  suggested  that  some  slight  improvement  over  the  yield-amplitude 
proportionality  law  might  he  obtained  if  one  took  the  amplitude  to  vary 
with  yield  Y  as 


V"  +  ‘W2'3' 


(3.13) 


where 


is  independent  of  yield  and  where 


is  given  by  (units  in 
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520  530  540  550  560 

TIME  AFTER  BLAST  (min) 


Figure  3.  Acoustic  pressure  per  megaton  of  the  source  yield  for  ground 
bursts  of  1,  10,  and  100  megatons  as  a  function  of  time.  The  observer 
is  on  the  ground  10,000  km  .north  of  the  source  in  a  subtropical  summer 
atmosphere. 
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kilotons  of  TNT) 


Y 

o 


(r  /p  f)(c/c  f)3/(«o  t  ,  ) 
1 o  'ref  ref  char  ref 


(3.14) 


where  p  and  c  are  evaluated  at  r.  and  where  u  ,  is  a  characteristic 
*o  o  char 

frequency  designating  the  early  portion  of  the  wave  form.  Typically,  Y 

should  he  of  the  order  of  5  X  10"  KT  for  sea-level  explosions.  For  u 

the  amplitudes  (first  principal  peak  to  trough)  of  the  computed  wave  forms 

shown  in  Fig.  3,  we  find  that  the  host  fit  to  Eq .  (3.13)  Rives  I)  »  0.0075 

dynes/cm^  KT  and  Y  -  4.4  K  lO^1  KT.  The  dogre-e  of  fit  is  Indicated  in 

®  Of' 

Ftp.  4.  For  contrast,  the  line  p  “  P*  Y  *  (derived  from  Harkridrr's 

KT 

computations)  is  also  plotted,  where  O'  is  selected  such  that  the  curve 
agrees  with  (3.13)  at  a  yield  of  5  MT.  The  discrepancy  shows  that  this 
model  more  nearly  exhibits  yield-amplitude  proportionality  than  does 
Harkrider's.  For  reasons  discussed  during  the  derivation  of(kl2)  we  are 
not  sure  that  our  source  model  is  adequate  in  the  quantitative  description 
of  the  deviation  from  yield-amplitude  proportionality.  However, the  fact 
that  the  computed  wave  forms  conform  to  this  proportionality  over  such  a 
wide  range  of  yields  tends  to  support  our  thesis  that  the  model  is  insensi¬ 
tive  to  the  choice  of  Rrejr  provided  that  the  constant  K  in  Eq.  (3.10)  is 
correct. 

In  lig.  5,  three  computed  waveforms  ate  shown  for  a  given  yield,  P 
given  observer  location,  and  a  given  model  atmosphere,  but  for  a  variable 
height  of  burst.  The  amplitude  of  the  entire  wave  train  is  seen  to  increase 
with  height  of  burst  below  some  altitude  (seen  in  Fig,  6  to  he  a  fuvict ion 
of  yield),  above  which  the  amplitude  decreases.  The  deci  in  amplitude 

is  greater  for  the  later  arriving  high-f requcncy  'acoustic'  portion  of 
the  wave  form. 

in  Fig.  6,  the  normalized  amplitude  p  /Y  (see  ig.  5)  is  plotted 
versus  height  of  burst  for  a  number  of  yie?cf§  Y.  A1  no ugh  the  details 
vary,  the  general  shape  of  all  such  curves  is  similar.  At  lower  altitudcj, 
the  amplitude  increases  slowly  with  height  of  hurst,  tenches  a  maximum  at 
a  height  somewhere  between  25  and  50  km,  then  decreases  at  a  i datively 
rapid  rate  with  subsequent  increasing  height, 

The^neral  shape  of  the  amplitude-burst  height  curves  may  he  interpreted 
with  the  aid  of  the  general  formulation  and  a  few  additional  considerations. 

For  the  most  part,  we  expect  that  the  earlier  low-frequency  portion  may  be 
considered  as  carried  by  a  single  composite  mode  [Bretherton,  1969;  Garrett, 
1969a,  b]  which  is  the  real  atmosphere's  counterpart  ol  the  edge  mode  discovered 
by  Lamb  [1910]  for  the  isothermal  atmosphere.  That  the  numerical  computations 
of  dispersion  curves  suggest  the  existence  of  such  a  mode  has  been  pointed 
out  by  Tress  and  Harkrlder  [1962],  who  use  the  descriptive  phrase  'pseudo 
mode'  to  denote  the  fact  that  its  dispersion  curve  coincides  with  those 
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NORMALIZED  PRESSURE  AMPLITUDE,  p  /D  <MT) 
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Figure  4.  Plot  of  the  antplitude  dependence  on  yield  found  in  the  present 

JomnuJ^r8  3  ^irVe  Sh°V'lng  thS  Y°'6  dePe^ce  derived  from  Harkrider 
computations  and  made  to  agree  with  the  first  curve  at  5  MT. 
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Figure  5.  Time  variation  of  ground-level  acoustic  pressure  10,000  km 
north  of  1  MT  explosions  at  altitudes  of  0,  50,  and  70  km  in  a  subtropical 
summer  atmosphere. 
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HEIGHT  OF  BURST,  (KM) 


of  mathematically  distinct  modes  over  different  frequency  ranges. 

In  the  paper  by  Pierce,  Posey,  and  Iliff  [1971]  it  is  argued  on  this 
basis  that  the  lower  frequency  portion  of  the  waveform  should  vary  with 
height  of  burst  zq  as 


amp 


(y-1)/y 


/U  +  (u 


char  Y 


T.,)2] 


(3.15) 


where  J  is  a  constant,  independent  of  hurst  height  and  of  yield.  Note 
that  does,  however,  depend  on  height. 

(y-1)/y 

Hie  factor  [1/p^]  in  (3.15)  increases  monotonically  with  height, 

although  at  a  relatively  slow  rate  compared  with  l/po  since  (y  -  l)/y  Is 

only  2/7  .  Thus,  were  it  not  for  the  denominator  [1  +  (ooTy)2],  ehe  quantity 

(Eq.  3.15)  would  increase  indefinitely  with  height.  As  it  is,  the  expression 
reaches  a  maximum  at  some  finite  height  ,  then  decreases  due  to  the  increase 
of  Ty  with  height  of  burst  as  1/p  .  The  altit-ude  at  which  this  maximum 

is  obtained  depends  on  yield  and  clearly  decreases  with  increasing  yield.  This 
is  in  accord  t/ith  the  results  shown  in  Fig.  6. 

To  check  to  what  extent  Eq .  (3.15)quantitatively  accounts  for  the 

variation  exhibited  in  Fig.  6.  we  plot  a  selected  portion  of  these  curves 

in  Fig.  7  along  with  corresponding  plots  of  Eq. (3,15)  The  values  w  ,  =  0.018 

—1  3  char 

sec  and  J  ~  3^0  ubars/MT  were  determined  taking  »  4.4XlOJ  MT 

and  D  -  7.5  ubars/MT  for  a  ground-level  burst.  Considering  the  limitations 
of  the  approximation,  we  consider  the  agreement  to  be  relatively  good. 

In  effect,  a  good  prediction  of  amplitude-height  dependence  was  made  using 
only  a  knowledge  of  amplitude-yield  dependence  for  a  ground-level  burst. 

It  should  be  pointed  out  that  the  validity  of  the  source  model  is  open 

to  question  (more  so  than  for  low-yield  low-altitude  hursts)  when  u>  ,  T„ 

char  1 

becomes  comparable  to  1.  Nevertheless,  it  is  gratifying  that  the  model 
does  predict  the  amplitude  to  eventually  decrease  when  the  burst  height  is 
raised  beyond  a  certain  point.  Any  other  behavior  would  certainly  be 
unacceptable  on  physical  grounds. 

We  are  unable  to  explain  the  discrepancy  between  our  derived  height  of 
burst  effect  and  those  presented  in  the  figures  in  Weston's  [1962]  paper. 

We  have  not  found  any  errors  In  Weston's  formulation  which  might  account 
for  this.  However,  In  view  of  the  fact  that  we  can  quantitatively  account 
for  discrepancies  between  our  results  and  those  obtained  independently  by 
Harkrider  [1964],  we  are  of  the  opinion  that  there  was  a  coding  error  in 
the  computation  of  Weston  b  wave  forms. 
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Figure  7.  Comparison  of  amplitude  dependence  upon  height  of  burst  as 
predicted  by  the  Lamb  edge  mode  hypothesis  and  ty  syntheses  based  on  the 


summation  of  fully-ducted  normal  modes. 
0.0185  sec-^ 
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Chapter  IV 


AN  APPROXIMATE  THEORY  BASED  ON 
LAMB'S  ATMOSPHERIC  EDGE  MODE 


4.1  INTRODUCTION 

The  Lamb  edge  mode  theory  has  previously  been  described  in  some 
detail  in  the  papers  by  Pierce  and  Posey  [1971],  Posey  and  Pierce  [1971 j, 
Posey  and  Pierce  [1972],  and  in  the  doctoral  thesis  by  Posey  [1971],  The 
present  chapter  summarizes  some  of  the  salient  aspects  of  this  theory 
and  its  implications. 

The  basic  idea  of  the  theory,  in  somewhat  simplified  terms,  is  that 
the  wave  may  be  regarded  as  traveling  along  the  ground  rather  than  through 
the  atmosphere.  The  situation  is  analogous  to  water  waves,  which  one 
regards  as  waves  on  the  surface  of  water  even  though  there  is  in  principle 
some  water  movement  at  arbitrary  depths  below  the  surface.  That  such  a 
description  is  possible  for  pressure  pulses  in  the  atmosphere  was  demon¬ 
strated  for  the  special  case  of  a  constant  temperature  atmosphere  hv  Lamh 
in  1910. 

However,  the  counterpart  of  this  type  of  motion  for  realistic 
(nonisothermal)  atmospheres  was  not  clearly  recognizable  in  theories  of 
wave  propagation  until  Bretherton  (1969)  and  Garrett  (1969a, h)  derived 
formulas  governing  the  propagation  of  constant  frequency  (sinusoidal)  waves 
of  this  type.  Of  key  importance  in  this  respect  is  the  wave  dispersion, 
whereby  waves  of  lower  frequency  travel  faster  than  waves  of  higher 
frequency.  The  Lamb  type  of  wave  (called  Lamb's  mode)  does  not  exhibit 
any  dispersion  for  the  isothermal  atmosphere  and  this  would  clearly  he 
inadequate  to  explain  the  data  since  the  shape  of  the  waveform  changes 
radically  with  increasing  distance  from  the  explosion.  Bretherton  and 
Garrett's  contribution  was  to  show  how  one  might  derive  formulas  govern¬ 
ing  this  dispersion  in  a  fairly  simple  manner, 

Clven  the  basic  idea  chat  the  wave  is  propagating  as  a  surface  wave 
in  the  real  atmosphere's  counterpart  of  Lamb's  guided  mode,  one  can 
consider  the  wave  arriving  at  a  given  point  to  have  travelled  along  a 
curved  path  (horizontal  ray  path)  which  connects  the  explosion  and  the 
recording  site.  [See  Fig.  8.]  The  determination  of  these  paths,  may  he 
made  using  the  same  principles  which  explain  why  light  rays  are  bent  when 
they  pass  from  air  to  water.  The  speed  of  waves  in  traveling  along  these 
paths  is  an  average  of  the  sound  speed  and  wind  velocity  profiles  above 
the  path.  The  distortion  of  the  waveform  as  it  moves  along  the  path  is 
governed  by  an  equation  similar  to  that  used  in  the  study  of  sonic  boom 
shock  waves  in  the  atmosphere  and  to  that  derived  bv  Korteweg  and  de  Vries 
in  1895  to  explain  why  water  waves  do  not  necessarily  break.  As  the  wave 
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Figure  8.  Model  of  edge  mode  propagation  from  an  explosive  source  in  an 
almost-stratified  atmosphere  above  a  spherical  earth. 


moves  from  point  to  point  along  the  path,  its  amplitude  continuously 
adjusts  itself  such  that  wave  action  (a  quantity  similar  to  energy) 
is  conserved.  The  pertinent  formulas  representing  just  how  this  guided 
wave  is  excited  by  the  detonation  of  a  nuclear  explosion  are  derived 
by  matching  the  guided  wave  onto  a  solution  valid  at  close-in  distances. 

An  important  result  of  the  theory  is  that  the  shape  of  the  waveform 
recorded  at  any  given  point  (providing  the  clock  recording  the  wave  is 
speeded  up  or  slowed  down  just  the  right  amount)  depends  on  just  one 
number,  the  ratio  of  a  characteristic  time  for  the  explosion  itself  to 
a  quantity  (also  having  the  units  of  time)  which  is  an  average  of  the 
properties  of  the  atmosphere  over  the  range  of  profiles  above  the  path. 

The  amplitude  of  the  waveform  depends  on  a  number  of  nuantities,  the  most 
important  probahly  being  the  spread  of  the  ray  paths.  If  a  number  of 
paths  should  be  caused  to  focus  (just  as  light  is  focused  by  a  magnifying 
glass)  larger  amplitudes  may  be  expected.  This  appears  to  be  the  explantion 
of  the  data  reported  by  Wexler  and  Hass  [1962] . 

In  regards  to  the  relation  of  the  Lamb  mode  theory  to  the  multimode 
theory  on  which  INFRASONLC  WAVEFORMS  is  based,  a  principal  virtue  of  the 
Lamb  mode  theory  is  its  simplicity.  One  may  note  that  the  development  of 
the  theory  of  long-range  atmospheric  pressure  pulse  propagation  since  the 
early  pioneering  work  of  Rayleigh  [1890]  ,  I.amh  [1910]  ,  Taylor  [  1  929]  , 

Solberg  [19116],  Pekeris  [1939,1948]  and  Scorer  [1  >0)  has  largely  been 
towards  progressively  more  complicated,  albeit  realistic,  models  [e.g. 

Yamamoto  1957;  Hunt,  Palmea  aid  Penney  1960;  Oazaryan  1961;  V.  H.  Weston  1961, 
1962  ;  Pfefferand  Xavichnv  1962;  Press  and  Harkrider  1962;  Weston  and  Van  Hulsteyn 
1962;  Clarke  1963;  Harkrider  1964:  Van  Hulstevn  1965;  Pierce  1965  ,  1967  ; 
MacKinnon  1967,  1968  ;  Balachandranand  Donn-1968;  Nelson,  Singhaus  and 
McNiel  1968;  Bal - -handran  1968,  1970;  Pierce,  Posey  and  Iliff  1971.]  Notable 
exceptions  to  t.  a  general  trend  are  the  works  of  Rcw  [1967],  Hines  [1967] 
and  Tolstoyand  Herron  [1970],  although  these  are  concerned  with  waves  which 
are  propagating  primarily  in  the  ionosphere  rather  than  near  the  Earth's 
surface . 

Even  with  the  aid  of  a  large  digital  computer,  agreement  with  data 
of  multimodal  calculations  tends  to  he  spasmodic.  The  reasons  for  this 
are  not  entirely  clear;  It  may  be  due  to  incorrect  assumptions  concerning 
the  atmospheric  profiles  or  concerning  the  source  mechanisms,  to  the  neglect 
of  horizontal  variations  ■*«  the  atmosphere,  or  of  dissipative  processes  in 
the  upper  atmosr '  "o,  <  '  accumulat i\ e  far-field  non-linear  effects.  Since 

the  implicatix  .a  of  existing  multi-mode  theories  such  as  those  of  Harkrider 
[1964],  MacKinnon  [1967]  or  Pierce  and  Posey  [1970]  are  difficult  to  explore 
without  the  aid  of  a  large  digital  computer,  the  inverse  problem  of  determining 
imperfectly  known  aggregate  properties  of  the  source  or  atmosphere  from 
empirical  waveforms  appears  formidable.  Even  if  the.  computer-assisted 
solution  of  the  inverse  .jem  were  feasible,  one's  confidence  in  its 
predictions  should  be  at  best  marginal  In  view  of  the  inherent  meagre  insight 
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affot^  -d  by  a  theory  which  relies  on  a  lengthy  sequence  of  computations, 
which  tends  to  obscure  the  causal  relationships  between  input  and  output 
variables.  One  does  not  know,  in  addition,  lust  what  a  priori  realistic 
extensions  to  the  existing  simplified  model  are  necessary,  without  first 
making  such  extensions  and  then  exploring  their  consequences  by  numerical 
experiment.  A  possible  trend  could  be  towards  an  overly  detailed  model, 
rather  than  towards  the  achievement  of  the  simplest  possible  model  adequate 
to  explain  the  existing  data.  That  the  latter  is  a  desirable  goal  is  a 
premise  which  with  we  believe  most  investigators  will  agree. 

One  theme  which  has  pervaded  much  of  the  discussions  of  the  empirical 
data  [see,  for  example,  Donuand  Ewing  1962]  is  that  the  early  portion  of 
the  waveform  may  he  considered  as  being  that  of  a  single  composite  mode. 

Press  and  Harkridcr  [1962]  noted  that  the  phase  velocity  versus  period  curves 
for  a  realistic  atmospheric  profile  tend  to  he  arranged  such  that  a  single 
curvt  of  nearly  constant  phase  velocity  extending  over  a  wide  range  of 
periods  may  be  constructed  by  the  connection  of  the  horizontal  segments 
of  a  number  of  individual  modal  dispersion  curves.  This  composite  mode 
is  evidently  the  real  satnoephere ' s  counterpart  of  the  edge  mode  (with 
constant  phase  velocity  equal  to  the  sound  speed)  predicted  by  Lamb  [1910] 
for  the  isothermal  atmosphere  model.  The  fact  that  Scorer  [1950]  found 
only  one  mode  for  a  profile  with  the  atmosphere  above  the  tropopause 
idealized  by  an  Isothermal  halfspace  suggests  that  the  other  portions  of 
the  Press  and  Harkrider  dispersion  curves  correspond  to  waves  guided  at 
relatively  high  altitudes.  This  has  been  supported  by  various  numerical 
experiments  [see,  for  example,  Harkrider  aid  Wei Is  1968]  as  well  as  by 
examination  of  the  height  profiles  associated  with  different  lrequencies 
and  different  modes  [Pfefferand  Zarichny  1963]. 

As  mentioned  above,  the  development  of  an  analytical  theory  for  the 
I ropagation  of  I.amb's  edge  mode  in  realistic  atmospheres  was  initiated  only 
relatively  recently  by  Bretherton  [1969]  and  by  Garrett  [1969a,  b] . 

Bretherton  showed  how  the  mode's  phase  velocity  may  be  found  to  first  order 
in  deviations  of  the  sound  speed  profile  from  that  of  an  isothermal  atmosphere 
and  to  first  order  in  the  wind  velocity,  while  Garrett  extended  the  method 
to  higher  orders  and  thus  succeeded  in  finding  the  dispersion  of  the  mode. 

The  latter  paper  by  Garrett  [1969b]  gives  a  detailed  discussion  of  the 
applications  of  the  theory  based  on  the  edge  mode  to  the  analysis  of  empirical 
waveforms,  as  well  as  of  the  limitations  of  the  theory.  In  spite  of  the 
relative  simplicity  of  the  Bretherton-Oarrett  model,  it  appears  to  be 
adequate  for  the  explanation  of  the  gross  propagation  speeds  and  dispersive 
characteristics  of  empirical  waveforms,  although  it  does  not  (at  least  to 
the  extent  which  the  perturbation  series  has  yet  been  carried)  explain 
more  subtle  features  such  as  the  inverse  dispersion  suggested  by  the  data 
analysis  of  Domand  Ewing  [1962].  , 

The  work  by  the  present  authors  in  regards  to  the  edge  mode  theory  has 
concentrated  on  the  development  of  a  relatively  simple  model  for  theoretical 
synthesis  of  transient  waveforms  recorded  at  large  distances  from  nuclear 
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explosions.  The  model  developed  would  appear  to  give  a  relatively  simple 
explanation  of  (i)  the  magnitude  of  the  amplitudes  of  the  first  few  peaks, 

(li)  the  periods  between  the  first  few  peaks,  (ill)  the  fact  that  records 
almost  always  begin  with  a  positive  pressure  rise,  and  (iv)  the  anomalies 
in  the  spatial  distribution  of  recorded  peak  overpressures  that  have  been 
reported  by  Wexlarand  Hass  [1962].  (Item  (iii)  above  is  unfortunately  not 
widely  appreciated,  due  to  the  fact  that  noise  often  obscures  the  early 
part  of  the  record  and  that  some  of  the  published  records  have  inadvertently 
been  published  with  the  positive  pressure  axis  extending  downwards.  Professor 
Honn  has  informed  us  that  this  was  the  case  in  the  paper  by  Donn  and  Shaw  [1967]. 
The  first  three  items  are  also  explained  by  the  more  elaborate  multimode 
theories  [see,  for  example,  Harkrider  1964]  and  we  accordingly  demonstrate 
that  the  edge  mode  model  gives  results  similar  to  the  multimode  calculations 
of  Pierce,  Posey,  and  Illff  [1971]  for  the  first  few  cycles.  The  advantage  of  the 
edge  mode  model  is  its  simplicity,  which  allows  an  exploration  ot  effects 
whose  consideration  may  have  been  prohibitively  difficult  within  the  context 
of  existing  multimode  theories.  Such  effects  include  accumulative  far-field 
non-linear  effects,  the  effects  of  terrain  variations  and  the  effects  of 
horizontal  atmospheric  variations. 

4..?  DISPERSION  OF  LAMB'S  MODE 

In  our  development  of  a  theory  for  the  propagation  of  Lamb's  edge  mode, 

we  considered  a  sequence  of  idealized  models.  The  mode's  dispersion  relation 

was  taken  for  a  temperature  and  wind-stratified  atmosphere  above  a  rigid  flat 

ground  at  z  "  z  ,  with  only  slight  and  unimportant  modifications  from  that 
fi 

previously  derived  by  Garrett  [1969a,  bj.  The  resulting  expression  relating 
angular  frequency  w  and  horizontal  wave  number  k  as  given  by  Pierce  and 
Posey  [1971]  is  to  first  order  in  e  given  by 


<o  ■  (Cj  +  vLk>k  +  0(r  ), 


(4.1) 


Here  c  is  any  convenient  parameter  character izing  the  order  of  magnitude 
of  (c(z)  -  Cj]/cj  where  c(r)  is  the  height  dependent  sound  speed  and 

(L  for  Lamb)  is  some  constant  value  representative  of  the  lower  atmosphere. 
One  may  choose  c.  and  also  a  mean  wind  velocity  in  such  a  manner  that 

L  l-t 

they  are  given  by  the  equations 


where,  for  any  f(z). 


1/Cj2  =  ^1  /c7};  Vj  -^v^, 

f  *  jfidz 


(4.2) 


(4.3) 
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Wltu  ft  ?/ 

fc  “  p  Y(z)/p  (z)  (4.4) 

representing  the  height  profile  of  what,  in  the  isothermal  case,  may  he 
construed  as  the  'netic  energy  density  for  l,amh*8  edge  wave.  Here  p  (z) 
and  p0(z)  represent  the  ambient  pressure  and  density  of  the  actual 
atmosphere.  The  integration  limits  should  he  and  some  jel^tively  large 
value  of  z  .  In  Fq .  (4.1),  v  ,  is  an  abbreviation  for  v  .0  where 

Im  li  K 

e^  is  the  unit  vector  pointing  in  the  direction  of  the  horizontal  wave 
*♦ 

number  k.  Thus,  to  first  order  in  c  ,  the  dispersion  relation  is  that 
appropriate  to  an  isothermal  atmosphere  (sound  speed  Cj )  with  constant 

horizontal  wind  velocity  v"j  . 

The  definitions  one  may  adopt  for  c^  and  v^  are  restricted  only 

to  first  order  in  r.  There  is  no  a  priori  reason  for  preferring  nnv  one 
of  alternate  definitions  which  differ  in  second  order.  The  definitions  we 
took  were  chosen  such  that  the  resulting  second  order  terms  (which  actually 
govern  the  dispersion)  arc  of  relatively  simple  form. 

To  second  order  in  e  ,  the  dispersion  relation  obtained  by  Garrett 
[1969a]  may  be  rewritten  (albeit  with  some  lengthy  manipulation)  in  the 
fora 


where 


with 


0) 

'  k('L  +  v..k  *  "kk>-,‘<k2  -  kBi>kk- 

(4.5) 

akk 

-  IV(2cl)]^[(v(z)  -  vL)  .  ek]2) 

(4.6a) 

1  5  V,  3,2\ 

hkk 

-  I  CL  VC  *  17 

(4 . 6h) 

C(z) 

7 

•  £,  1  (z)  |<£(z)(cl  2  -c  2)dz 

(4.6c) 

P 


z 

D(z)  1(z)  |  £(z)  ($  -  "vj  )dz.  (4 . 6d) 

£ 


It  should  be  noted  that  the  latter  two  quantities  both  remain  bounded  at 
high  altitudes  by  virtue  of  the  definitions  in  F.qs.  (4.2-4).  The  quantity 
kBL  whlch  aPPears  in  F.q  .  (4.5)  is  <drl/cj  ,  where  w  -  (y  -  l)*'^g/c 


is 
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the  Brunt -Vaieala  frequency  for  an  isothermal  atmosphere  with  sound 

speed  Cj  .  The  form  of  the  second-order  terms  in  the  dispersion  relation 

would  tend  to  suggest  that  the  approximation  is  best  for  frequencies  of 

the  order  of  the  Brunt-Vaisal a  frequency,  rather  than  in  the  limit  of  low 

frequency.  Fortunately,  the  dominant  frequencies  in  pressure  waveforms 

at  large  distances  from  nuclear  explosions  are  of  the  same  order  of  magnitude. 

[A  detailed  derivation  of  the  above  equations  is  given  by  Posey  (1971).] 

One  may  note  that  Fq.  (A. 5)  may  also  be  written  in  the  form 


where 


w  “  CFk  '  hkkk 


CL  +  VLk  +  akk  +  hkkkBL 


(A.  7) 


(A. 8) 


In  order  to  facilitate  comparison  of  F,q .  (A.  7)  with  the 
determined  multi-mode  dispersion  curves,  it  may  be  expressed 
expressing  frequency  <o  as  a  function  of  phase  velocity  v  ** 


numerically 
ir>  a  form 
m/k. 


v_ f (v  - 


<v>/t 


.1/2 


(A.  9) 


Figure  9  is  a  representative  comparison  between  this  relation  and  the 
multi-mode  curves  for  the  same  atmosphere.  The  lower  225  km  of  the  model 
atmosphere  used,  as  listed  in  Table  1  ,  was  selected  to  represent  average 

conditions  between  Johnston  Island  and  Berkeley,  California  for  the  month 
of  October,  based  on  data  from  the  Hflnf|bnnk  of  Cenphyslra  .and.  Spaco- 
F.nvi ronments  [Valley,  1965]  and  from  the  1965  COSPAR  International  Reference 
Atmosphere  [COSFAR,  1965].  The  direction  of  propagation  is  along  the  great 
circle  path  from  Johnston  Island  to  Berkeley,  35°  north  of  east.  Since  the 
Lamb  mode  dispersion  is  only  weaklv  affected  by  conditions  above  115  km 
[Posey,  1971],  the  computation  of  Fq .  (A. 9)  uses  only  the  model  atmosphere 
below  this  altitude.  The  edge  mode  propagation  parameters  are  found  to 
be  -  0.318  km/ sec  and  h^  "  0.132  kmVsec.  Notice  that  the  agreement 


is  within  0.27,  in  the  flat  portions  of  GRC,  S-  and  S.  for  frequencies 
up  to  0.07  see"^  (periods  down  to  about  90  sec),  which,  as  discussed  by 
Posey  [1971]  ,are  relatively  Insensitive  to  the  atmosphere  above  110  km. 


Garrett  [1969b]  lias  discussed  the  effect  of 
boundary  impedance  at  a  given  altitude  z  -  H  . 

variation  in  the  phase  velocity  of  the  edge  mode 


imposing  an  arbitrary 
He  finds  that  the  relative 


is  cf  the  order 
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Figure  9.  Comparison  of  approximate  edge  mode  dispersion  curve  (dashed 
line)  with  guided  mode  dispersion  curves  obtained  by  numerical  solution 
of  the  acoustic-gravity  wave  eigenvalue  problem  (solid  lines).  The 
nomenclature  labelling  the  acoustic-gravity  wave  modes  is  that  of  Press  and 
Harkrider  (1962). 


(where  Hfi  is  a  scale  height)  for  plausible  variations  in  the  boundary 
impedance.  The  agreement  between  the  independently  derived  dispersion 
relations  showi'.  in  Fig.  9  seems  to  substantiate  this. 

A  relation  for  the  group  velocity,  v  »  dw/dk,  can  also  be  found 
from  Eq.  (4.7),  ^ 

Vg  “  CE  ~  3hkk(w/cE)2  (4’1Q) 


and  a  comparison  of  this  with  the  multi-mode  curves  for  the  same  model 
ae  described  above  is  shown  in  Fig.  10. 

A  FORTRAN  computer  program  which,  among  other  quantities,  evaluates 
the  various  Lamb  mode  propagation  constants  Cj  ,  v^  ,  hj,v  Cg,  etc.,  for 

an  arbitrary  multilayered  atmosphere  was  written  during  the  contract  and 
is  given  in  Posey's  doctoral  thesis  [  Posey,  1971  ], 


4.3  THEORETICAL  EXPRESSIONS  FOR  EDGE  MODE  WAVEFORMS 

The  general  results  of  the  theoretical  development  by  Pierce  and 
Posey  [1971]  of  the  excitation  and  propagation  of  Lamb's  atmospheric  edge 
mode  from  nuclear  explosions  Is  summarized  here.  The  general  expression 
foe  the  acoustic  pressure  signature  at  a  point  x  at  times  t  is  given  by 


p(x,t)  ~ 


C  DV(x,x0)gE  ? 

>  2 j2  2  *  Y/2  \  ♦(t*s,0), 

C<cL0V  cL0  Cdb/de>  J 


where  D  is  a  numerical  constant  given  by  (Kd0,50) 

D  -  K#<y  -  1)  (2  -  y)4irSC  0.013 
while  V(ic,x  q)  Is  a  dimensionless  quantity  given  by 


(4.11) 


(4.12) 


v(+  •*  ,  _  A(x)A(x0)cI0  1  (cL  +  vIk)0 _ 

0  ^(*o)cL1/2|CL?k  +^1/2(cL  f  VLk)1/2 


(4.13) 


with 


AC*)  -  v01,y<.t>ha(i/Uy>  - <1/2) 


(4.14) 
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Figure  10.  Comparison  of  multi-mode  and  edge  mode  group  velocity  curves. 
These  curves  are  appropriate  for  propagation  35°  north  of  east  in  the 
atmosphere  of  Table  I. 


Here  x  denotes  the  observation  location,  the  source  location, 
x^(or  xQg)  the  point  on  the  ground  directly  below  the  observer  (or 

source).  The  subscript  0  implies  the  corresponding  quantity  is  evaluated 
for  the  atmosphere  at  the  source  location.  The  quantity  g  in  F.q.  (4.11) 
is  the  acceleration  of  gravity,  E  is  the  total  energy  released  by  the 
explosion,  T^  is  the  source  characteristic  time  defined  by  Eq.  (3.3)  in 

the  present  report.  The  quantity  dt/dd  is  a  parameter  which  allows  for 
focusing  [Fig,  8]  or  defocusing  of  horizontal  ray  paths  caused  by 
departures  of  the  atmospheie  from  perfect  stratification  (i.e.  ,  horizontal 
variations)  and  formally  is  defined  as  the  ray  channel  width  per  unit 
angle  9  with  which  horizontal  ray  paths  exit  from  the  source  location. 

For  a  perfectly  stratified  atmosphere  but  with  the  earth  curvature  correction 
(bee  Sec.  3.1.1  of  the  present  report)  included,  it  is  given  by 

db/d0  "  sin(r/»-e)  (4.15) 

The  function  VKs.t,©)  which  appears  in  Eq.  (4.11)  is  given  by 
T 

«  - — 

w*- 

D 

where  Ai(x)  is  the  Airy  function,  defined  by 


(4.17) 


(4.18) 


gives  the  characteristic,  waveform  time  dependence  for  the  Lamb  edge  mode 
in  the  intermediate  field  (before  dispersion  has  had  an  appreciable 
accumulative  effect).  A  plot  of  this  function  is  given  in  Fig.  11  . 

In  Eq.  (4.11),  we  have  also  made  the  abbreviations 


Ai(x) 


-  f, 


cos(v  /3  +  xv)dv 


and  where 


•1, 


M(u)  -  ]  (1  -  £)exp(-  £)[y  -  £)"  1/2  d?  U(u) 
0 

ft 


\l\i  +  (1  -  2p)exp(-  y)  exp(y2)dy  AU(u) 


Ai  *  t  f  +  p  M(p)dp  (4.16) 
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(4.19a) 


8 


for  the  fwo  characteristic  functions  s  which  appear  in  the  integral. 

Both  of  these  quantities  have  the  uni  of  time  and  are  here  referred  to 

as  the  waveform  arrival  time  and  the  characteristic  dispersion  time* 
respectively.  The  parameter  s  is  related  to  distance  l  along  a 
horizontal  ray  path  (see  Fig.  8)  by  the  equation 


ds 

dt 


C  4-  v 

L  Lk 

rrr+in 

Vv  L 


C  4-  v 
I,  Lk 


+  2c  v 

L  Lk  L  J 


(4.20) 


For  typical  winds  in  the  atmosphere  it  would  appear  sufficient  to  take 
s  -  l  . 


In  the  often  realized  limit  of  propagation  distances  sufficiently 
large  that  /t^  <<  1  and  for  arrivals  sufficiently  early  that 

t  -  Ta  <<  t^/Ty,  Pierce  and  Posey  [l971J  argue  that  the  function  6 

may  be  suitably  approximated  by 


*  "  (2/m<TY/Tr)3/2PP([t  -  ^1/td)  (4.21) 

where  PP(x)  is  the  function 

Ot 

PP(x)  -  jo  Ai'(a2  -  x)da  .  (4.22) 

In  the  limit  represented  by  Eq.  (4.21),  the  acoustic  overpressure  p, 
as  given  by  Eq.  (4.11),  is  directly  proportional  to  the  energy  release  E  . 
This  is  clear  since  the  only  other  parameter  in  the  formulation  which  varies 
with  energy  yield  is  Ty  and  since  the  T^  '2  factor  in  Eq.  (4.21) 

cancels  the  corresponding  factor  in  the  denominator  in  Eq.  (4. 111.  Note 
that  this  is  in  accord  with  the  numerical  study  of  Pierce,  °osey,  and  II iff 
[l97lj  ,  the  results  of  which  are  summarized  in  Sec.  3.3  of  the  present 
report.  The  result  quoted  there  that,  for  smaller  yields  or  lower  heights 
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of  burst,  the  amplitude  of  ihe  first  fe<  peaks  varies  with  source  height 
(v_ I  }  ft 

Zq  as  1 1/p^ <Zq ) J  v  is  also  Indicated  by  Eq.  (4.11). 

The  value  0.50  of  the  dimensionless  constant  K  quoted  in  Eq.  (4.12) 

turns  out  to  be  applicable  to  both  above  ground  and  near  surface  explosions 

since  the  alternative  choices  of  Ap  ,  and  t  c  discussed  in  Secs.  3.2.1 

ret  ref 

and  3.2.2  are  such  that  APieftref^  is  t*ie  same  in  both  cases.  [See  Eq.(3  10).] 
However,  in  the  computation  of  as  given  by  Eq.  (3.3),  one  should  use 

t  ^  °  0.33  sec  and  cref  “  0.416  sec  for  above  ground  and  near  surface 
explosions, respectively. 

The  general  shape  of  the  Lamb  mode  waveform  is  fully  described  apart 
from  a  multiplicative  time  independent  quantity  by  the  function  ip  .  The 
numerical  plots  in  Figs.  12  and  13  should  give  a  general  indication  of  the 
predictions  of  Eq.  (4.16)  for  various  ty^td  •  The  1 curve *n  Fig-12  is  close 

to  the  limiting  case  represented  by  Eq.(4.21).  It  may  be  nooed  from  Fig.  1 3  that  the  theory 
predicts  the  first  peak  to  be  smaller  than  the  second  as  long  as  T  /r^  <  0.34. 

The  ratios  of  the  periods  to  in  the  first  few  cycles  would  appear  from 

Fig.  12  to  be  relatively  insensitive  to  the  value  of  TY^TD  *  Thus  the 

period,  first  peak  to  second  peak,  is  about  h.Gi^  while  the  period,  second 

peak  to  third  peak,  is  about  2.6td 

4.4  COMPARISON  WITH  MULTIMODAL  COMPUTATIONS 

A  computer  program  encompassing  the  edge  mode  theory  was  written  during 

the  reporting  period  by  Posey  fl971J .  The  program  is,  strictly  speaking, 

in  its  present  form  only  applicable  to  atmospheres «  '.thout  horizontal 

variations,  although  it  allows  the  user  the  option  of  specifying  three 

separate  model  atmospheres  for  the  computation.  One  atmosphere,  the  "source 

atmosphere"  could  be  utilized  in  computing  all  quantities  (0  subscript)  in 

the  coefficient  of  >p  in  Eq.  (4.11)  which  correspond  to  the  atmosphere  near 

the  source  and  in  the  computation  of  T^  .  Another  atmosphere  (the  "propagation 

atmosphere")  could  be  utilized  in  finding  those  quantities  appropriate  to  the 

atmosphere  near  and  above  the  observer  location.  Finally,  a  third  atmosphere 

could  be  utilized  in  the  computation  of  t  and  t  and  in  the  computation 

D  3 

of  i^(s,t,9).  The  quantity  db/d9  would  either  have  to  be  estimated  independently 
or  taken  to  be  given  by  Eq.  (4.15)  or  by  r  if  r  <<  . 

The  computer  program  mentioned  above  enabled  us  to  check,  for  the  case  of 
a  stratified  atmosphere,  the  extent  of  validity  of  the  single  mode  hypothesis 
(where  all  modes  other  than  the  Lamb  edge  mode  are  neglected)  by  comparison 
of  waveforms  calculated  according  to  the  Lamb  edge  mode  theory  with  those 
computed  by  the  multimode  theory  as  represented  by  the  computer  program 
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f/(ry/r0) 


Figure  12  Variation  of  the  waveform  factor  i  (t)  with  the  parameter 
Ty/ Tp  .  Here  ^  i  i  normalized  with  respect  to  Ty/ and  the  time 

is  plotted  relative  to  t  in  units  of  t_ 

a  D 
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Figure  13,  Predicted  ratios  a  3 /a ^  and  32^sl  of  yuccesslve  Peak  to 
trough  amplitudes  *3,  and  a3  in  far-field  edge  mode  waveforms 

produced  by  nuclear  explosions. 
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INFRASONIC  WAVEFORMS.  An  example  we  chose  for  comparison  is  that  where  the 
source  was  a  hypothetical  10MT  burst  at  3  km  above  Johnson  Island  and  the 
observer  was  at  ground  level  at  r  ■  5600  kra  at  Berkeley,  California. 

The  model  atmosphere  is  as  described  previously,  corresponding  to  the 
dispersion  curves  shown  in  Fig.  9.  The  numerical  calculations  give 
c^  •  0.312  km  s“  ,  v^  -  0.005  km  s“^  ,  the  source  characteristic  time 

T  as  12  s,  the  effective  low-frequency  sound  speed  c  as  0.318  km  s  ^  , 

Y  3  -1  e 

the  dispersion  constant  h^  as  0.13  km  s  ,  the  arrival  time  xa  as 

17600  a  (4.9  hr),  while  the  characteristic  dispersion  time  t  is 

60  s.  Thus  td 

In  Fig.  14,  the  first  15  min  of  (1)  an  experimental  waveform,  (2)  a 
waveform  computed  according  to  the  multimode  theory,  and  (3)  a  waveform 
computed  according  to  the  edge  mode  theory  are  simultaneously  plotted. 

The  experimental  waveform  is  that  recorded  at  Berkeley,  California  following 
the  megaton  class  Housatonic  shot  of  1962  October  30  and  which  is  given 
by  Donnaad  Stnw  11967]  in  an  inverted  form  (Professor  Donn  has  informed 
us  that  the  maximum  peak-to-t rough  amplitude  of  this  record  is  about 
350  wbar).  The  second  waveform,  based  on  the  multimode  theory, has  what 
we  regard,  al!  eit  subjectively,  as  a  relatively  good  agreement  in  waveform 
shape  and  period-time  variation  with  the  experimental  record.  The  comparison 
of  central  concern  here,  however,  is  that  between  the  second  and  third 
records.  One  should  note  that  over  the  first  three  cycles  the  agreement 
is  substantial.  The  zero  crossings  coincide  almost  exactly,  while  the  heights 
of  individual  peaks  or  depths  of  individual  troughs  agree  well  except  for 
a  discrepancy  of  30  per  cent  in  the  first  peak  and  first  trough.  The 
greater  disagreement  at  later  times  is  to  be  expected,  in  view  of  the  fact 
that  the  higher  frequencies  arriving  at  these  times  may  be  carried  in  other 
modes  corresponding  to  natural  channeling  within  an  altitude  region  of 
lower  effective  sound  speed. 

As  for  those  discrepancies  which  do  exist  at  earlier  times,  particularly 
the  'precursor'  dip  in  the  multimode  synthesis,  it  is  difficult  a  priori 
to  say  which  theory  gives  the  more  nearly  correct  prediction.  The  edge 
mode  theory  presented  here  is  sufficiently  simple  that  negligible  numerical 
err  s  arise  in  the  computation.  On  the  other  hand,  the  multimode  synthesis 
is  based  on  numerical  computations  of  Inverse  Fourier  transforms  (the 
integration  variable  being  frequency),  where  the  limits  of  integration  arc 
truncated.  The  truncations  Laving  the  most  significant  effect  are  those  at 
the  low  frequency  cut-offs  shown  in  Fig.  9  which  correspond  to  the  transition 
from  fully-ducted  modes  to  leaky  modes  [Harkrider  &  Wells  1968],  The  neglect 
of  these  leaky  modes  may  place  spurious  precursors  in  the  synthesized  waveforms 
with  periods  of  the  order  of  the  lower  cutoff  frequency,  although  hopefully 
this  neglect  will  have  negligible  effect  on  the  dominant  portion  of  the 
waveform  (the  high-frequency  part  of  the  apparent  precursor  is  due  to  the 
choice  of  upper  integration  limit).  In  contrast,  within  the  context  of  the 
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Time  after  blast  (min) 


Figure  14.  Observed  and  theoretical  pressure  waveforms  at  Berkeley, 
California,  following  the  Housatonlc  detonation  at  Johnson  Island  on 
October  30,  1962.  The  observed  waveform  Is  taken  from  Donn  and  Shaw  [1967] 
The  energy  yield  assumed  in  the  theoretical  computations  was  10  MT.  The 
maximum  peak  to  trough  amplitude  of  the  experimental  waveform  is  about 
350  ubar. 
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edge  mode  theory  presented  here,  the  edge  mode  is  formally  fully  ducting 
down  to  xero  frequency,  and  there  is  no  occasion  to  decide  whether  or 
not  one  should  Include  leaky  modes  in  his  computation  scheme. 


A. 5  IMPLICATIONS  OK  THE  EDGE  MODE  THEORY 

The  edge  mode  model,  being  a  simplified  approximate  model,  certainly 
cannot  give  a  detailed  explanation  of  the  existing  data,  and  in  many 
Instances  may  lead  to  predictions  which  disagree  significantly  with  the 
data.  However,  if  one  adopts  the  view  that  the  edge  mode  model  provides 
a  framework  for  discussion  of  the  data,  then  there  is  perhaps  some  merit 
in  listing  a  number  of  the  implications  of.  the  model,  with  various 
idealizations  as  to  the  atmospheric  structure.  If  the  predictions  do 
not  agree  with  the  data  in  some  respect,  then  we  may  seek  to  find  what 
modifications  or  extensions  are  necessary  to  the  theory,  rather  than  to 
abandon  the  edge  mode,  model  completely. 


4.5.1  Variation  of  waveform  period  with  distance 

The  theory,  in  the  normal  limit  when  Ty/t^cl,  predicts  that  the 

period,  first  peak  to  second  peak,  should  increase  with  distance  r  from 
the  source  as  fox  the  stratified  atmosphere  model.  This  follows  from 

Fig.  12  and  Eq.  (4.19b).  The  data  does  show  some  tiend  towards  longer 
periods  at  greater  distance^,  hut  the  spatial  variations  with  azimuthal 
angle  are  substantial.  Since  this  is  a  relatively  weak  dependence  on  r, 
it  is  difficult  to  check  quantitatively. 


4.5.2  Ratios  of  successive  periods 

The  waveform  factors  plotted  in  Fig.  12  suggest  that  the  ratio  of  the 
first  two  periods  (first  peak  to  second  peak,  second  peak  to  third  peak) 
should  be  4/2'.  6*1;  5.  On  the  whole,  the  agreement  with  data  of  this 
prediction  is  fair.  Periods  tend  to  decrease  with  increasing  time  within 
a  waveform,  and  in  those  instances  where,  the  peaks  may  be  clearly  Identified, 
the  ratio  varies  between  1.3  and  1*8,  For  example  the  New  Orlea  ,La., 

record  of  1961  October  30,  the  ratio  (depending  on  what  one  interprets  as 
the  first  and  third  peaks)  might  be  evaluated  as  1.5. 


4,5.3  Record  begins  with  pressure  rise 

That  the  theory  predicts  this  is  clear  from  the  curves  in  Fig.  12.  This, 
while  admittedly  a  qualitative  rather  than  quantitative  prediction,  is 
reasonably  well  supported  by  the  data.  For  example,  in  those  records 
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given  by  Donn  ct  al.  [1963]  where  the  positive  pressure  direction  is 
indicated,  and  where  one  can  distinguish  the  signal  onset  from  the  noise, 
the  first  positive  pressure  peak  precedes  the  first  trough. 


4.3.4  Ratio  of  successive  peaks 

If  T„/t^  is  small  (which,  as  may  be  noted  from  the  numerical  example 

in  the  preceding  section,  is  the  normal  case)  then  the  first  peak  should  be 
lower  than  the  second  peak  and  the  ratio,  first  trough  to  second  peak 
divided  by  first  peak  to  first  trough,  should  be  about  1.2.  The  agreement 
of  this  prediction  witli  the  data  is  negligible.  On  many  records  the  second 
peak  is  actually  -.mailer,  and  the  ratios  of  the  two  seem  random.  The 
disagreement  may  be  due  to  a  variety  of  causes:  instrument  response, 
interference  from  other  modes,  higher-order  terms  in  the  dispersion  relation, 
interference  from  ambient  noise,  etc. 


4.5.5  Variation  of  waveform  amplitude  with  distance 

In  the  case  when  there  is  negligible  focusing  of  horizontal  edge  mode 
rays,  the  theory,  Eqs.  (4.11)  and  (4.15)  ,  predicts  that  wave  amplitude  varies 
with  great  circle  distance  r  as  [sln(r/r^)]  while  individual  peaks  of 

the  waveform  factor  t// ,  for  the  strati^ed  atmosplicr^mode ,  decrease  with 

r  in  the  limit  of  small  7  /t  as  t  or  as  r"-*-'  *•  from  Eqs.  (4.19) 

’  u  u  -1/2 
and  (4.21).  Thus  amplitude  varies  in  the  idealized  case  as  r  (sin(r/re)J 

This  variation,  as  previously  discussed  by  Harkrider  f 1964} ,  agrees 
relatively  well  with  the  data  of  Wexler  and  H'-ss  [1962],  although  there  are 
substantial  variations  of  amplitude  with  azimuth  angle.  In  Fig.  15,  maximum 
amplitudes  of  synthesized  Lamb  mode  waveforms  [Posey,  1971]  are  compared 
with  similar  results  presented  by  Harkrider  [1964]  for  synthesized  multi- 
mode  waveforms  [multiplied  by  po(z^)/po(z^)  as  suggested  by  Pierce,  1965  ] 

and  with  the  empirical  data  presented  by  Wexler  and  Hass.  The  atmospheric 
model  and  the  source  parameters  used  in  the  Lamb  mode  syntheses  are  the 
same  as  used  by  Harkrider;  i.e.,  a  58  MT  explosion  at  an  altitude  of  3.66  km 
in  the  ARDC  standard  atmosphere,  as  shown  in  Harkrider's  Fig.  1,  with  no 
winds . 


4.5.6  Determination  of  Kxplosion  Energy  from  Pressure  Vaveforms 

One  of  the  more  intriguing  consequences  of  the  Lamb  edge  mode  theory  is 
that  it  affords  a  simple  method  which  is  relatively  Insensitive  to 
atmospheric  structure  for  estimating  an  explosion's  energy  yield  from 
experimental  pressure  records.  As  noted  in  Sec.  4.3,  the  shape  of  the 
waveform  factor  \ p  is  determined  solely  by  the  value  of  T^/x^  ’ 


1/2 
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O  Average  Observed  (Waxier  and  llabs) 
A  Lamb  Mode  Synthesis 
□  Multi-Mode  Synthesis  (llarkrider) 
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Figure  15.  Comparison  of  observed  and  theoretical  variation  of  maximum 
amplitude  with  distance  from  the  source.  The  results  of  a  series  of 
Lamb  mode  syntheses  are  compared  with  observations  reported  by  Wexler 
and  Hass  [1962]  and  with  the  results  of  waveform  syntheses  performed 
by  llarkrider  [1964]  using  a  multi-mode  formulation.  (Harkrider's  re¬ 
ported  amplitudes  are  here  corrected  by  the  factor  p^(zo)/pc(r  )  sug¬ 
gested  by  Pierce  [1965],  which  in  this  case  is  about  0.76.) 
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Moreover,  the  ratios  of  the  periods  to  in  the  first  few  cycles  would 

appear  from  Fig.  12  to  be  relatively  insensitive  to  this  ratio.  One  finds, 
in  particular,  the  period  Tp  ^  between  first  and  second  peaks  is  given 

by  4.0  .  Thus  one  can  replace  the  largely  aiiriori  unknown  quantity 

in  the  theoretical  expressions  by  _/4  .  The  plots  in  Fig.  12  also  show 

that  the  quantity  ,  representing  Che  first  peak  to  first  trough  amplitude 
of  if,  is  a  slowly  varying  function  of  T^/t^  '  A  relatively  good  approximation 

(within  10%)  for  most  cases  of  interest  is  to  set 

.y<vy3/2  *  °- 9  (4.23) 

With  the  above  approximations  in  mind, one  may  note  that  if  Eq.  (4.11)  is 
written 

p  -  Gil  (4,24) 


where  C  is  an  abbreviation  for  the  quantity  in  braces  in  Eq.  (4.11),  then 
the  absolute  pressure  variation  ppp^  from  the  first  peak  to  the  first: 

trough  should  be  given  by 

PFPT  “  G  aX  (4.25) 

or  ,  from  Eq.  (4,23)  and  the  approximation  j/4 »  as 

pppT  -  (G)(0.9)(4TY/T12)3/2  (4.26) 


When  the  expression  for  G  is  inserted  and  one  solves  for  E,  one 
finds 


E  - 


c  7/2(db/dG)1/2 

LO _ _ _ 

7.2  D  V  g 


(4.27) 


where  the  various  quantities  are  as  defined  in  Sec.  4.3. 

This  can  be  further  simplified  if  a  number  of  plausible  assumptions 
are  made:  (1)  the  atmosphere  is  assumed  to  be  horizontally  stratified; 

(2)  the  ground  level  pressure  at  source  and  observer  is  taken  to  be  one 
atmosphere;  (3)  the  observer  is  considered  to  be  on  the  ground;  and  (4)  wind 
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velocities  are  considered  small  compared  to  the  sound  speed.  In  this 
event,  one  finds  the  encigy  yield  to  be  given  by  (in  equivalent  kilotons 
of  TNT) 

Y  -  B  sin1/2(r/r  )  [p  (x  )  ]1_aA)  )  T.  JV2p_T  (4.28) 

e  O  O  O  1 ,  l  i'rl 


where  b  *  0.199  if  p  is  measured  in  atmospheres,  p_,D,-  is  measured  in 

O  rrl 

pbar,  c  in  km/sec,  T.  Q  in  seconds,  and  Y  in  KT. 


Equation  (4.28)  would  enable  one  to  estimate  tbe  yield  Y  if  the 
height  of  burst  were  known.  Usually,  this  is  not  the  case.  However,  the 
dependence  ot;  ambient  pressure  or  burst  height  is  relatively  weak.  The 
variation  only  amounts  to  30%  even  for  heights  as  great  as  10  km.  The  speed 
of  sound  at  the  source  con  also  vary,  but  generally  stays  within  15a  of 
310  m/'aec  for  the  first  60  km.  Thus,  would  appear  expedient  and  not 
too  gross  an  approximation  to  set  p  (x  )  equal  to  1  atm  and  to  set 

~~)r 

c(xf)  equal  to  310  m/sec.  When  this  is  dona  on~  finds 

Y  -  B‘  ain1/2(r/re)  ^lf212  PppT  (4.29) 

where  8'  has  c  value  of 

-1  -1/2 

B'  -  0.034  pbar  sic  KT  (4.30) 


Because  of  the  nature  of  the  approximation  outlined  above,  this  prediction 
of  energy  yield  may  be  too  large  for  high  altitude  explosions  or  for  records 
at  long  ranges  from  small  yield  explosions  (i.e.,  Ty/t^  <  0.2)  and  too 

small  for  records  at  short  ranges  from  large  yield  sources  (i.e., 

"VV  0.2). 

The  waveform  parameters  T.  0  and  p  read  from  an  observed 

1  9 &  rrl 

microbarogram  should  properly  be  used  in  the  above  yield  estimation  formula 
only  when  the  record  prior  to  the  second  peak  represents  almost  solely  the 
Lamb  mode  pressure  perturbation.  In  principle,  this  can  only  be  the  case 
when  there  is  negligible  background  noise,  when  there  is  not  yet  appreciable 
contribution  to  the  record  from  other  modes,  and  when  there  is  negligible 
phase  shift  due  to  instrument  response. 

In  order  to  check  the  agreement  between  the  approximate  Lamb  mode 
relation  and  the  empirical  waveforms,  the  quantities  P-p-pj  and  ^  ware 

measured  on  a  number  of  records  presented  by  Harkrider  [1964]  and  by  Donn 
and  Shaw  [1967].  The  explosion  sources  for  the  records  (generally  recorded 
at  more  than  one  geographic  location)  are  labeled  a  through  n  (see  Table  II). 
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event 

Bath’s  Yield 

date  Estimate  (MT)  R(km) 

Tj  2(sec) 

Pm(wbar) 

Source  of 
Microbarogram 

a 

9/10/61 

10 

3676 

193 

100 

DS 

6644 

280 

70 

DS 

b 

9/11/61 

9 

8000 

229 

125 

11 

c 

9/14/61 

7 

6569 

258 

61 

DS 

d 

10/4/61 

8 

8000 

220 

80 

H 

e 

10/6/61 

11 

6674 

300 

59 

DS 

8000 

390 

45 

H 

f 

10/20/61 

5 

8000 

310 

21 

11 

g 

10/23/61 

25 

6677 

376 

280 

DS 

h 

10/30/61 

58 

6688 

400 

500 

DS 

5631 

540 

133 

DS 

8000 

690 

140 

H 

i 

10/31/61 

e 

5717 

386 

25 

DS 

j 

5/4/62 

3 

5375 

200 

36 

DS 

k 

6/10/62 

9 

2177 

187 

125 

DS 

1 

6/12/62 

6 

2172 

200 

55 

DS 

m 

6/27/62 

24 

2192 

187 

264 

DS 

5393 

433 

83 

DS 

9307 

500 

90 

DS 

9350 

467 

90 

DS 

n 

7/11/62 

12 

2185 

200 

180 

DS 

Table  II.  Observed  amplitudes  and  periods.  This  data  corresponds  to 
empirical  microbarographs  published  by  Donn  and  Shaw  [1967]  and  by 
Harkrider  [1964],  refered  to  as  DS  and  H,  respectively,  (a-i,  U.S.S.R. 
explosions;  j-n,  U.S.  explosions). 
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The  records  used  include  all  of  those  in  these  two  data  collections  for 
which  Bath  [1962]  gives  a  source  yield  estimate,  for  which  the  pressure 
scale  could  be  determined  and  for  which  the  first  cycle  of  the  signal  is 
recognizable  above  the  background  noise.  Pressure  amplitudes  for  Harkrider's 
records  were  computed  using  his  microbarograph  response  data.  Pressure 
amplitudes  for  the  Donn  and  Shaw  records  were  determined  according  to  the 
premises  (W.  Donn,  private  communication)  that  (a)  all  records  recorded 
by  Lamont  type  A  microbarographs  are  to  the  same  scale  and  (b)  the  clip 
to  clip  amplitude  of  off  scale  oscillations  was  360  pbars.  The  ordinate 
in  Fig.  16  gives  fsin-^' ^ (r/r^)}  in  pbar  MT-^  where  Y  is  the 

explosion  yield  in  MT.  The  abscissa  gives  the  period  j  in  s  •  Note 

that  the  plot  is  full  logarithmic.  The  solid  line  represents  the 
theoretical  relation  of  Eqs.  (4.29)  and  (4.30). 

Although  Posey  and  Pierce  [l972]  have  demonstrated  that  the  value 
of  the  constant  B’  in  Eq .  (4.29)  is  somewhat  sensitive  to  the  assumed 
form  of  the  negative  phase  of  the  blast  wave  at  close-in  distances,  it 
ippears  that  the  value  actually  derived  gives  the  optimum  fit  to  the  data. 

The  scatter  about  the  theoretical  curve  could  be  due  to  various  causes 
other  than  interference  with  the  Lamb  mode  by  other  modes  and  by  noise. 

One  possibility  which  seems  especially  likely  is  the  variation  in  amplitude 
due  to  the  horizontal  refraction  and  subsequent  focusing  and  defocusing 
caused  by  departures  of  the  atmosphere  from  perfect  stratification. 

One  may  properly  argue  that  the  plot  in  Fig.  16  at  best  only  demonstrates 
that  estimates  of  yield  based  or.  the  theory  outlined  here  are  consistent 
with  Bath's  estimates.  Thus  it  may  be  of  some  interest  here  to  briefly 
review  the  basis  for  Bath's  estimates. 

According  to  Bath,  the  yields  of  the  U.S.S.R.  atmospheric  tests  of 
October  23  and  October  30,  1961,  were  announced  as  25  and  58  MT,  respectively, 
and  Bath  says  that  American  estimates  set  the  height  of  burst  for  each  at 
about  3600  meters.  By  ignoring  height  of  burst  effects  and  assuming  a 
linear  relation  between  yield  and  amplitude  of  the  seismic  surface  wave  as 
recorded  at  any  given  station  in  Sweden  for  all  air  blasts  above  the 
Novaya  Zemlya  test  site,  Bath  estimates  the  yields  of  the  remaining 
atmospheric  explosions  in  the  fall,  1961,  Soviet  test  series.  Each  estimate 
is  an  average  based  on  from  3  to  11  observations.  Ba-.h  feels  that  these 
estimates  (shown  in  his  Table  XII)  are  "fairly  reliable".  He  also  gives 
estimates  for  the  yields  of  some  of  the  atmospht  ric  tests  conducted  by  the 
U.S.  near  Christmas  Island  In  1962,  apparently  based  on  seismic  energy 
estimates  from  seismic  surface  waves  recorded  at  Uppsala,  Sweden,  and 
assuming  a  constant  seismic  coupling  factor  a  [seismic  energy/yield]  of 

2X10  \  One  might  tend  to  believe  that  Bath's  yield  estimates  for  the 
Novaya  Zemlya  shots  are  probably  more  reliable  than  for  the  Pacific  shots 
because  results  from  several  observations  were  averaged  to  get  the  former 
and  because  the  latter  were  much  further  from  Sweden.  One  may  note  also 
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Figure  16.  Comparison  of  data  with  the  theoretical  relation  between 
amplitude  and  period  of  infrasonic  waveforms  generated  by  nuclear  explosions. 
The  data  points  are  lettered  a  to  n  corresponding  to  particular  events 
defined  in  the  text. 
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that  Bath's  estimates  are  based  on  the  following  general  relation 
betweer  the  yield  and  the  Richter  scale  magnitude  M: 

M  -  5.1  +  0.7  log  a  Y 

where  a  is  the  seismic  coupling  factor. 


4.5.7  The  anomalous  spatial  variations  in  amplitude  may  be  explained 
in  terms  of  focusing  or  defocusing  of  horizontal  ray  paths 

-1/2 

This  prediction  is  based  on  the  presence  of  the  factor  (db/d0) 
in  the  denominator  of  Eq.  (4.11)  where  (db/dQ)d9  is  the  ray  channel  width 
separating  rays  9  and  0  +  d0.  Such  anomalous  amplitude  variations  have  been 
reported  in  some  detail  for  pressure  waves  recorded  following  the  Soviet 
explosion  of  1961  October  30,  by  Wexler  and  Haas  [1962].  For  example,  the 
amplitudes  recorded  in  Illinois  tended  to  be  2  to  4  times  larger  than’those 
re~oi^ded  in  Nebraska  or  Kansas.  The  ray  equations  suggest  that, 
if  |vL|  <<  c  ,  the  ray  paths  are  nearly  perpendicular  to  lines  connecting 
stations  which  receive  the  pressure  wave  simultaneously.  Thus  the  model 
would  suggest  that  one  might  be  able  to  predict  regions  where  anomalously 
high  or  low  amplitude  waveforms  are  recorded  from  an  examination  of  lines 
of  constant  arrival  time.  This  is  in  fair  accord  with  the  data  exhibited 
by  Wexler  &  Hass. 

Wexlei  and  Hass  [1962]  give  lines  of  constant  arrival  times  over  the 
entire  globe  and  lines  of  constant  amplitude  over  the  northern  hemisphere 
for  the  wave  generated  by  the  Soviet  explosion  of  October  30,  1961,  at 
Novaya  Zemlya.  Records  from  about  1000  barographs  were  used  In  determining 
these  lines.  They  also  show  a  plot  of  average  amplitude  vs.  colatitude 
from  the  source.  For  the  most  part,  these  were  meteorological  barographs 
designed  to  record  much  slower  pressure  variations  than  the  1  to  10  minute 
periods  characteristic  of  the  strongest  portion  of  the  explosively  generated 
acoustic-gravity  waves.  Wexler  and  Hass  also  state  that  most  of  the  charts 
had  a  recording  speed  of  only  9  hour/ in.  This  suggests  that  they  might 
be  relatively  insensitive  to  the  details  of  oscillations  with  periods  of 
15  minutes  or  less  such  as  are  of  interest  here  and  that  one  could  introduce 
significant  time  errors  through  inaccurate  mounting  of  the  paper  (15  min 
is  only  1/36  inch  on  the  chart).  While  the  above  considerations  suggest 
that  one  not  place  too  great  a  reliance  on  the  quantitative  results  presented, 
some  valuable  qualitative  information  on  large  scale  anomalies  in  waveform 
arrival  time  and  amplitude  is  certainly  present.  (One  may  note  that  the 
greatest  rippling  in  the  lines  of  constant  arrival  occur  in  the  U.S.  and 
Western  Europe,  where  the  density  of  reporting  stations  is  large.) 

For  the  reasons  discussed  above,  we  decided  for  the  present  that 
detailed  or  elaborate  computations  of  Lamb  mode  propagation  parameters  from 


-74- 


existing  atmospheric  data  were  not  warranted  by  the  accuracy  of  Wexler 
and  Hass'  arrival  time  and  amplitude  daca  (  a  portion  of  which  is 
reproduced  in  Figs.  17  and  18  here).  Instead,  it  seemed  appropriate 
to  conduct  a  simple  qualitative  comparison  between  the  Wexler  and  Hass 
data  and  the  portions  of  the  sea  level  and  500  mb  (about  6  km)  northern 
hemisphere  weather  maps  for  1200  GMT,  October  30,  1961,  which  are 
reproduced  in  Figs.  19  and  20,  respectively.  That  is,  by  an  examination 
of  the  weather  maps,  assuming  that  the  average  is  about  the  same  for 

rays  which  are  never  separated  by  more  than  a  few  hundred  kilometers,  and 
taking  winds  to  be  geostrophic  (parallel  to  the  isobars,  clockwise  around 
a  high  pressure  area  and  counterclockwise  around  a  low  in  the  northern 
hemisphere,  and  monotonically  increasing  in  strength  with  the  pressure 
gradient) ,  one  should  be  able  to  judge  qualitatively  how  c£,  which  is 
approximately  +  VLk’  var:*-es  whether  or  not  the  undulations  in  the 

lines  of  constant  arrival  and/or  the  amplitude  anomalies  reported  by  Wexler 
and  Hass  are  due  to  this  variation.  Of  course,  implicit  in  such  a  comparison 
is  the  assumption,  which  is  not  necessarily  true,  that  both  the  earliest 
observation  ar.d  the  wave's  largest  amplitude  are  due  to  the  Lamb  edge  mode. 

The  former  is  more  likely  true  than  the  second. 

For  the  sake  of  simplicity,  it  is  assume^  that  the  general  atmospheric 

circulation  pattern,  i.e.,  the  variation  of  v  ,  is  qualitatively  represented 

L 

by  the  winds  at  500  mh;  In  Fig.  21,  the  major  centers  of  low  ui  high  ambient 
pressure  indicated  in  Fig.  20  are  shown  and  numbered  (L^  through  L^  for  the 

lows,  and  H2  and  for  the  highs).  At  each  location  the  direction  of 

the  geostrophic  wind  is  indicated  by  the  arrows.  The  location  of  the  58  MT 
blast  at  Novaya  Zemlya  (73.82on,  53.57<>E  according  to  Donn  and  Shaw  [1962]) 

Is  shown  by  a  star.  As  a  ray  propagates  through  the  center  of  a  low  pressure 
system,  the  rays  to  its  right  should  propagate  more  rapidly  and  those  to 
its  left  more  slowly  because  of  the  counterclockwise  circulation.  The  opposite 
is  true  at  the  center  of  a  high. 

An  unrefracted  ray  (great  circle  path  from  source)  through  low  pressure 
center  L,  would  pass  through  Alaska,  so  that  the  Lamb  wave  should  be  retarded 
in  its  arrival  along  the  northern  coast  of  Alaska  and  western  Canada  relative 
to  its  arrival  just  to  the  west.  The  presence  of  the  high  H.  indicates  a 
retarded  c£  over  the  northern  Canadian  islands  relative  to  the  western 
coast  of  Greenland  and  an  increased  c  down  the  eastern  coast.  All  of  these 
contentions  are  supported  by  the  undulations  of  the  1200  GMT  arrival  line 
shown  in  Fig.  17. 

The  low  pressure  center  L^  is  just  about  on  the  1300  GMT  line,  so 
that  its  effect  should  be  evident  in  the  changes  from  1200  to  1400  GMT,  and, 
indeed,  while  the  former  is  fairly  straight  north  of  ,  the  latter  displays 
a  wiggle  consistent  witli  increasing  c^.  from  east  to  west  just  below  L^. 
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Fig.  17.  Observed  Arrival  Tine  of  an  Infrasonic  Wave.  The  source  of 
the  wave  is  the  Soviet  thermonuclear  blast  over  Novaya  Zemlya,  October  30 
1961,  at  about  0830  GMT,  and  the  recording  instruments  are  barographs  at 
weather  stations  or  on  ships.  The  dashed  rectangle  indicates  the  portion 
of  the  northern  hemispheric  map  shown  in  Figs.  19-21  .  (Extracted  from 
paper  by  Waxier  and  Hass  [1962].} 


Fig.  18.  Maximum  Observed  Amplitude  of  an  Infrasonic  Wave.  This 
map  is  for  the  same  disturbance  as  In  Fig.  17,  and  is  derived  from  the 
same  records.  The  amplitude  is  given  in  mbars.  (Extracted  from  paper 
by  Wexler  and  Hass  [1962].) 


Similarly,  the  effect  of  L0  is  seen  to  he  consistent  with  expectation 
from  a  comparison  of  tne  1200  and  1400  GMT  lines  in  its  vicinity.  The 
effect  of  L~  can  be  seen  by  comparing  the  1400  and  1600  GMT  lines  between 
about  150°E  and  170°W.  Changes  in  c  due  to  winds  about  H_,  Lc ,  H.  and 

are  not  as  clear,  but,  as  can  be  seen  in  Fig.  20,  each  of  these  except 

L,.  is  a  weak  pressure  center  surrounded  by  relatively  small  pressure 
gradients  and  winds.  In  the  portion  of  the  northern  hemisphere  not  shown 
in  the  500  mbar  map  of  Fig.  20,  there  were  only  three  strong  pressure  systems, 
of  which  two,  a  low  just  east  of  Iceland  ard  a  high  over  the  Caspian  Sea, 
seem  to  produce  undulations  in  the  time  lines  consistent  with  the  assumed 
variation  of  v^. 

The  above  comparison  shows  that  many  of  the  irregularities  of  the 
order  of  1000  km  in  the  lines  o£  constant  arrival  time  can  be  qualitatively 
attributed  to  the  variation  of  v  .  Some  other  large  scale  irregularities 

L 

cannot  be  explained  by  the  assumed  circulation  pattern  and  are  probably  due 
to  strong  winds  above  the  500  mb  level  which  do  not  conform  to  that  pattern 
or  to  variations  in  c.  ,  which  are  ignored.  A  portion  of  the  smaller  scale 
irregularities  could  Be  due  to  instrumentation  errors. 

While  the  earliest  portion  of  an  explosively  generated  acoustic-gravity 
wavetrain  is  most  likely  due  to  the  Lamb  mode,  the  largest  amplitude  may 
or  may  not  be  due  solely  to  the  Lamb  mode,  depending  primarily  upon  ducting 
conditions  above  the  observer.  Thus,  it  would  not  be  suprising  should  the 
amplitude  anomalies  of  Fig.  18  not  be  as  easily  explained  by  the  Lamb  mode 
theory  as  are  the  arrival  time  irregularities. 

In  order  to  determine  whether  the  observed  large  scale  variations  in 
overpressure  amplitude  (of  the  order  of  1000  km)  are  the  result  of 
horizontal  refraction  due  to  variations  in  cE  ,  the  empirical  lines  of 
constant  arrival  are  assumed  to  be  lines  of  constant  phase  for  the  Lamb 
mode.  That  is,  rays  should  be  nearly  perpendicular  to  these  lines,  so 
that  a  concave  line  would  lead  to  focusing  of  the  Lamb  mode  rays  while 
a  convex  line  would  lead  to  defocuoing.  This  effect  is  witnessed  in  some 
instances,  such  as  in  the  Atlantic,  where  the  convex  lines  of  constant 
phase  in  the  east  are  accompanied  by  small  observed  amplitudes  and  the 
concave  lines  in  the  middle  Atlantic  are  accompanied  by  larger  amplitudes. 
However,  some  of  the  pronounced  amplitude  variations,  such  as  the  increase 
over  Indochina,  cannot  be  explained  by  the  assumed  Lamb  mode  refraction 
pattern. 
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Fig.  19.  Sea  Level  Synoptic  Weather  Map  for  the  Northern  Pacific  Ocean 
and  North  America.  The  portion  of  the  northern  hemispheric  map  indicated 
by  the  dashed  box  in  Fig.  17  is  shown  here, 
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Fig,  20.  500  robar  Synoptic  Weather  Map  for  the  Northern  Pacific  Ocean 

and  North  America.  The  portion  of  the  northern  hemispheric  map  indicated 
by  the  dashed  box  in  Fig.  17  is  shown  here. 
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CHAPTER  V 


A  SHORT  PERIOD  MODLL  FOR 
WAVES  GENERATED  BY  EXPLOSIONS 


5.1  INTRODUCTION 

The  theories  and  computational  techniques  described  In  the  previous 
two  chapters  apply  best  for  interpretation  of  the  lower  frequency  (2  to 
15  minute  periods)  recorded  at  very  large  distance  from  megaton  class 
nuclear  explosions.  We  have  recently  started  to  develop  an  alternative 
theory  applicable  to  shorter  period  signals  generated  by  lower  yield 
explosions  or  observed  at  closer-in  distances.  In  the  present  chapter, 
the.  salient  features  of  this  theory  are  briefly  reviewed. 


5.2  RAY  PATHS 


The  basic  mathematical  model  presumes  that  the  pressure  wave 
disturbance  at  any  given  point  somewhat  removed  from  the  immediate  vicinity 
of  the  point  of  detonation  can  be  expressed  in  the  form 


n  fv  .  t- 1 
1  »—  »  -  / 


T 

rays 


ft-  *  * 

i  ->.i) 


where  the  sum  includes  one  term  for  each  of  the  geometric  acoustic  ray 
paths  which  connect  the  point  of  detonation  with  the  observer  location, 
including  those  rays  which  may  have  been  reflected  one  or  more  times  at 
the  ground.  (In  order  to  facilitate  bookkeeping,  it  is  assumed  that  the 
obse.  -ation  point  is  never  identically  on  the  ground,  although  it  may  be 
a  very  small  distance  above  it.) 

In  the  simpler  version  of  the  theory,  presently  under  study,  it  Is 
presumed  that  these  ray  paths  may  he  taken  to  be  those  appropriate  for 
higher  frequency  low  amplitude  disturbances.  In  the  limit,  the  paths  are 
independent  of  frequency  and  wave  amplitude.  This  approximation  does  not 
necessarily  imply  a  neglect  of  nonlinear  effects  or  of  frequency  dispersion 
and  has  some  precedence  in  the  contemporary  theories  of  sonic  boom  propagation 
['fliitham,  1956;  Hayes,  Haefei{,and  Kulsrud  1969]  where  nonlinear  affects 
are  neglected  In  the  determination  of  ray  paths  hut  included  In  the 
consideration  of  propagation  along  rr.y  tubes. 

The  method  of  determination  of  such  ray  paths  for  a  temperature  and 
wind  stratified  atmosphere  ha3  previously  been  described  by  cne  of  the 
authors  in  some  detail  [Tierce,  1966]  and  a  computer  program  based  on  that 
nethod  is  presently  in  use  by  W.  Ponn  and  coworkers  at  Lament -Doherty 


Preceding  page  blank 
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Geological  Observatory,  (see,  for  example,  Donn  and  Rind  [1971]).  A  sample 
plot  of  such  ray  paths  Is  shown  In  Kig.  22.  The  appropriate  ray  tracing 
equations  on  which  such  plots  are  based  are 


-*■ 

d* 

dx 


e  +  v 

K 


(5.2a) 


“  -kVc-(k-V)v  -  K  x(Vxv) 
at 


(5.2b) 


where,  for  variable  t,  x(t)  traces  ou^  the  ray  path.  Here  the  ambient 
sound  spee<J  c  and  the  wind^. velocity  v  are  considered  as  functions  of 
position  x.  The  quantity  *  is  an  auxiliary  quantity  termed  the  wave 
slowness  vector  (units  of  reciprocal  velc  ity)  which  also  is  a  function 
of  the  ray  time  parameter  T  .  At  all  points  if  is  related  to  *  by  the 
equation 


(1  -  k-v)2  -  c2x2  (5.3) 

-I* 

In  Eq.  (5.2a)  ,  e  is  the  unit  vector  in  the  direction  of  <.  Note 
that  in  the  absence  of  winds,  k  -  1/c  and  is  the  ray  direction. 

Initial  conditions  on  the  integration  of  Eq.  (5.2)  are  thar  x  at  T  »  P 
represent  the  point  of  detonation  and  that  the  ray  direction  point  in  some 
given  direction  specified  by  parameters  (Oo,$o)which  give  the  polar  and 

azimuthal  angles  which  dx/dx  initially  makes  in  a  coordinate  system 
with  tlie  z  axis  passing  vertically  through  the  point  of  detonation.  Thus 
and  form  a  means  of  labeling  all  rays  which  emanate  from  the  source. 

When  a  ray  meets  the  ground  it  is  assumed  that  it  reflects  specularly 
such  that  incident  and  reflected  rays  locally  near  the  ground  make  the  same 
angle  with  the  ground.  The  horizontal  components  of  k  do  not  change  on 
reflection  but  the  vertical  component  changes  sign.  It  is  of  course  assumed 
that  there  is  no  vertical  component  of  wind  velocity  at  the  ground. 

In  the  computation  scheme  we  envision,  one  of  the  initial  steps  would 
be  to  find  all  initial  ray  parameter  pairs  (Aq,$o)  which  correspond  to  rays 

which  are  governed  by  the  equations  described  above  and  which  connect  the 
source  and  observer  positions.  Each  such  ray  would  then  correspond  to  a  term 
in  Eq.  (5.1). 
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Figure  22.  Computer  generated  plots  of  ray  paths  from  a  point  source  to 
the  east  in  summer.  (Extracted  from  paper  by  Donn  ar i  Rind  1 197 1 ] . ) 
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5.3  PROPAGATION  ALONG  A  RAY  SEGMENT 


For  simplicity  in  the  description  of  the  overall  computation  scheme 
we  speak  of  a  ray  segment  as  a  portion  of  a  ray  path  which  contains  no 
ground  reflection  and  along  which  the  ray  does  not  intersect  any  of  it3 
neighboring  rays  (i.e.,,  the  ray  does  not  touch  a  caustic  or,  equivalently, 
the  ray  tube  area  does  not  vanish).  While  the  existence  of  caustics  (i.e., 
the  surfaces  in  space  composed  of  points  at  which  such  intersections  occur) 
must  be  taken  into  account  in  any  theory  seeking  to  describe  waveforms 
at  all  but  very  close-in  distances)  the  consideration  of  these  is 
postponed  until  Sec.  5.5. 

Along  a  given  segment,  the  pressure  contribution  due  to  that  ray  is 
taken  to  be  of  the  following  form 


Pray  “  J(*)iHt,*) 


(5.4) 


where  a  is  a  parameter  characterizing  distance  along  a  ray  path,  specifically 
related  to  distance  f  along  the  path  by  a  relation  analogous  to  that  of 
Eq.  (4.20)  ,  i.e. 


dl 

ds 


c  +  v.c 

K 

T__ 

fee  +  v' 

K 


(5.5) 


such  that  o  -  t  for  a  windless  atmosphere.  Tne  parameter  s  is  defined 
such  that  it  is  zero  at  the  point  of  detonation.  Thus  any  given  ray  segment 
will  in  general  correspond  to  a  finite  range  of  s  . 


The  parameter  J(s)  in  Eq.  (5.4)  is  independent  of  time  and  defined 
to  vary  with  s  in  such  a  manner  that  wave  action  (Bretherton  and  Garrett,  1968] 
or,  equivalently,  the  Blokhintzev  invariant  [Blokhlntze* ,  1946  ]  is  conserved. 
Thus 


_d 

ds 


J2|d$/dr|A _ )  Q 

PqC2(1  “  K  .  v)J 


(5.6) 


where  pq  is  ambient  density,  |dx/dT|  i3  the  ray  speed  of  Eq.fo.2a) 

and  A  is  ray  tube  area.  Here  A  is  defined  as  the  (s  variable) 
cross-sectional  area  of  any  given  (small)  fixed  >be  of  rays  which  includes 
the  ray  in  question.  Note  that  the  quantity  in  braces  in  Eq.  (5.6)  is 
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independent  of  s  .  Thus,  if  J  is  known  at  any  given  value  of  ,  it 
can  be  found  easily  for  any  other  value  of  S. 


In  the  simplest  limit  of  geometrical  acoustics  with  nonlinear  effects 
neglected,  with  the  neglect  of  dissipation  and  with  the  neglect  of  any 
dispersion  induced  by  the  effects  of  gravity,  the  remaining  factor  ij’  in 
Eq.  (5.4)  satisfies  the  one  dimensional  wave  equation 


dip 

3t 


+  (c  +  v 


3J 

as 


o 


(5.7) 


which  has  the  general  solution 


ip  =  f  (t  - 


r 


ds  j 

c 

eq 


(5.8) 


where  f(t)  is  any  arbitrary  function  and  where  we  have  abbreviated 


c 

eq 


c  +  v-e 


(5.9) 


as  the  effective  wave  speed.  However,  the  approximation  represented 
by  the  above  equations  may  not  be  entirely  sufficient,  so  we  have  sought 
a  generalization  of  the  Eq.  (5.7)  which  incorporates  other  effects. 


5.3.1  Dispersion  Effects 

To  incorporate  dispersion  into  Eq.  (5.7)  we  recall  that  acoustic-gravity 
waves  in  an  isothermal  atmosphere  with  constant  winds  JPierce,  1966] 
satisfy  a  dispersion  relation  of  the  general  form 

kz  "  ^  *  “A)/c2  -  ~  wB)/n2J(kx  +  kv)  (5'10) 

z  a  b  x  y 

where  0  -  w  -  k.v  is  the  Doppler  shifted  angular  frequency;  k  ,  k  ,  k 

x  y  z 

are  the  components  of  the  wave  propagation  vector.  Here  u>^  and  u>g  are  the 
we’l  known  characteristic  atmospheric  frequencies  given  by 


-  (y/2)g/c  (5.11a) 

uig  -  (y-l)l^2g/c  (5.11b) 
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if  this  above  dispersion  relation  is  solved  for  to  ,  one  finds  to 
lowest  order  in  1/k  that  u  ■  c  k  »  corresponding  to  nondispersive 

propagation.  If  the  next  nonzero  term  is  included,  one  finds 


,  .  1  2, 
«  ■  C  k  + 

eq  /  u 


(5.12) 


where 


2  2  //.  2  ^  „2W.,2  .  ,2  ^  .2.,  2 

“  w.  -  {(k  +  k  /  k  +  k  )}  w 

0  A  x  y  x  y  z  B 


(5.13) 


is  a  positive  quantify  which  depends  on  the  direction  of  the  wave 
propagation  vector  k,  but  not  on  its  magnitude. 

Suppose  now  that  i p(£,t)  in  Eq.  (5.4)  is  specified  as  a  function  of  t  at 
some  point  s  on  a  ray  segment  .  Then  if  ip  propagates  according  to  the 
dispersion  reJation  (5.12)  ,  its  time  variation  at  any  other  point  s  on 
the  same  ray  segment  should  be  given  by 


<P(t,s)  “Re  f  }  (w)e-1^“T  + 

'  ft 


where 


(5.14) 


T  -  t 


(5.15) 


rs  2 

I  w 
-2 

Js  c 


0  eq 


(5.16) 


and  where  the  Fourier  transform  *  feat)  1b*  gusli  that 


*(t,«0)  ■  Fe  *  (•d«”i<a£3« 

'o 


(5.17) 
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Note  that  the  resulting  expression  for  iHt.s)  represented  by  Eq.  (5.14) 
would  reduce  to  Eq.  (5,8)  were  it  not  for  the  factor  exp [-i  (1/2)  Q/u]  in  the 
Fourier  transform.  It  should  be  understood,  of  course,  that  (w)  is 

considered  to  be  largest  for  higher  frequencies  than  u>  .  Consistent  with 
this,  one  may  rewrite  (5.14)  in  the  form 


*(t,s)  -  R 


*  n 


(<u)e 


-iwt 


i  t 


prop 


(io2  - 


2  ,1/2 

“eff* 


du  (5.18) 


where 


Cprop 


(5.19a) 


u  „  “  Q/t 

eff  prop 


(5.19b) 


may  be  interpreted  as  the  propagation  tine  and  an  averaged  value  of  u 
along  the  ray  segment,  respectively. 

One  may  next  reexpress  Eq.  (5.18)  as  a  convolution  integral  in  Cue  form 


<»(t 


-  f 


G(t  -  t  ,  9  ,S  )D»(t  ,£  )dt 
o  o  o  >  o 


(5.20) 


where  the  Green's  function  is  given  by 
G(t  -  to,s,so)  *  ^  | 


2  2  a/2 

-iu>(t-t  )  it  (w  -  u>  ) 

o  prop  eff  . 

e  e  (5.21) 


Here  the  integration  contour  extends  slightly  above  the  real  axis  and 

the  branch  cut  extends  from  -w  ,,  to  +  u>  ,,  along  the  real  axis.  The 

eff  eff 

integral  with  some  effort  may  be  evaluated  as 


»  ¥ 


X  H(t-t  -t  ) 
o  prop 


(5.22) 
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'l 

vhere  J^(D)  is  the  Bessel  function  of 
abbreviated  D  as  the  positive  square 


first  order  and  where  we  have 
root  of 


D2  - 


eff 


?  2 

[(t-t  )  -  t  ] 

o  prop 


(5.23) 


The  symbols  <5  and  H  in  the  above  formulas  represent  the  familiar 
Dirac  delta  function  and  Heaviside  step  function,  respectively. 


With  the  substitution  of  Eq.  (5.22)  into  Eq . 

t-t 


(5.20)  one  finds  that 


the  dispersed  waveform  is  given  by 

2 


prop 


<Ht,s) 


ij/  (t-t  ,s  ) 
prop  o 


Oi  r,t 

eff  prop 


(J,  (D)/D)iKt  ,s  )dt 
1  o  o  o 


(5.24) 


One  may  also  show  that  the  quantity  D  may  be  further  approximated  as 


“eff  (2 


t  ) 1/2  (t-t  -t  ) 
prop  prop  o 


1/2 


(5.25) 


for  all  practical  purposes  if  ij/  is  composed  primarily  of  higher 
frequencies  than  • 

Although  we  have  at  present  no  computations  of  Eq.  (5.24)  to 

exhibit,  it  would  appear  (see  rig.  23)  that  what  is  being  described  by  the 

mathematics  in  the  case  of,  say,  a  N  shaped  wave  is  the  development  of 

an  oscillatory  tall  whose  amplitude  and  length  increase  with  propagation 

distance  and  whose  period  at  any  given  point  tends  to  increase  with  time 

at  any  given  point.  The  latter  follows  since  the  zeros  of  J. (D)  occur 

at  successively  larger  Intervals  of  t  for  fixed  t  and  f  .  Also, 

}  °  prop  o 

it  follows  from  the  fact  that  the  group  velocity  according  to  the  dispersion 

relation  of  Eq .  (5.12)  tends  to  decrease  with  increasing  period  (decreasing  u>). 


5.3.2  Nonlinear  Effects 

If  Eq.  (5.7)  were  to  he  modified  to  incorporate  nonlinear  effects,  a 
simple  approach  (which  we  in  essence  adopt  here  )  would  be  for  one  to  utilize 
Whitham's  rule  [Whithara,  1952,  1956]  and  to  replace  the  travel  speed 
•+  ■+ 

c  +  v.e^  by  that  characteristic  of  the  perturbed  rather  than  the  ambient 
medium,  such  that  Eq.  (5.7)  becomes 

<S?>  +  cnl  <ir>  ’  0  <5-J6) 
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Fig.  23.  Sketch  of  the  effects  of  gravitationally  Induced  dispersion  on 
an  originally  N- shaped  waveform.  Successive  plots  correspond  to  later 
values  of  ties.  The  parameter  s  is  equivalent  to  distance  along  a  ray  path. 
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Fig.  24.  Sketch  of  the  application  of  the  equal  area  rale  to  construct 
shock  locations  froet  a  suit  1-valued  vavefora. 


with  the  corrected  wave  speed 


cNL  -  c  +  v.c^  +  c  [  (Y  U)  /  (2Y  )  lp/pQ 


(5.27) 


The  third  term  represents  the  first  order  correction  to  the  sum  of  the 
sound  speed  and  the  fluid  velocity  for  the  case  of  a  plane  wave  propagating 
In  a  homogeneous  ideal  gas.  Here  p  is  the  acoustic  pressure,  in 

Kq.  (5.4)  ,  and  pq  is  the  ambient  pressure.  Equation  (5.26)  is  sometimes 
referred  to  as  the  inviscid  Burger's  equation. 

The  above  derived  nonlinear  equation  for  y  has  the  approximate 
solution  which  can  be  given  in  a  parametric  form  as 


t(s,t) 


(5.28) 


1  "  *+  ‘prop  “  <5,29) 


where 

s 

u  -  [  (V H)/ (2  Y)  ]  I  (J/p  )(l/c  >ds  (5.30) 

J  °  eq 

4 

0 


Here  $  is  an  auxiliary  parameter  defined  such  that  it  equals  time 

t  when  s=  s  .  The  function  <p  (t)  gives  i/<(s,t)  when  $ 

o  o  o 

Thus,  were  s  specified,  one  could  let  <i>  run  through  a  sequence  of 

values  and  list  corresponding  values  of  and  t  computed  from  the 

Eqs.  (5.28)  and  (5.29).  If  these  values  were  plotted  to  give  p  versus 

t,  one  would  have  a  plot  of  ip(s,t)  for  fixed  s  and  variable  t  .  (The 

complication  of  shocks  is  considered  below.) 

In  the  event  that  the  quantity  ip  constructed  for  given  s  becomes 
multiple  valued.it  is  postulated  that  discontinuities  in  the  actual  waveform 
(i.e.,  shocks)  form.  In  the  time  interval  between  shocks,  p  versus  t  is 
still  described  parametrically  by  the  above  equations,  only  with  a  restricted 
range  of  4>.  The  approximate  waveform  which  incorporates  weak  shocks  may 
be  found  graphically  by  first  plotting  the  multivalued  function  described 
parametrically  by  ((>),*  (#)J  for  variable  <p  and  then  placing  shocks 
(vertical  lines  on  a  ip  versus  t  graph)  in  such  a  way  that  the  equal 
area  rule  is  satisfied.  Each  such  shock  hits  the  graph  at  three  points 
and  delimits  an  upper  loop  on  its  left  and  a  lower  loop  on  its  right  (wc 
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assume  time  runs  to  the  right).  The  location  of  the  shock  is  such  thaL 
the  area  of  the  upper  loop  just  equals  that  of  the  lower  loop  (see  Fig.  24) 
and  the  shock  replaces  the  two  loops  in  the  approximate  waveform. 

One  may  show,  albeit  with  some  effort,  that  t’’e  scheme  described  above 
is  consistent  with  the  Kankine-llugoniot  relations  [see,  fot  example,  Landau 
and  Llfshitz,  1959)  in  the  limit  of  sufficiently  weak  shocks.  Also,  if  the 
above  method  is  used  twice  in  succession;  waveform  at  determined  from 

waveform  at  which  was  previously  found  from  waveform  at  sq;  then  one 

can  show  that  the  resulting  waveform  is  just  the  same  as  that  which  would  be 

found  if  the  waveform  at  s_  were  determined  directly  from  that  at  s  „ 

2  o 


5.3.3  Composite  Numerical  Method 

In  order  to  take  nonlinear  effects  and  dispersion  effects  into  account 
simultaneously,  one  may  adopt  the  following  method.  Let  a  given  ray  segment 
be  divided  into  a  small  number  of  segments,  the  length  of  which  may  be  of 
the  order  of  100  km  or  less.  The  optimal  length  remains  a  topic  for  analytical 
investigation,  although  we  doubt  that  one  need  take  it  less  than  25  km.  Let 
the  distance  parameter  at  the  near  end  of  one  of  these  segments  be  s  ,  that 
at  the  further  end  be  s^.  Given  ^(t.s^)  as  a  function  of  time,  one  seeks 

to  determine  <i(t,s  )  as  a  function  of  time. 

TliU  o  F  hr>  <  C  f  /n  doF  niMn  d  n  o  a  f  tin  F  -I  rtn  tL  f  *-  Q  A  fK  F  Vin  ttn.nl  a«F  n  f 
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nonlinear  effects,  but  with  the  consideration  of  dispersion.  This 
intermediate  waveform  is  found  from  Eqa,  (5.24)  and  (5.25),  as  is  discussed 
in  Sec,  5.3.1,  Next  one  defines 


^(t.S  )  -  *  (t  +  t  ,s  )  (5.31) 

o  ProP  1 

as  that  waveform  which  would  have  resulted  in  i^(t,s^)  were  the  propagation 

from  sq  to  to  take  place  without  dispersion.  This  dispersed  but  time 

shifted  waveform  is  then  used  to  define  <|i0(4>)  in  Eq.  (5.28).  The  nonlinear 
propagation  equations  (5.28)  and  (5,29),  along  with  the  method  for  constructing 
shocks  described  in  Sec.  5.3.2,  are  then  used  to  find  the  function  ^>(t,s). 

Alternately,  one  could  consider  nonlinear  effects  first,  then  time 
shift  back,  and  construct  the  dispersed  waveform  from  the  nonlinearly  distorted 
waveform.  Both  methods  should  give  virtually  the  same  result  if  the  distance 
from  sq  to  is  sufficiently  short.  What  is  assumed  in  that  over  short 

distances  the  two  forms  of  distort  n,  nonlinear  and  dispersive,  are  additive. 
Over  large  distances  the  distortion  due  to  either  may  be  large,  but  the 
distortion  is  an  accumulative  effect  which  is  small  over  short  distances. 
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The  above  scheme  can  be  used  repeatedly  over  a  sequence  of  subsegments 
such  that  one  can  in  principle  with,  hopefully,  only  a  modest  amount  of 
computing  find  the  waveform  at  tin  further  end  of  a  chafn  of  subsegments 
from  a  knowledge  of  the  waveform  at  the  nearer  end. 

Similar  considerations  may  be  employed  to  incorporate  viscous  and 
thermal  dissipation  into  the  calculation,  although  we  have  not  yet  derived 
the  appropriate  waveform  distortion  equations. 


5. A  INITIAL  CONDITIONS 


We  here  consider  that  ray  segment  of  a  given  ray  which  begins  at 

the  point  of  detonation.  It  should  be  clear  from  the  general  formulation 

given  in  the  preceding  sections  that  one  needs  to  know  the  waveform  factor 

ij>(t,s)  at  some  relatively  close-in  value  of  s  In  order  to  find  the 

waveform  at  all  larger  values  of  s  .  This  should  in  principle  be  accomplished 

by  matching  Eq .  (5.4)  to  the  blast  wave  at  some  relatively  close-ln  distance 

yet  sufficiently  far  out  that  the  weak  shock  theory  should  be  applicable. 

Since  this  leaves  considerable  lattitude  in  the  starting  value  9  of  s  1 

s  t  ar  ti 

the  initial  conditions  arc  not  completely  clear.  The.  situation  is  further 
complicated  by  the  fact  that,  if  the  explosion  is  moderately  strong  or  if 
it  takes  place  at  higher  altitudes,  then  one  may  expect  the  effects  of  the 
atmosphere's  density  stratification  to  be  appreciable  before  the  shock  has 
become  a  weak  shock  (say,  less  than  10X  ambient  pressure).  This  would 


1  mr*  1  < 
- r 
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or  conform  to  the  idealized  case  of  a  blast  wave  propagating  with  spherical 
symmetry  in  a  homogeneous  atmosphere. 


In  order  to  circumvent  such  difficulties,  we  sought  to  find  just  what 
pioperties  of  the  blast  wave  were  at  least  asymptotically  independent  of 
distance  a  along  a  ray  path.  In  this  regard. one  may  note  that,  among  other 
quantities,  the  blast  wave  passing  a  given  point  may  be  characterized  by  an 
overpressure  P  and  a  positive  phase  derivation  T  .  Both  P  and  T 

8  S  S  S 

vary  with  s,  P  decreasing  and  T  increasing  with  distance  from  the 

8  S 

explosion,  but  two  aggregate  quantities,  even  with  nonlinear  effects  taken 
into  account,  tend  to  remain  invariant.  These  are 


I 


1 


P  T  /J 
s  s 


l2  -  TgJ/Ps  -[<Y+1)/(2Y)] 


{alJ/(ceqPo),d8 


(5.32a) 

(5.32b) 


where  J  is  the  quantity  satisfying  Eq.  (5.6)  and  where  a  Is  any  convenient 
nonzero  constant.  Although  J  is  arbitrary  to  the  extent  of  a  constant 
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multiplicative  factor,  the  choice  of  that  factor  will  not  affect  the 
invariance  of  1^  and  I2  with  a.  (The  proof  of  the  invariance  of  the 

two  quantities  1^  and  Is  based  on  the  Idealization  that,  in  the 

positive  phase  of  the  blast  wave,  the  plot  of  pressure  versus  time  at 
a  fixed  point  is  a  straight  line  which  goes  from  Fg  to  0  in  a  time 

interval  T^.  For  brevity,  the  proof  is  omitted  here.) 


It  is  possibi  to  choose  a  such  that  1^  is  identically  zero.  If  this 

is  done  for  some  value  of  s ,  then  will  be  zero  for  all  a  for  the 

same  choice  of  a  ,  Let  us  denote  this  choice  of  a  big  a*.  Then  a*  may 
also  be  regarded  as  an  invariant. 


Once  a*  and  I,  are  known,  then  the  positive  phase  durvation  T 
l  8 

and  the  shock  overpressure  at  any  value  of  s  (assumed  sufficiently  large 

that  the  weak  shock  theory  is  applicable)  may  be  found  from  the  relations 

a 

Ts  -  {[(Y+1)/(2Y))I1  |  [j/(ce<iPo)]da)1/2  (5.33a) 

a* 


P 

s 

J 


<jf  f  . 


(5.33b) 


Thus  a*  may  be  considered  as  that  point  on  the  ray  where  the  weak  shock 
theory  formally  predicts  Tg  to  be  zero  and  Pg  to  be  infinite.  Even  though 

the  weak  shock  theory  is  clearly  not  applicable  for  s  near  a  a,  we  may 
atill  regard  a*  as  a  convenient  parameter  characterizing  the  explosion. 

The  values  of  a*  and  for  a  given  nuclear  explosion  may  be  estimated 

from  Brode's  [1954]  computations  of  the  blast  wave  generated  by  a  point 
energy  source  in  a  homogeneous  atmosphere.  From  examination  of  his 
figures  which  show  the  blast  wave  at  larger  distances  from  the  source  we 


estimate  (taking  Y  -  1.4  and  requiring  that 
small  values  of  s) 

J  ■  1/s  at  sufficiently 

■*  •  °-7  <v.)1/3 

(5.34a) 

h  '  »-o  Vp/'3  p0/c 

(5.34b) 
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Here  is  the  total  hydrodynamic  energy  (taken  to  be  half  the  total 

energy  E  released)  and  p  and  c  are  the  ambient  pressure  and  sound 
speed  at  the  source. 

The  above  assumes  the  explosion  to  be  some  distance  above  the  ground. 
For  an  explosion  very  close  to  the  ground,  it  would  appear  to  be  more 
accurate  to  replace  by  2E^  in  the  above  and  to  only  consider  rays 

which  propagate  initially  horizontally  or  obliquely  upwards. 

In  so  far  as  the  source  model  is  concerned,  our  choice  for  the  present 

is  that:  (1'  J  as  a  function  of  s  he  computed  from  Eq.  (5.6)  subject  to 

the  requirement  that  J  reduce  to  1/s  in  the  limit  of  small  g,  an 

initial  at  s*s  ■  5a*  be  taken  of  the  general  form  of  a  "Classtone 

,  „  start 


-At/T 

*"<“>  <5-35> 


where  P  and  T  are  found  from  bqs.  (5.33)  and  where  At  is 
ss 

t-t  .  Here  t  is  the  time  that  the  blast  wave  first  arrives 

stare  start 

at  the  point  sstart  on  ray  in  question.  A  suitable  approximate 
expression  foi  t'nis  quantity  would  appear  to  be 


tscart 


2a*  Ceq 


|8start 

2a* 


(v  +1) 
(2>  ) 


J 

c 

eq 


(5.36) 


where  is  given  as  a  function  of  s  by  Eq.  (5.33b).  Note  that  the 

time  origin  (t  ■  0)  corresponds  to  the  time  the  blast  first  readies  the 
point  2a*  .  It  is  tacitly  assumed  that  departures  of  the  blast  from 
spherical  symmetry  are  negligible  before  that  time. 


5.5  PROPAGATION  PAST  CAUSTICS  AND  GROUND  REFLECTION 

Let  us  suppose  that  a  given  ray  under  consideration  touches  a  caustic 
(see  Fig.  25)  at  a  point  where  c«wt. .  The  computational  procedure  described 

in  the  previous  sections  is  assumed  to  have  been  carried  to  the  point  where 
at  some  value  of  »  slightly  less  that  Sc>  say  a  -e  ,  one  knows  that  the 

two  factors  J(s  -e)  and  <J'(t,3  -e)  which,  according  to  F.q.  (5.4),  describe 
c  c 

the  waveform  shortly  before  it  reaches  the  caustic.  The  task  which  now 
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CAUSTIC 


Fig.  25.  Sketch  of  rays  in  the  vicinity  of  t  caustic  (which  fonu 
the  locus  of  points  of  Intersection  of  adjacent  rays).  Here,,  for 
simplicity,  the  radius  of  curvature  of  the  caustic  is  considered  to  be 
substantially  lass  than  chat  of  the  rays,  although  this  is  not  always  the 


presents  itself  is  the  determination  of  J  and  tfi  at.  some  point,  say 
8  +c ,  on  the  ray  past  the  caustic. 


As  is  well  known  (see,  for  example,  Maglieri,  Hilton,  Huckel ,  Henderson, 
and  McLeod,  1970J,  waveforms  on  passing  caustics  undergo  a  considerable 
alteration  in  shape.  This  alteration  takes  place  within  a  relatively 
short  distance  of  several  wavelengths;  and  subsequently  the  waveform 
propagates  along  the  ray  with  relatively  little  additional  distortion. 

Thus  we  can  adopt  the  idealization  that  the  alteration  is  abrupt  and 
the  choice  of  c  becomes  largely  immaterial. 


For  simplicity,  we  ignore  nonlinear  effects  in  the  vicinity  of  the 
caustic  and  take  the  alteration  of  waveform  to  be  as  predicted  by  linear 
acoustics.  The  linear  acoustic  theory  which  includes  diffraction  effects 
of  waveform  alteration  at  a  caustic  is  now  relatively  well  known  and  has 
been  discussed  in  the  literature  by,  among  others,  Friedlander  [1958], 
Tolstoy  [1965,  1968],  Ludwig  [1966]  and  by  Sachs  and  SiJbiger  [1971]. 

The  general  prediction  of  the  linear  theory  as  applied  to  the  problem 
at  hand  is  that  J  at  Sc+E  may  be  found  by  requiring  that  the  quantity 

in  braces  in  Eq.  (5.6)  is  the  same  at  s  +c  and  s  -c  .  The  choice  of 

c  c 

J  is  understood  to  be  a  positive  quantity.  Then  the  function  if(t,  s^+r) 


is  predicted  to  be  the  Hilbert  transform  of 


il>(t--t  ,s  -e)  such  that 
prop  c 


.I(t,sc+E) 


fit-i  . B-r; 

o  prop  c 


dt 


(5.37) 


where 


Cprop 


s  +e 


c 


c 


(5.38) 


is  tlie  small  increment  of  time  required  for  a  wave  to  propagate  from 

s  -e  to  s  +c  .  In  Eq.  (5.37), the  capital  P  implies  that  the  principal  val 
c  c 

should  be  taken,  i.e.  ,  for  any  function  f(t) 


-iff  1  dt 

t  -t  O  X  ]  £ 

n  '  n 


(5.39) 


ue 


-OB 


Note  that,  as  long  as 
side  remains  bounded. 


f(t)  is  continuous,  the  integrand  on  the  right  hand 
even  in  the  limit  of  small  £  . 
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ijt 


In  order  to  give  the  reader  some  additional  insight  into  the  waveform 
alteration  represented  by  a  Hilbert  transform,  vie  show  in  Fig.  26  some 
plots  of  functionsand  their  transforms  which  wc:  extracted  from  the  paper 
by  Sachs  and  Silbiger  [1971], 


One  of  the  complications  which  must  be  considered  in  reference  to 
the  above  described  theory  is  that,  if  tire  incident  waveform  lias  any 
discontinuities  (corresponding  to  weak  shocks),  then  its  Hilbert  transform 
will  have  a  logarithmic  singularity.  11  ■«.•  existence  of  such  singularities  should  be  evident  in  tK' 
plots  in  Fig- 26.  In  particular ,  if  the  function  f(t)  in  Eq .  (5.39)  has  a 
positive  jump  upwards  of  A  at  a  value  t.  of  t,  then  the  Hilbert  transform 


(i.e.,  the  integral)  will  be  dominated  at  times  near 


t . 
3 


by  a  term  of  the  form 


log  sing. 


7 


(5.40) 


While  the  presence  of  such  singularities  would  add  a  certain  unreclism 
to  the  theory,  a  formal  application  of  the  nonlinear  propagation  theory 
described  in  Sec.  (5.32)  using  the  equal  area  rule  to  locate  and  define 
shocks  would  indicate  that  each  such  singularity  is  immediately  transformed 
into  a  finite  strength  shock  as  it  propagates  past  the  caustic. 


A  given  logarithmic  singularity  described  by,  say,  Eq.  (5.40)  after 
propagating  some  distance  according  to  Eqs.  (5.28)  and  (5.29)  will  be 
described  locally  by  the  two  lines 


t 


t 


prop 


±e~^’/& 


-H<l> 


(5.41) 


obtained  by  eliminating  (f  from  Eqs.  (5.28)  and  (5.29).  These  two  lines 

describe  a  drooped  singularity  such  as  sketched  in  Fig.  27a.  The  multivalued 

configuration  is  reduced  to  a  finite  amplitude  single  valued  waveform  by 

placing  a  snock  at  t  which  intersects  the  unmodified  waveform  in  Fig.  27a 

at  ,  and  >p  ,  such  as  is  sketched.  The  equal  area  rule  requires  areas 

a  b  c 

Aj  and  sketched  in  the  figure  to  equal  area  A^.  One  may  easil'" 

check  that,  even  though  the  unmodified  waveform  is  singular,  the  area  A^  is 
finite  and,  moreover,  given  by 


A 


1 


2A  -Trip  /  A 
—  e  a 
x 


(5.42) 


Once  one  works  out  expressions  for  the  other  areas,  he  finds  the  equal 
area  rule  implies  that 

*  -  *  «  (1.20) (2A/x)  (5.43) 

a  c 


i 


* 
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F(  t  )=»,  0"  t-M0  | 
F(t)=0,  t«O.I>-i0 


FH(t)4loq|-^- 


Figure  26.  Three  simple  pulse  shapes  and  their  Hilbert  transforms. 
(Extracted  from  paper  by  Sachs  and  Silbiger  [1971].) 
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Fig.  27.  Sketch  of  hov  nonlinear  effects  reduce  a  lcgarithmic 
singularity  to  a  shock  of  finite  amplitude.  (a)  Unmodified  multivalued 
waveform,  illustrating  application  of  equal  area  rule.  (b)  Modified 
waveform  vlth  constructed  shock. 
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where  the  factor  1.20  represent r  the  root  of  the  trancendental  equation 


x  “  coth  x 


(5.44) 


Thus,  in  so  far  as  the  singularity  dominates  the  Hilbert  transformed 
waveform,  the  net  discontinuity  is  independent  of  r  and  hence  of 
propagation  distance.  The  peak  value  ij'  (  of  iji  depends,  however,  on 
li,  being  given  by 


A 

« 


,  ..20it, 

tn  ( — — ) 
pA 


(5.4b' 


which  reduces  rapidly  as  n  increases. 

If  we  use  the  prescription  outlined  above,  we  may  take  the  waveform 
which  is  received  just  past  the  caustic  to  he  finite  and  we  may  therefore 
analyze  its  subsequent  propagation  on  the  ray  segment  beyond  the  caustic 
according  to  the  method  outlined  in  Sec.  5.3. 

The  only  other  case  which  requires  special  consideration  is  reflection 
of  a  wave  at  the  ground.  Consistent  with  the  other  approximations  we  have 
made  it  would  appear  that  we  should  take  J  and  if(s,t)  to  he  the  same  for 
both  incident  and  reflected  rays  on  the  ground. 


5 . 6  REMARKS 


The  general  method  for  finding  a  given  ray's  conti ihution  to  the 


overall  pressure  waveform  is  to  take  p 


ray 


in  the  form  Jifi  as  in  Eq .  (5.4) 


where  the  values  of  J(s)  and  <{i(t,s)  at  a  relatively  close-in  value  of 
3  arc  found  according  io  the  procedure  outlined  in  Sec.  5.4.  The  distortion 
of  the  waveform  propagating  along  the  ray  is  then  computed,  ray  segment 
by  ray  segment,  with  suitable  alterations  at  caustics  and  upon  ground 
reflection,  according  to  the  method  outlined  in  Sec.B.  5,3  and  5.5  .  The 
total  waveform  is  then  found  by  a  summation  of  all  such  ray  waveforms  which 
correspond  to  ray  paths  which  connect  the  point  of  detonation  witli  the 
observation  point. 


The  method  lias  a  number  of  limitations  which  should  he  pointed  out. 
First  we  are  tacitly  ignoring  all  nonlinear  interactions  between  rays  which 
pass  through  the  same  point.  In  particular,  nonlinear  interaction  between 
incident  and  ground  reflected  waves  is  ignored.  Also,  we  have  not  taken 
into  account  any  nonlinear  contributions  to  ray  refraction  since  ray  paths 
are  taken  from  linear  acoustics.  Diffraction,  which,  among  other  things, 
is  largely  responsible  for  the  excitation  of  the  Lamb  edge  mode  discussed 
in  Chapter  4,  is  handled  only  in  a  gross  way  with  its  effects  considered 


only  when  wt  tukc  into  account  the  distortion  aL  a  caustic.  The  general 
method  breaks  down  when  the  observation  point  is  very  near  a  caustic  since 
the  nonlinear  effects  at  a  caustic  have  been  neglected.  Also,  the  method 
of  incorporating  frequency  dispersion  into  the  theory  is  such  that  one 
cannot  have  any  confidence  in  the  lower  frequency  contribution  to  the 
overall  waveform,  especially  lot  frequencies  much  less  than  a  representative 
lower  atmosphere  (corresponding  to  a  period  of  the  order  of  4  minutes). 

Nevertheless,  if  one  specifically  wishes  to  consider  the  higher  frequency 
portion  oi  the  waveform  corresponding  to  periods  less  than,  say,  two  minutes, 
the  method  outlined  "hove  could  very  possibly  give  theoretical  waveforms 
which  have  more  than  token  -esemblancc  to  actual  waveforms  recorded  in  the 
field.  The  overall  method  would  appear  tc  he  somewhat  more  sophisticated 
than  any  existing  computational  scheme  based  primarily  on  geometrical  acoustic 
concepts  in  that  it  considers  frequency  dispersion  and  distortion  at  caustics 
in  addition  to  nonlinear  effects.  It  does  seem  feasible,  however,  to  develop 
a  general  computer  program  encompassing  the  method  which  could  complement 
tiie  program  INFRASONIC  WAVhi'OKMS  described  in  Chapter  3  in  regards  to  the 
domain  of  application. 
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chapter  VI 


A  Horn  EOK  ACOUSTI C-ERAV  ITV 
EXCITATION  BY  1IKEBALL  RISE 


t>.l  INTRODUCTION 

The  models  of  acoustic-gravity  vnvr-  excitation  By  nuclear  explosions 
discussed  in  t  ho  preceding  chapter;;  ol  this  report  are  based  on  the  general 
tenet  that  the  near  field  blast  wave  is  the  source  of  acoustic  gravity  waves 
observed  at  large  distances.  Tin  authors  believe,  however,  that  the  processes 
associated  with  the  rise  of  the  fireball  subsequent  to  the  general  Ion  of 
the  blast  wave  may  be  an  important  source  of  waves  in  some  instances.  These 
instances  include  the  later  arriving  lowei  frequency  waves  which  might  be 
observed  dirocLly  below  the  source  or  at  high  altitudes.  Fireball  rise  is 
probably  not  a  significant  sodreo  of  the  earlier  arriving  waves  observed  at 
(’round  level  at  large  horizontal  distance;;,  however. 

Wo  may  note  that  fireball  rise  was  previously  considered  and  subsequently 
dismissed  as  a  significant  source  of  far  fiel  radiation  by  Hunt,  Ta Inter , 
and  Penney  [i960].  We  believe,  however,  that  these  authors’  arguments  in 
support  of  this  tenet  are  fallacious  since  they  are  based  on  Lhc  assumption 
that  the  wave  amplitudes  generated  by  fireball  rise  should  fall  off  with 
increasing  distance  faster  than  t  hose  generated  by  the  blast  wave  and  since 
there  is  no  upriui  i  icason  to  believe  the  fireball  rise  generated  waves  to 
bo  loss  susceptible  .0  channeling  or  guiding  by  the  inhomogeneous  atmosphere 
than  those  arising  from  the  blast  wave.  They  would,  however,  arrive  later. 

More  recently,  the  question  of  fireball  rise  as  a  possible  source  mechanism 
was  reopened  by  Tolstoy  and  Eau  11971 J .  These  authors  argued  that  t lie  rise 
mechanism  may  well  be  the  dominant  source  of  wave  generation,  especially  at 
larger  periods,  if  the  fireball  is  of  substantially  large  radius.  The  present 
authors  found  Tolstoy  and  Eau's  excellent  discussion  of  this  topic  to  be 
most  intriguing  and.  consequently  sought  to  refine  the  Tolstoy-I.au  model  to 
Include  a  more  realistic  and  less  phenomenological  source  mechanism.  The 
source  model  used  by  Tolstoy  and  Eau  consisted,  among  other  features,  of  a 
point  vertical  force  term  on  the  rigtit  hand  side  of  the  Euler's  equation 
(i.c.,  Newton's  second  law  for  a  fluid)  in  its  linearized  form.  The  point 
of  application  of  the  force  was  presumed  to  be  rising  nr  constant  speed 
through  the  atmosphere.  The  magnitude  of  the  force  was  inferred  from  Warren's 
| 19(>0]  calculations  of  the  wave  drag  on  a  rigid  sphere  rising  at  constant 
speed  in  a  density  stratified  medium.  Thus  the  force  depended,  among  other 
quantities,  on  the  radius  of  the  sphere,  the  density  of  the  fluid,  and  the 
velocity.  The  radius  of  the  sphere  and  its  velocity  were  estimated  independent ly 
from  other  known  features  of  a  nuclear  explosion. 
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Iii  the  present  chapter,  we  propose  a  somewhat  illfferent  model  which 
takes  into  account  (lie  entrainment  of  air  hy  the  rising  fireball  and  which 
incorporates  a  dynamical  model  oi  iirehall  rise  into  the  for  null  at  ion.  We 
do  not  explicitly  use  Warren's  [19<>0]  wave  ding  expression  as  it  is  difficult 
to  ascertain  the  extent  to  which  a  rising  turbulent  1  trehall  may  he  simulated 
hy  a  rip, id  sphere.  An  earlier  version  oi  the  theory  given  here  is  reported 
in  t'  ■  paper  hy  Pierce  (1972J. 

It  should  he  noted  that  a  somewhat  similai  model  oi  wave  generation 
hy  fireball  rise  is  presented  in  the  recent  paper  hy  Murphy  and  Kalinins  [1972]. 
[Unfortunately,  at.  the  time  of  this  writing  we  hav  not  yet  had  sufficient 
time  to  make  a  detailed  comparison  of  the  two  models.] 


(>.  2  DYNAMICAL  MODKI.  OK  K1KKBALI  RISK 

According  to  the  computations  of  Brodo  [1955],  the  blast  wave  propagnt ing 
out  from  the  point  o f  detonation  carries  away  a  finite  amount  of  mass,  such 
that,  in  the  absence  of  pravitational  effects,  there  remains  a  spherically 
shaped  rep,ion  of  neplip,ible  density.  The  radius  of  (Iris  repion  appears 
[Pierce,  10(>8]  to  he  o*  the  general  order  of  map.nituJe  of 

Ro  “  °-/,(EH/lo)1/3 


where  is  the  hydrodynamic  energy  release  (taken  ns  one  half  of  tlie  total 

energy  released)  and  p  is  the  ambient  pressure  at  the  source  locatiorr. 


Wc  niay  expect  thnt  such  a  bubble  will  rise  under  the  influence  of  gravity 
in  such  a  manner  that  it  retains  its  identity,  if  not  its  spherical  shape. 
Taylor  [1950]  argues  on  the  basis  of  potential  How  in  the  air  directly 
above  the  bubble  and  assuming  constant  pressure  within  the  bubble  thnt  the 
init  tal  rate  of  rise  should  be 


v 

init. 


(2/3) 


(fi  K  ) 

o 


i  n 


(6.2) 


.Scorer  [1957],  on  the  basis  of  similitude  considerations  and  from  experimental 
data  on  the  motion  of  globular  masses  of  different  density  under  the  influence 
of  buoyancy  in  a  incompressible  fluid,  finds  that  the  velocity  of  "bubble" 
rise  may  be  represented  by  the  semi-empirical  equation 

1.2(r[(po-  Pi)/fic]R)in 


V  - 


(6.3) 


wlidi’  K  Is  n  1 enpt  li  ana  lupous  to  radLus  which  may  bo  inlet  |>i  ct  ed  as  t  lie 
taditis  of  a  sphere  of  the  same  volutin-  anil  whole  o^  is  the  ileus  ily  of  the 

substance  Inside  the  "bubble".  This  apices  with  Taylor's  formula  except 
fin  the  magnitude,  of  the  coefficient.  Since  some  del  Inite  choice  should  be 
made  at  this  point,  we  choose  to  use  l!q.  (ti.3)  in  the  discussion  that 
f ol 1 ows . 

One  may  note  that  hq.  (0.3)  may  he  interpreted  as  the  velocity  ol  a 
bubble  rising  with  negligible  acceleration  such  that  the  net  buoyancy  force 
Oxer ted  on  it  is  counterbalanced  by  a  drag  force,  i.e. 

r„  “  h(n0- 


The  relation  of  this  drag  force  to  1  lie  velocity  may  be  inferred  from  an 


elimination  of 


i 


from  i'.qs.  (0,3)  and  (0,4)  with  the  net  result  that 


(4/3)fl 

(1.2)2 


„2  2 
K  p  v 
o 


(c.5) 


One  may  note  that  this  diap  force  is  proportional  to  the  square  of 
the  velocity,  to  the  apparent  surface  men  of  the  bubble,  and  to  the  density 
of  the  surrouiullnp,  medium.  Apart  from  the  numerical  coefficient,  Lhis 
formula  mipht  also  he  inferred  from  similitude  considerations. 


L  t  OO  l-t'll 


iniCr  uwi'y  cOtl.jt'eiuiC  L  lUil  <i 


suitable  form  of  Newton's  sc  com!  law  for  the  motion  of  the  fireball  niiftlit  be 


y|—  {  (M  +  K  M  )  (H,.  -  M)p  -  F„ 

dt  disp  dt  disp  -1  1) 


(fi.6) 


where  the  sipn  of  F  is  opposite  to  that  of  the  velocity  and  where  It  is 
the  mass  of  pas  within  the  fireball,  M(Jl  is  the  mass  displaced, 


M..  “  (4/3) ii  IT  P 

disp  o 


((>.7) 


and  K  has  a  magnitude  such  that  K  M  represents  the  effective  Inertia 

di-'P 

of  the  surrounding  air.  Alternately,  one  may  interpret  K  to  he  of  the 
magnitude  such  that 
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A  K  nJUpw«/jl)2 


(ki;) 


out.  side 


«>.H) 


equals  tlie  kinetic  energy  oi  the  ail  exleilor  to  the  fireball.  One  possible 
choice  for  K  is  the  value  of  1/2  which  would  he  appl  icahle  to  the  case  where 
tl  ie  fiieball  is  a  rigid  sphere  end  the  sun  omul  inj’.  medium  is  incomp  t  ossibl  e 
and  homopeueous .  (See,  (or  example,  l.amb  11879,  1943],  p.  i24) .  (hit  present 

opinion  is  that,  foi  the  actual  i  nhomop.eiieuus  and  comptessihle  atmosphere, 

K  may  vary  1  rom  0  to  1/2.  Its  value  durinp,  the  early  stup.es  of  fiiehail  rise 
anJ  growth  would  piohahly  he  more  neatly  1/2,  while  that  at  the  later  stapes 
is  conjectured  to  he  closer  to  zero.  (These  remarks  are,  howevet  ,  larpely 
unsubstantiated  at  present.)  In  order  !o  have  a  definitive  model,  we  take 
K  “  1/2  in  what  follows. 

The  entrainment  of  mass  hy  the  risinp  fireball  may  he  formally  taken 
into  account  by  considerinp  M  as  well  at  z  to  be  a  dependent  variable  in 
Fq .  ((>.6).  An  additional  equation  describing  the  entrainment  of  mass  may 
he  taken  as 


where  the  entrainment  constant  «  is  inferred  from  Scorer's  [1957]  experiments 
to  have  a  value  of  1/4.  Note  that  the  above  assumes  the  entrainment  rate 
to  he  proportional  to  Lhe  area,  the  velocity,  and  the  density  of  the  outside 
medium. 

To  complete  the  dynamical  model,  a  third  equation  is  needed,  since  K 
and  p ^  cannot  he  aptiori  specified  as  functions  of  time.  For  this  purpose, 

a  suitable  approximate  equation  may  he  derived  hy  considerinp  the  net  change 
in  density  p ^  in  some  time  interval  At  as  formally  equivalent  to  a  two 

step  process.  First  R  expands  hy  (AK)^  with  no  net  chanpe  in  density 

distribution  such  that  a  net  mass  increment 

AM  -  (4nR2)p  AR  (6.10) 

o  a 


is  added.  This  mass  then  mixes  irreversibly  with  Lhe  mass  inside  to  give 
a  new  averaged  density 

3  (AR) 

pi  +  APj  -  Pi  +  — -  Pj)  (6.11) 
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Thus  the  coni  1  ibut  Ion  ol  the  entrainment  process  to  the  net  rate  of 
change  of  density  Is 


3"l 

‘dT* 


ent t  alnmeut 


) 


1 

V 


JV 

.It 


(0.12) 


wliei'e  V  Is  the  volume  enclosed.  ThrouRbout  the  piocess  just  described 
It  Is  presumed  that  the  pressure  inside  the  bubble  remains  always  e.ptal 
to  that  ol  the  surround  lap  medium.  When  the  bubble  rises  to  a  layer  of 
different  ambient  pressure,  the  dc.isity  Is  assumed  to  readjust  itself 
ad ial at  leal ly ,  such  that  the  contribution  o'  the  external  pressure  variatfoir 
to  tin-  net  late  of  change  of  internal  density  is 


(JT> 


P  res  sui  t? 


ni  a'’o 


(6.13) 


Thus,  after  add  i  tip,  Kqs.  (<> .  12)  and  ((>,13),  wc  find 


dt  I  f  lit 

O 


+  (0 

O 


av 

.it' 


t !.  -\  f.\ 


or,  with  the  use  of 


0  .V 


M  and  o  V  -  M  and,  after  some  alnebra,  we  find 
o  disp 


_d 

dt 


Unfa  VY)) 
o 


_  y_  dM 
P  V  dt 


(f>.i-r>) 


l.et  ns  next  note  that  f'.qs.  ((>.6),  ((>,‘0,  and  (b.lb)  may  be '  considered  as 
three  simultaneous  ordinary  differential  equations  for  z ,  M,  and  V  (or  R)  as 
iunctlons  of  time,  Suitable  initial  conditions  would  appear  to  be  that  , 
at.  t  “(':(!)  R  equal  the  R  Riven  by  iq .  ((>,1);  (2)  z  equal  the  lu-ip.ht  of 

°  ]  /  ? 

detonation  7  ;  (3)  d//dt  be  Riven  by  1 , 2  (r  K)  and  (‘'i)  the  initial  mass 

M  be  zero.  Once  the  ambient  pressure  and  density  are  specified  as  functions 
of  z,  it  would  appear  to  he  a  st  ralRht  forward  matter  to  nun  rically  intop, rate 
these  equations  on  a  computer  and  thus  fully  describe  M,  V,  K,  or  as 
a  function  of  time  t. 
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(i.3  AlTKOXIMATl:  SOLUTION  ot  KIKLHALL  KIM!  lOUATIONS 

Sonic  remarks  iti.iv  In-  In  older  at  this  point  iMmi-rn  1  up.  tin1  uatuie  * » 1 
the  solution  ol  the  eqwaL  Ions  presented  In  t  lie  pievloits  Sih  t  Inn.  Up  in 
a  l  line  somewhat  before  tin-  velocity  fitsl  l>o  comes  zeto,  It  would  appeal  that 
a  reasonable  appi  ox  l  mat  ion  mlpht  he  to  ncp.lect  the  Inertial  term  In  l.q.  (6.0). 
In  this  event.  it  is  convenient  to  cons  tiler  lather  than  t  as  the 

I  lulepciulent  variable,  such  that  lq.  becomes 


dj_ 

dz 


1 


(1.2  ft  K) 


1/2 


M 

*1oV 


1 


(t..  1(.) 


where 

V  -  (4/3) nR1 


Ot.17) 


I'.quntion  ((>.9)  may  also  he  rewritten  ns 


dll 

dz 


1  (4  Ilk  )p 


while  Kq.  ( (i .  13)  becomes 


_d 

dz 


It  i(pov*  ) 


-3.-  ill 

p  V  dz 


(6.1H) 


(6.19) 


If  one  eliminates  dM/dz.  from  the  lnttct  two  of  the  above  equations, 
he  finds  after  so.ne  alpehra  that 


..  1  /  (3  >  ), 

k  >’o  1 


"  1 


!/(«>) 


which  integrates  to 


R  «  R  in 
o 


-1/(33  ) 


r  p,i,<3y>- 

z 

o 


(6.20) 


(6.21) 
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This  may  in  turn  be  inserted  into  Eq.  (6.18),  yielding  M  as  the  integral 
over  z  from  tc  z  of  the  right  hand  side.  Then  Eq.  6.16  mi.y  be 

integrated  directly.  Although  the  integrations  are  cumbersome,  they  are 
straightforward . 

In  the  case  of  an  isothermal  atmosphere,  Eq.  (6.21)  integrates  to 


R 


R 

o 


+  Y  3  Y  H 


,  CA3  Y  H) 
le 


1J 


(6.22) 


where 


H  »  c2/(Yg) 


(6.23) 


is  the  ecale  height  and  we  have  abbreviated  t  for  z  ~  z0  ‘  Thus,  K 

tends  to  first  increase  linearly  with  z  but,  once  4  becomes  comparable  to 
several  scale  heights  (if  the  sphere  actually  rises  that  far),  the  growth  is 
exponential . 

An  estimate  of  the  peak  height  of  rise  may  be  taken  as  that  value  of  z 
for  which  dt/dz  in  Eq.  (6.16)  becomes  singular,  i.e.  where  M  *  V,  or 

equivalenLly ;  where 

(z  2  3 

3a  p  R  dz  *  p  R  (6.24) 

)  o  o 

z 

o 


In  the  case  of  an  isothermal  atmosphere,  this  requires  that  Rq  and  z 
jointly  satisfy  the  transcendental  equation 

Y  I (B  +  N  -  l)3-  (N  -  l)3]  -  B  (2  -  r)(N3Y  -  1) 


6  Y  B 
3  Y-l 


3  Y 

3  Y  -2 


(N 


3  Y 


N2)  -  0 


(6.25) 


where  we  have  abbreviated 

B  -  (R  /H)/(3  Ya  ) 
o 


<6. 26a) 


(z-z  )/ (3  YH) 

N  *  e 


(6.26b) 
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In  the  nominal  limit  (which  is  not  necessarily  of  greatest  Interest)  of 
vecy  small  explosions  at  lower  altitudes,  where  R  <<  H,  the  solution  of 
Eq.  (6.25)  is  found  to  be  J 

z  -  z  -  y H  B3  (6.27) 

o 


Since  R  is  proportional  to  the  cube  root  of  yield,  this  would  suggest 
that  the  net  rise  is  proportional  to  the  yield  for  smaller  yield  bursts.  We 
at  present  have  no  information  as  to  whether  this  is  consistent  with  what 
has  been  observed. 


The  presence  of  the  inertia  term  in  Eq.  (6.6)  implies  that  the  fireball 

does  not  simply  rise  to  the  equilibrium  altitude,  but,  instead,  overshoots 

and  then  oscillates,  such  as  is  sketched  in  Fig.  28,  with  a  period  of 

roughly  the  Brunt-Vaisala  period  at  the  stabilization  altitude.  The  oscillatory 

phase  of  the  motion  may  be  described  approximately  if  one  formally  considers 

z  -  2.  -  Az  to  be  small,  where  is  the  stabilization  altitude  or, 

n  M 

equivalently,  that  where  M  first  equals  p^V  .  Thus,  we  set  M  «  +  AM 

and  V  ■=  +  AV,  where  and  V  are  the  corresponding  values  of  mass 

and  volume  when  this  altitude  is  reached.  The  "linearired"  versions  of 
Eqs.  (6.9)  and  (6.15)  then  require 

-4  (AM)  -  [4*R2p  «  (Az)  |  (6.28) 

Ol  O  It  UL 


o  dt 


(AV)  - 


dt 


(AM)  + 


V 


dt 


(At) 


(6.29) 


Note  that  the  assumed  nature  of  the  entrainment  process  is  such  that  the  mass 
continues  to  grow.  However  PqAV  -  AM  is  independent  of  entrainment. 

In  Eq.  (6.6),  if  we  discard  terms  of  higher  than  first  order  (with  the 
exception  of  the  drag  force,  which  we  retain),  we  find  that 

a?  x1  +  a£«‘»  +  H„“S 4Z  -  *  fd  <6-30> 


where 


<i> 


2 

B 


(6.31) 


should  be  positive  and 


is  recognized  as  the  square  of  the  Brunt 


-Vaisala 
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ALTITUDE  OF  BUBBLE 


TIME 


Fig.  28.  Sketch  of  the  prediction  of  the  theory  outlined  in  the  text 
In  regards  to  the  height  of  a  rising  fireball  (or  hot  gaseous  bubble)  as 
a  function  of  tiae. 
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frequency  at  the  stabilization  altitude.  It  is  interesting  that  Kq.  (6.30) 
is  independent  of  the  entrainment  coefficient. 

Preparatory  to  solving  Eq.  (6.30),  it  is  convenient  to  rewrite  it  in  the 

form 

2 

(1  +  K  )M  ~  (Az)  +  M  o>2A  z  «  +  M  (R_)-1(d(Az)/dt]2  (6.3?) 

o  o  B  o  V 


wher*»  from  the  expression,  Eq.  (6.5),  for  we  identify  the  constant  as 


%  •  (i-2>2  «« 


(6.33) 


Next,  one  may  note  that  Eq.  (6.3£)  has  the  energy  integral 
{(1  +  K)v2  -  |  Ry2  <^[1  +  2(Az)/RD]}e_2(A2)/RD 


(6.34) 


where  the  upper  or  lower  signs  are  to  be  used  depending  on  whether  the 
"bubble"  is  rising  or  descending.  This  quantity  should  remain  constant  for 
any  phase  of  the  motion  during  which  the  sign  of  v  remains  unchanged. 

During  the  initial  rise  phase,  it  is  clear  that  the  constant  must  be  identically 
zero  since,  otherwise,  the  solution  would  not  natch  on  to  that  represented 
by  Eq ,  (6.16).  Thus,  up  until  v  first  becomes  zero,  we  have 


v 


1/2  “r’S)  1/? 

^  [Kp  -  2(Az)]1/Z 


(6.35) 


For  Az  negative  end  |Az|>>  this  should  reduce  to  Eq.  (6.16). 
Note  that  the  net  amount  of  overshoot  is  Rp/2. 


After  v  first  reverses  sign,  one  finds  in  a  similar  fashion  that 


r  T1^2 

where  we  have  abbreviated 

jf  *>  exp  [  -1  +  2Az/Rjj] 


The  above  expression  for 
rise  recowences)  when 

Az  “  -  .  39 


v  indicates  a  second  turning  point  (  such  that 


*D 


(6.3g) 
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Expressions  describing  subsequent  oscillations  of  the  fireball  can  be 
worked  out  in  analogous  fashion.  The  above  would  indicate  that  the 
amplitude  of  the  oscillation  is  initially  of  the  same  order  of  magnitude 
as  the  fireball  radius.  The  period  would  appear  also  (taking  K  to  be 
between  0  and  1/2)  to  be  of  the.  order,  but  slightly  in  excess  of,  the 
Brunl-Vaisala  period. 


6.4  DERIVATION  OF  SOURCE  TERMS 

The  presence  of  the  rising  and  subsequently  oscillating  fireball  in 
the  atmosphere  can  be  formally  taken  into  account  in  so  far  as  wave 
generation  is  concerned  by  adding  source  terms  to  the  right  hand  sides  of 
the  linearized  equations  of  atmospheric  fluid  dynamics.  For  simplicity, 
in  the  present  discussion  we  neglect  winds,  such  that  these  equations 
become 


3pi 

at"  +  V*  (poui)  "  sm 


po  aSj/at  +  vPl  +  RPi22  -  tv 


(6.39a) 

(6.39b) 


3p,/3t  4  yV  •  (p^)  +  (y  -  l)poguly  *  (y  -  1)S 


o°  lz 


(6.39c) 


where  the  zeroth  order  equations  require 


dpQ/dz  «-gp  .  (6.40) 

Here  p  ,  u^,  and  are  the  acoustic  perturbations  (first  ord^r  quantities) 

to  density,  fluid  velocity,  and  pressure.  The  quantities  S^,  and 

are  source  terms,  whose  detailed  forms  are  discussed  below. 

In  the  paper  by  Tierce  [1972],  a  general  derivation  is  given  for  these 
source  terms  which  may  be  applied  heiv.  Let  r^(t)  be  the  position  of  the 

renter  of  the  rising  fireball  which  is  specified  as  a  function  of  time,  such 
as  would  be  obtained  from  a  solution  of  the  equations  presented  in  the  previous 
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sections.  Then  one  may  formally  consider  each  source  term  to  be  formally 
expandable  as  a  ^um  gf  monopoles,  dipoles,  quadrupoles,  eLc.  ,  where  each 
is  centered  at  x  ■  r  .  Thus,  for  example, 

SM  “  SMo  6(*  “  ?c)  *■  SM1  3 i  6(*  '  V  +  •••  <6‘A1) 


where  the  first  term  is  identified  as  the  mass  monopole  source  and  the 
second  term  as  the  mass  dipole  source.  (Note  that  we  have  here  defined 
our  symbols  somewhat  differently  from  that  of,  the  paper  by  Pierce  (1972].) 
The  reason  that  the  mass  dipole  source  term  involves  only  a  z  derivative 
is  because  of  the  assumed  axial  symmetry  of  the  fireball. 


6.4.1  Monopole  Source  Terms 

As  long  as  we  restrict  ourselves  to  the  consideration  of  acoustic- 
gravity  radiation  with  wavelengths  somewhat  larger  than  the  diameter  of 
the  rising  fireball,  the  importance  of  the  source  terms  should  diminish 
rapidly  with  increasing  order.  Thus,  a  first  approximation  would  be  to 
take  just  the  monopole  terms (l.e. ,  those  containing  a  delta  function  rather 
than  a  derivative  of  a  delta  function)  as  sources.  The  corresponding  coefficients 
for  these  monopole  sources  turn  out  to  be 


SMo  '  -  dT  -  "disp’ 


im  -  -  (d/dL )(t  -  tdlsp)  -  <« -  Hdlsp)8 

SEO  ‘  -  <d/dt)(f '  fdl.p’  -  -  W  ' 


(6.42a) 

(6.42b) 

(6.42c) 


where  M,  x,  aijd  £  are  the  mass,  mora’ntum,  and  energy  of  the  fireball, 
while  ^jisp’  ndisp’  an<*  ^disp  are  t*,e  corresPc,ndfn8  quantities  which  would 


hold  for  a  volume  of  the  same  size  of  the  ambient  atmosphere  which  would 
occupy  the  same  region  of  space  if  the  fireball  were  not  present.  Since  the 


ambient  atmosphere  is  quiescent,  x 


is  zero,  although  M 


is  not. 


The  quantity 


disp  *  - -  disp 

may  be  considered  as  just  the  thermal  e\  rgy  of  the 


disp 

ambient  air.  However,  since  the  fireball  is  assumed  to  have  the  same  pressure 
as  the  ambient  air,  the  thermal  energy  of  the  fireball  per  unit  volume  (not 
per  unit  mass)  will  be  the  same  as  that  of  the  air.  Thus  the  difference 

£  ^disp  should  be  the  kinetic  energy  of  the  gas  inside  the  fireball.  [Note 
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that  in  the  aforementioned  paper  by  Pierce  [ 1972]  ,  there  are  some  obvious 
misprints  and  that  the  correction  of  the  source  terms  to  account  for  the 
presence  of  th"*  ambient  atmosphere  was  erroneously  omitted  in  Eqs.  (18) 
of  that  paper. ] 


^Since  the  fireball  is  rising  vertically,  the  only  nonzero  component 
of  is  the  vertical  component,  which,  with  reference  to  Eq.  (6.6),  is 


just 


S 

iroz 


_<1 

dt 


(KM 


disp 


(6.43) 


Also,  F.q.  (6.9)  allows  us  to  write  the  mass  source  coefficient  in  the  form 

dz 
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Mo 


o  dz  dp 

-•‘larh-v  ^ 


dt 


(6.44) 


The  remaining  source  term  coefficient  S  may  he  estimated  by  taking 
the  kinetic  energy  of  the  gas  inside  the  °fireball  to  be  the  same  as  that 
for  a  Hill  spherical  vertex  lLamb,  1879,  1945]  ,  which  turns  out  to  be  just 
twice  that  for  a  rigid  sphere  of  the  same  mass  rising  at  speed  dz^/dt.  Thus 


M..  (dz  /dt)‘ 
disp  c 


(6.45) 


and  one  accordingly  obtains 
d  ? 

Ssv.  “  -  TT  (  M(dz  /dt)  }  -  g  M  dz  /dt 
to  at  c  c 


(6.46) 


In  order  to  find  these  three  coefficients  S„  ,  S  ,  S_  explicitly 

Mo  hqz’  Eo  r 

as  functions  of  time,  one  must  in  principle  first  find  M,  i^,  R,  and  V 

from  the  solution  of  the  fireball  rise  equations,  then  insert  these  into 
the  equations  derived  above. 

It  may  be  of  some  interest  at  th?s  point  to  compare  our  expression, 

Eq.  6.43,  for  S  with  that  used  in  the  aforementioned  paper  by  Tolstoy 

Tl  OZ 

and  I.au  [1971],  These  authors  in  effect  take 

S  *=  +  K  (6.47) 

it  oz  2 
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•or  a  rising  fireball  where 


Is  approximated  by 


K(Vo 

2  4  ,  ,  ,  ,  .2 

fl  C  (dz  /dt) 

c 


(6.48) 


Here  we  have  made  appropriate  insertions  into  Tolstoy  and  Lau’s  Kqs.  (38) 
and  (71)  to  put  It  into  the  notation  of  the  present  report.  The  above 
expression  should  be  contrasted  with  the  Kq .  (6.5)  p,ivcn  previously  which  wan 
extracted  from  Scorer's  experiments.  Tiro  discrepancy  appears  to  he 
considerable  when  one  considers  the  disparate  velocity  dependences  of  the  drag 
forces  in  the  two  expressions.  There  seems  at  present  to  he  no  clear  cut 
means  of  deciding  which  of  the  two  is  more  nearly  correct,  although  a  force 
which  increases  with  increasing  velocity  would  seem  to  be  more  physically 
plausible. 


6.4.2  Mass  Dipole  Source 

Although  we  do  expect  the  influence  of  higher  order  multipole  source  terms 
to  diminish  rapidly  as  the  multipole  order  increases,  the  possibility  may  be 
considered  that  the  mass  dipole  terra  may  lead  to  significant  wave  generation. 
Here  we  estimate  the  magnitude  of  the  mass  dipole  source  term  following  t^e 
general  formulation  of  Pierce  [1972]. 

The  coefficient  S  which  appears  in  F.q.  (6.41)  for  the  mass  dipole 
coefficient  is  given  hv  an  integral  of  the  form 

SM1  *  "  j  <*M  >  <z  '  zc>  d  A  (6.49) 

where  the  integration  Is  over  an  area  A  enclosing  the  fireball  and  where 
.  given  hy, 

(f»M  *  PQ  Uj  .  n  (6.50) 


represents  the  mass  flowing  out  through  the  hounding  surface  per  unit  time 
per  unit  area  to  first  order.  As  a  rough  approximation  to  this,  we  let 
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and  disregard  the  z  dependence  of  p 


Doing  this  gives 


SK1  ~  "  Po 
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where  M  is  the  mass  idsplaced  by  the  rising  fireball, 

diep 
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0.5  RLMARKS 


The  solution  of  Fqs.  (0.39)  for  the  special  ■  asc  of  an  isothermal 
atmosphere  is  easily  derived  (I, in  and  Yoh,  1972  ;  Pierce,  1972]  and  may 
he  used  to  gain  insight  into  the  generation  of  waves  in  more  realistic 
atmospheres.  For  brevity,  this  solution  is  not  reproduced  here  although 
some  of  its  features  should  he  evident  without  reference  to  the  detailed 
derivat ion . 

One  niay  note  that,  during  the  rise  phase  of  the  fireball  motion,  the  time 
dependence  of  the  various  source  terms  is  not  oscillatory.  The  dominant 
property  of  the  source  terms  (other  than  magnitude)  is  that  their  point  of 
application  is  rising.  Thus  one  would  expect  a  wave  wake  trailing  the 
fireball  which  is  composed  primarily  of  waves  whose  phase  velocity  is  less 
than  the  vertical  source  velocity.  Since  lower  frequency  waves  (depending 
mi  their  frequency)  may  have  phase  velocities  considerably  less  than  the 
sound  speed,  it  is  not  necessary  that  the  source  rise  at  supersonic  speeds 
in  order  to  have  appreciable  wave  generation.  However,  a  subsonic  source  will 
predominantly  excite  waves  of  frequency  less  than  the  Brunt -Vai sal  a  frequency. 

A  second  feature  of  interest  is  the  fact  that,  after  the  rise  phase, 
the  analysis  indicates  that  the  fireball,  and  lienee  the  source  term  point  of 
nppl ication, should  oscillate  with  a  frequency  somewhat  less  than  .  To 
first  order  in  z  ,  the  mass  dipole  source  term  would  appear  to  he  equivalent 
to  a  stationary  source  whose  strength  is  oscillating  in  the  same  manner  as  the 
velocity.  Since  the  governing  equations  are  linear  one  would  expect  a  radiated 
wave  which  oscillates  with  the  same  period  as  the  fireball  itself.  Source 
terms  which  are  of  second  order  in  z  may  he  expected  to  generate  waves  of 
Lwice  the  frequency  or  half  the  period.  Since  waves  with  periods  near  the  Brunt 
period  travel  extremely  slowly  (the  group  velocity  is  nearly  zero)  one  might 
expect  that  at  larger  distances  and  earlier  times  the  dominant  frequency  of  the 
wave  generated  by  the  fireball  oscillation  is  more  nearly  twice  the  Brunt- Vaisnla 
f  requency . 

At  larger  distances,  it  may  he  more  nearly  appropriate  to  solve  Fqs.  ((>,39) 
by  the  guided  mode  method.  In  this  event,  it  would  appear  to  be  a  straightforward 
task  to  extend  the  existing  program  INFKASON1C  WAVF. FORMS  discussed  in  Chapter  3 
of  the  present  report  to  include  rising  and  oscillating  sources,  and  to  include 
the  source  model  introduced  hete  in  particular.  Additional  comments  concerning  the 
implementation  of  the  theory  may  he  found  in  Sec.  7. A. 
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Chapter  Vll 


•’(INCLUSIONS  AND  RECOMMENDATIONS 


7.1  MULTIMODE  SYNTHESIS 

The  computer  program  1NFRASON1C  WAVEFORMS  appears  to  be  a  relatively 
useful  tool  for  the  analysis  ami  Interpretation  of  infrasonic  pressure 
signatures  detected  at  largo  distances  from  nuclear  explosions.  The  program 
is  now  in  use  in  a  number  of  laboratories  and  has  already  given  us  some 
interesting  insights  into  the  dependence  of  waveforms  on  atmospheric  temperature 
and  wind  profiles,  and  on  energy  yield  and  height  of  burst. 

The  limitations  of  the  present  program  ore,  nevertheless,  somewhat 
stringent  and  it  is  to  he  hoped  that  future  work  may  he  devoted  to  removing 
these  limitations.  It  would  appear  entirely  feasible,  for  example,  to  modify 
the  program  such  as  to  include  leaking  modes  in  the  synthesis.  Also,  the 
method  of  integration  over  frequency  might  be  improved  such  that  analytic 
asymptotic  expressions  for  the  Fourier  transforms  are  utilized  to  extend 
the  range  of  integration  over  a  wider  range  of  frequencies,  especially  higher 
frequencies  and  very  low  frequencies. 

One  also  should  he  able  to  extend  the  model  such  that  it  would  be 
applicable  to  atmospheres  where  the  temperature  and  wind-velocity  profiles 
vary  with  horizontal  coordinates. 

The  present  source  model  might  also  be  improved.  In  the  present  version, 
a  distinction  must  he  made  between  near  surface  bursts  and  air  hursts.  It 
would  seem  feasible  to  have  a  somewhat  generalized  model  which  takes  the 
nonlinear  reflection  of  shocks  at  the  ground  into  account  and  which  gives  a 
source  model  whose  parameters  depend  continuously  on  the  height  of  burst. 

The  development  of  a  source  model  which  adequately  describes  higher 
altitude  hursts  appears  to  be  a  relatively  difficult  problem,  and  it  is 
clear  that  additional  efforts  are  required  to  devise  such  a  model.  Perhaps 
the  departures  from  spherical  symmetry  of  shocks  generated  by  high  altitude 
hursts  may  he  taken  into  account  with  the  incorporation  of  dipole  and  higher 
order  multipole  terms  in  the  source  term  expressions  which  we  add  to  the 
linearized  equations  of  atmospheric  fluid  dynamics. 

Finally,  we  may  note  that  the  wide  latitude  in  choice  of  input  parameters 
(number  of  modes,  range  of  frequencies  and  phase  velocities,  etc.)  in  the 
present  version  of  the  program  may  lead  to  some  ambiguities  in  its  numerical 
predictions  of  waveforms.  It  would  appear  feasible  to  modify  the  program 
such  that  the  user  not  be  forced  to  make  too  many  seemingly  arbitrary 
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decisions  regarding  Input  pai ameter s .  Thus,  for  example,  the  computer 
may  decide  Internally,  subject  to  requirements  of  numerical  accuracy,  just 
how  many  and  which  modes  are  required  for  an  adequate  synthesis  and  over 
Just  which  ranges  of  frequencies  the  integration  should  extend.  In  this 
manner,  the  output  of  the  program  will  become  more  standardized. 


7.2  LAMB  KbCi:  M01)£  SYNTHKSIS 

The  hypothesis  that  the  dominant  portion  of  the  early  arriving  lower 
frequency  parts  of  the  waveforms  observed  at  large  distances  is  carried 
by  the  real  atmosphere's  counterpart  of  Lamb's  edge  mode  for  the  isothermal 
atmosphere  seems  to  he  amply  substantiated  by  the  numerical  studies  carried 
out  during  the  contract.  The  hypothesis  leads  to  a  number  of  Intereating 
implications  which  are  mentioned  in  Sec.  A, 5  of  the  present  report.  One 
oi  tile  most  intriguing  of  these  is  that  an  estimate  of  energy  yield  which 
is  relatively  insensitve  to  details  of  atmospheric  structure  may  lie  obtained 
from  the  period  and  amplitude  of  the  earlier  portion  of  the  waveform.  Our 
analysis  indicates  that  such  yield  estimates  agree  fairly  well  with  Bath's 
estimates  based  on  seismic  records. 

The  principal  advantage  of  the  edge  mode  theory  is  its  simplicity. 
Numerical  synthesis  based  on  this  theory  requires  only  a  small  fraction  of 
the  computer  time  required  by  the  multimode  synthesis  program  INI’RASONIC 
WAVKFORMS .  Also,  this  simplicity  makes  the  theory  much  more  readily  amenable  to 
the  consideration  of  horizontal  variations  ir.  atmospheric  structure. 

The  present  edge  mode  synthesis  program,  written  by  Posey  (1971), 
as  mentioned  in  Chapter  A,  strictly  holds  only  for  perfectly  stratified 
atmospheres.  However,  the  incorporation  of  horizontal  variations  into 
the  program  seems  to  be  a  task  of  considerably  smaller  magnitude  than  that 
required  for  the  present  multimode  synthesis  program.  Thus,  one  would 
recommend  that  a  suitable  modification  to  the  Lamb  edge  mode  program  be  made 
within  the  near  future. 

A  possible  improvement  to  the  edge  mode  model  might  be  the  inclusion 
of  the  next  higher  order  term  in  the  dispersion  relation.  Also,  the  question 
of  the  truncation  of  integration  over  height  discussed  in  Sec.  A. 2  might 
be  examined  in  some  more  depth. 

Undoubtedly,  any  theory  based  on  the  propagation  of  just  one  mode,  even 
though  it  may  be  a  composite  mode,  Is  going  to  have  considerable  limitations. 
Thus  it  would  appear  to  be  of  interest  to  investigate  ways  in  which  the  model 
might  be  augmented  to  include  waves  otfier  than  the  edge  mode.  There  are  two 
ways  in  which  this  might  be  done.  Dm  s  simply  to  Identify  one  or  more  other 
guided  modes  which  may  contribute  to  the  waveform  and  to  develop  a  theory 
of  comparable  simplicity  for  their  excitation  and  propagation.  Iht  other  way 
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would  be  to  add  on  n  higher  frequency  wave  which  propagates  out  from  the 
burst  according  to  a  geometrical  acoustics'  theory.  Thus  one  might,  for 
example,  simply  add  the  wave  given  by  the  edge  mode  theory  to  that  which 
would  be  lound  by  the  method  outlined  in  Chapter  V. 


7.3  HIGHER  FREQUENCY  MODEL 

The  theoretical  mode]  discussed  in  Chapter  V  of  the  present  report 
appears  to  have  some  promise  for  the  Interpretation  of  portions  of  wave 
trains  which  have  periods  somewhat  less  than  2  minutes.  The  question  of 
the  range  of  frequencies  whose  propagation  is  adequately  described  by  this 
model  remains  somewhat  open,  hut  it  would  appear  that  the  range  should  extend 
up  to  periods  as  long  as  30  sec,  and  perhaps  as  long  as  2  minutes.  Numerical 
experimentation  would  undoubtedly  help  in  resolving  rhis  question. 

The  model  contains  three  principal  innovations  wuich  should  help  to 
offset  the  limitations  which  one  normally  associates  with  a  ray  theory 
computation.  One  of  these  is  the  incorporation  of  nonlinear  effects  in 
the  propagation  along  ray  paths,  similar  to  what  is  currently  done  in  sonic 
boom  propagation  computations.  Another  is  the  incorporation  of  gravitationally 
induced  dispersion  into  the  propagation.  The  third  innovation  is  the  use 
of  the  Hilbert  transform  to  describe  the  alteration  in  shape  of  the  waveform 
when  it  passes  a  caustic.  [Although  the  use  of  such  a  technique  appears  to 
lie  relatively  well  known  to  workers  in  underwater  acoustics,  it  lias  evidently 
not  yec  been  incorporated  in  a  general  purpose  computer  program  for  arbitrarily 
stratified  oceans.] 

The  method  we  suggest  for  incorporation  of  the  source  into  the  higher 
frequency,  propagation  model  appears  to  be  somewhat  more  satisfactory  t ban 
those  presently  used  for  the  multimode  and  edge  mode  syntheses  in  that  it 
takes  into  account  t he  interaction  of  near  field  nonlinear  effects  with  the 
variation  of  atmospheric  density.  Thus  the  model  may  be  applicable  for  a 
somewhat  higher  range  of  burst  altitudes. 

While  there  are  a  number  of  limitations  to  the  model,  towards  the 
removal  of  which  future  theoretical  efforts  might  be  directed,  our  present 
view  is  that  the  model  is  sufficiently  promising  in  its  present  form  to 
warrant  its  implementation  in  the  form  of  a  digital  computer  program.  Once 
such  a  program  is  written,  it  would  appear  to  be  of  considerable  interest  to 
compare  its  predictions  with  data.  In  this  respect,  one  might  note  that  the 
data  for  such  a  comparison  need  not  be  restricted  to  just  nuclear  explosions, 
but  could  include  high  energy  (HE)  explosion  data  or  data  taken  following 
accidental  natural  gas  explosions. 


7.4  WAVES  EXCITED  BY  FIREBALL  RISE 

The  model  described  in  Chapter  VI  for  the  generation  of  acoustic-gravity 
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waves  by  rising  and  subsequently  oscillating,  fireballs  has  two  principal 
features.  First,  one  has  a  model  for  the  dynamics  of  the  fireball  motion 
which  tacitly  ignores  wove  generation.  Next,  one  has  a  scheme  [Fierce,  1972) 
lor  deriving  effective  source  terms  for  the  linearized  equations  of 
atmospheric  fluid  mechanics  from  the  previously  derived  properties  of  the 
fireball.  The  actual  solution  of  t lie  resulting  linearized  equations  can, 
ut  least  for  larger  horizontal  distances,  he  found  by  using  n  multimode 
synthesis  program  only  slightly  different  from  the  present  1NFKAS0N1C  WAVEFORMS. 

The  model  suggested  for  the  fireball  dynamics  is  probably  still  too 
simplistic  for  actual  events,  although  there  appears  to  be  very  lit.tle 
Information  in  this  regard  pertaining  directly  to  nuclear  explosions  in  the 
open  literature.  Thu  model  suggested  here  has  the  otic  principal  advantage 
that  it  relies  on  a  minimum  of  empirical  parameters.  Given  the  yield  and 
height  of  burst,  the  model  will  give  a  relatively  complete  (albeit,  perhaps 
not  entirely  correct)  description  of  the  fireball  motion. 

It  is  important  that  one  realize  that  the  technique  for  extracting 
source  terms  from  properties  of  fireball  dynamics  as  described  here  and  in 
the  paper  by  Fierce  [1972]  is  not  necessarily  restricted  to  the  particular 
model  of  firebill  dynamics  described  in  Secs.  6.2  and  6.3  of  the  present 
report.  Thus,  should  one  devise  a  more  accurate  or  detailed  model  of  fireball 
rise,  he  still  should  be  able  to  use  the  same  technique. 

As  yet,  we  have  not  carried  the  fireball  rise  wave*  generation  model  to 
the  point  of  making  detailed  numerical  predictions.  While  this  would  appear 
to  be  a  relatively  uncomplicated  task,  it  might  be  appropriate  to  iirst  review 
all  other  sources  of  information  pertinent  to  the  phenomenology  of  fireball 
rise.  In  t he  absence  of  one's  finding  any  more  realistic  model  in  the  near 
future  which  is  not  so  complicated  as  to  be  excessively  difficult  to  use  in 
the  source  term  formulation,  the  implications  of  the  present  model  should 
be  explored,  if  only  to  give  some  clearer  indication  of  the  magnitude  and 
general  qualitative  nature  of  the  contribution  of  fireball  rise  to  the 
far  field  waveforms.  In  particular,  one  might  examine  some  of  the  assertions 
made  in  the  article  by  Tolstoy  and  Lsu  [1972]. 


7.5  OTHER  TOFICS 

The  work  done  under  the  contract  on  the  generation  of  Infrasouud  by 
natural  sources  (such  as  turbulence  and  severe  weather)  is  only  briefly 
described  in  the  present  report.  The  paper  by  Moo  and  Fierce  11972],  the 
abstract  of  which  is  included  in  Chapter  II  oi  the  present  report,  gives 
the  status  oi  tills  research  as  of  April  1972.  Further  work  has  since  been 
done  towards  the  prediction  of  the  frequency  spectra  of  infrasonic  waves 
propagated  at  large  distances  from  localized  tropospheric  disturbances  with 
given  statistical  properties.  Since  this  work  is  still  in  progress  and  will 
shortly  be  described  In  some  depth  in  a  doctoral  thesis  by  C.  A.  Moo,  it  was 
decided  that  a  detailed  description  of  the  work  accomplished  to  date  in  the 
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picscnt  report  would  be  unwarranted.  Some  further  Indication  of  the  work 
done  in  tills  respecL  may  be  found  in  the  following  abstract  of  a  papi r 
currently  In  preparation  lm  presentation  at  the  April  197J  meeting  oi  the 
Acoustical  Society  of  America. 

Nonlinear  Theory  ol  Infrasouic  l-'luctuation  Spectra.  C,  A.  Moo  and 
A.  D.  Pierce,  Massachusetts  Institute  of  Technology,.  Cambridge,  Mass.  02139. 

In  a  previous  paper  (AGAKI)  Conference  Proceedings  No.  115  (1972)]  the  authors 
have  suggested  that  the  2-5  minute  period  oscillations  observed  in  the 
ionosphere  during  periods  of  thunderstorm  activity  In  the  troposphere  may 
be  explained  in  terms  oi  the  interaction  oi  the  velocity  fluctuation  fields 
in  thunderstorm  regions.  This  analysis  is  here  extended  to  the  invest  ip, at  ion 
of  tlie  nonlinear  interactions  of  a  system  oi  propagating  and  evanescent  internal 
acoust ic-pravi ty  waves.  The  development  rests  on  Hamilton's  principle  for 
nonlinear  acoustic-gravity  waves  with  the  Lagrangian  approximated  to  third 
order.  The  coherent  nonlinear  Interaction  problem  is  discussed  with  the 
artifice  of  two  primary  waves  which  interact  in  a  restricted  replon  of  space. 
Although  two  interact  inp  prnvity  waves  cannot  produce  a  third  propagating 
pravity  wave  of  second  order  they  may  produce  an  acoustic  wave  in  the 
atmosphere,  'i  Ire  pcncral  problem  of  internet  iup  random  wave  fields  is 
formulated  and  the  evolution  of  an  initially  monotonic  pravity  wave  spectrum 
due  to  noulineur  interactions  is  discussed.  In  addition  to  affording  a 
possible  explanation  of  the  anomalous  ionospheric  oscillations,  the  nonlinear 
interaction  theory  may  explain  the  veil  known  spectral  pap  which  is  found 
in  the  kinetic  energy  spectrum  of  atmospheric  turbulence  during  periods  of 
convective  activity. 

7 . b  CONCLUDING  RKMARKS 

The  subject  matter  of  the  present  report  has  primarily  been  concerned 
with  the  prediction  ol  acoustic-gravity  waveforms  at  large  distances  following 
a  nuclear  explosion  in  the  atmosphere.  The  emphasis  of  this  research  has 
essentially  been  on  the  development  of  one  or  more  theoretical  models  which 
give  a  sufficiently  detailed  numerical  piedlct’on  to  enable  either  a  qualitative 
or  qualitative  comparison  with  data.  Our  development  of  such  models  has  been 
guided  largely  by  what  has  been  done  by  other  workers  in  the  past,  bv  a 
consideration  of  those  physical  phenomena  which  appeat  most  important  ,  and  by 
a  desire  that  the  model  he  the  simplest  possible  (at  least  in  conceptual  terms) 
to  explain  or  interpret  those  features  ol  the  data  which  appear  pertinent. 

Since  the  atmosphere  is  always  changing,  it  is  always  imperfectly  known 
at  any  given  instant.  Since  the  waveform  one  receives  at  a  given  point 
is  a  functional  of  this  atmosphere,  it  appears  that  detailed  agreement  of  theory 
with  experiment  will  never  be  possible.  However,  what  does  appear  to  he  the 
case  is  that  certain  features  of  the  waveform  may  be  less  sensitive  to  small 
details  of  atmospheric  structure  than  others.  This  became  vividly  apparent 
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to  the  authors  when  they  found  that  the  I.amb  edge  mode  modal  predicted 
the  product  of  period  to  the  three-halves  power  and  amplitude  to 
correlate  better  with  energy  yield  than  either  period  or  amplitude,  separately, 
for  the  case  of  the  early  arriving  low  frequency  signal.  This  type  of 
result  exemplifies  what  we  hope  should  emerge  from  further  studies  of  this 
type.  Computer  programs  such  as  INFRASONIC  WAVEFORMS  may  be  useful  in  thi6 
respect  primarily  because  they  allow  one  to  carry  out  numerical  experimentation, 
from  which  simple  trends  (such  as  yield-amplitude  proportionality)  might  be 
discovered.  Once  such  trends  are  discovered,  they  may  guide  one  in  the 
development  of  simpler  models  which  emphasize  the  essential  physics  of  the 
phenomena  without  unnecessary  complications. 
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